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On the Theory of Cooperative Phenomena in Crystals 
[Continued from Vol. 9, No. 34, p. 244, April 1960] 


By C. Domb 
King’s College, Strand, W.C.2 


§4. APPROXIMATE METHODS 
4.1. Introduction 


The approach which we have used so far in the theory of cooperative 
phenomena has not followed the historical development of the subject. 
Approximate methods had been widely used for several years before any 
exact information on two-dimensional models, which we have discussed in 
the previous chapter, became available. The simplest approximation, 
usually referred to as zeroth order, was used by Bragg and Williams (1934) 
in a discussion of order—disorder transitions in alloys. It introduced for 
the first time the concept of long-range order, and with relatively simple 
mathematics seemed to account for the main features of the experimental 
results. The Curie temperature was characterized by the disappearance 
of long-range order, and was associated with a specific heat singularity. 
The form of the predicted singularity was not in complete accord with 
experiment, the most unsatisfactory feature being the drop of the specific 
heat to zero above the Curie temperature. 

In 1935 Bethe introduced an improved approximation (usually referred 
to as a first-order approximation) which took account of short-range order 
as well as long-range order; the specific heat curve derived from this 
approximation did not fall to zero above the Curie temperature, but a 
small ‘tail’ remained, and this was in better agreement with experiment. 
The approximation was applied by Peierls (1936) to the Ising model, 
and independently Guggenheim (1935) developed an approximation of 
comparable order which he applied to the theory of regular solutions. 

The exact approach to the two-dimensional Ising model by Kramers 
and Wannier in 1941, and the complete solution in the absence of a magnetic 
field by Onsager in 1944 showed clearly the limitations of the above 
approximations in predicting behaviour in the neighbourhood of the Curie 
point. In fact the exact solution depended on an infinity of order para- 
meters, whereas the approximations took account of only one or two of 
them. The form of the specific heat curve in Onsager’s solution differed 
substantially from those obtained by the approximations, but differed 
equally markedly from experimentally observed specific heat curves. 
Perhaps the most striking difference was the large ‘tail’ of the Onsager 
curve, which showed that a substantial fraction of the total entropy 
change of the system took place in the temperature region above the 
Curie point. It was therefore important to determine whether this 
behaviour arose from the two-dimensional nature of the model, and how it 
was modified for a three-dimensional model 


P.M:S. : 
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We have already pointed out that the methods used in the previous 
chapter for deriving exact solutions depend on particular properties of 
two-dimensional nets, and cannot be generalized to three dimensions. 
In the years which elapsed after the appearance of Onsager’s solution, 
despite several claims, no real progress was made towards an exact solution 
of a three-dimensional model; the problem seems to be one of outstanding 
difficulty. However, higher order approximations in closed form were 
developed and systematic methods were devised of obtaining terms of the 
series expansions discussed in § 3.6. 

The literature on such approximations is very extensive, and it has 
not been practicable to give a detailed account of all the different approaches 
which have been used. We have, instead, selected those which seemed to 
us to be the most significant, for detailed consideration, and have en- 
deavoured to summarize the remainder. Reliable information regarding 
the true thermodynamic properties of a given model can best be obtained 
from a series of successive approximations, which enable one to assess 
the general trend of a particular property and hence to estimate the error 
involved in stopping at a particular stage. Unfortunately the labour 
arising in such approximations increases rapidly with the order of the 
approximation, and the process cannot usually be followed for more than 
2 or 3 stages. 

Exact series expansions have been developed quite extensively; when 
they converge might up to the critical region they can usually provide 
useful and reliable information. When they do not converge in the 
critical region some re-arrangement of the terms must be sought, and 
approximations in closed form can be a useful guide to the asymptotic 
behaviour of the series. 

By combining information obtained from series expansions and higher 
order approximations in closed form, we believe that most of the important 
properties of the three-dimensional model can be reliably established, and 
we give an account of these properties in §4.7. We should then be in a 
better position to compare the results with experiment, and to see whether 
a serious refinement of the model is needed to secure effective agreement. 


4.2. Zero Order Approximations 


The first attempt to account in detail for the properties of a ferromagnetic 
was made many years ago by Pierre Weiss (1907) who postulated the 
existence of a large ‘internal field’. Using this hypothesis Weiss was 
able to reproduce the most important physical features of a ferromagnetic, 
the existence of a Curie temperature, of a spontaneous magnetization 
below the Curie temperature, and the form of the susceptibility curve 
above the Curie temperature. However, the origin of the ‘internal 
field’ was not discussed in detail, and thus a statistical formulation in 
terms of atomic interactions was not possible. In 1925 Ising attempted 
such a formulation, although the form of interaction which he took was. 
somewhat empirical. In 1928 Heisenberg gave a quantum mechanical 
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treatment of atomic forces, and hence explained clearly the origin of the 
‘internal field’. He showed that in the first approximation the atomic 
interaction is proportional to the scalar product of the vector spin operators 
of the atoms. 

Both on the Ising and on the Heisenberg models it is correct to state that 
any atom of the lattice is subject to an internal field arising from the 
interactions with neighbouring atoms; the field however is not steady but 
fluctuating as these neighbouring atoms change their orientations. The 
Weiss theory corresponds to replacing this fluctuating field by its mean 
value; we refer to an approximation of this kind as ‘zero order’ since its 
properties do not depend in any way on the detailed structure of the lattice. 
We might expect this approximation to get better as the number of neigh- 
bours interacting with a given atom increased, hence diminishing the 
significance of fluctuations. 

Early attempts to take account of molecular interactions in solutions 
assumed as a first approximation that the molecules mixed randomly. 
The conclusions resulting from this assumption are mathematically 
identical with those arising from the mean field approximation in ferro- 
magnetism. Subsequently the Bragg—Williams theory (1934) produced 
equivalent results for the order—disorder problem. 

The use of a mean field and the assumption of random mixing are 
independent approaches to the same statistical problem, each casting 
its own light on the nature of the approximation. For more complex 
problems sometimes one and sometimes the other can be generalized, 
and we have therefore considered it worth while to give an account of both 
approaches. 


4.2.1. The mean field approximation 


The energy of interaction of two neighbouring spins on the Ising model is 
—Jo,,0,;. Hence the total action of neighbours on any spin 7 of the lattice 
can be replaced by an equivalent magnetic field of magnitude. 


eT ee en ements? 4%. '( L) 


the sum being taken over the q neighbours of the spin. We replace this 
fluctuating field by its mean value 


Ho= A %4)» ge oS) a, een 4 


where «,, a are the fractions of spins pointing in directions | and 2 
respectively. But a, «, can be determined from Boltzmann’s Principle : 


Serexpi iw (iH EA kilo os os =). i (8) 
OT 

Combining (2) and (3) we have an equation for Hy: 
mH {qJ =tanh [w(H4Hy)/kT], . . . . - (4) 
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which is equivalent to the equation derived by Weiss (1907). 
In the absence of an external field equation (4) may be written 


tanh Sues . (€=mH, |kT). 5 eae ek pa PRE EP) 
qd 


From fig. 1 it will be seen that when k7'/qJ >1 the only solution is €=0; 
but when k7'/qJ <1 there are two additional equal and opposite non-zero 
solutions for €. The magnetization is given by 
T= (ag 04) = “5 Ho 5. assish op Matec mnr ey) 
and hence a non-zero value of € in zero external field corresponds to the 
existence of a spontaneous magnetization. We must first prove that this 
non-zero value of € gives rise to the lowest free energy; this can be shown 
more readily by the approach in the next section. 
The Curie temperature corresponds to the appearance of a non-zero 
solution for € in (5), and is given by 


On Reet em ay 


Condition for the existence of a spontaneous magnetization in the mean field 
approximation. When C>1 the only solution of tanh €=Cé is €=0; 
when C'<1 there are two additional solutions (C=kT'/qJ). 


The form of the spontaneous magnetization can be calculated from (5) and 
(6). For low temperatures (& large) 


Iy~m[1—2exp(—2qJ/kT)]; ot fre OH 
but near the Curie temperature (€ small) 


Iy=nl 3(1— 7) | Pe ait aos i) 


— 
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The susceptibility is given by (@//0H) and can readily be calculated 
from (4) and (6). We have, 


(0H ,/0H)/qJ =sech? [(H + H,)/kT][1+ (0H, /dH) |/kT 
= [1+ (0H,/0H)\/kT[1—(w#H,/qJ)?]. . . . (10) 
The susceptibility in zero field for 7’> 7’, is obtained by putting H,=0 in 
(10); we derive the Curie-Weiss law 
0H, 
AST, a 5 (oe dH 
for which the inverse susceptibility is exactly linear. 

It should be noted that thermodynamic quantities in this approximation 
depend only on 7'/7'c, and hence are the same for all lattices. The approxi- 
mation is completely wrong for a linear chain (q¢=2) since we know that 
there is no singularity in this case (§3.1). It is somewhat better for the 
two-dimensional lattices improving steadily with increasing coordination 
number; from §3 eqn. (183) the values of k7'./q¢J for the honeycomb, 
quadratic, and triangular lattices are 0-506, 0-567, 0-607 respectively. 

There is no difficulty in generalizing the approximation to take account 
of distant neighbour interactions. The only change is in the Curie 
temperature (7) which is now given by 

Remy Goda teeGudpy is ioes- 0s 1 2 (12) 


4.2.2. The approximation of random mixing 


‘) = n2|k(T — To), acer) 
=0 


Let us consider all configurations with N, spins pointing in direction 1 
and NV, indirection 2. Thereare N !/N ,!N, ! such possible configurations, 
and when the positions of these spins in the lattice are completely random, 
the mean number of 1—2 nearest neighbour pairs is 


(Ny,)=258 52 (gN/2)=GNNAIN. . . . . (13) 


Our approximation is to assume that formula (13) holds even when there 
are interactions between the spins, so that we ignore correlations between 
neighbouring atoms. We can now write down the free energy for the 
Ising model (§ 2 eqns. (5), (9)), replacing N,, by its mean value from (13), 


Fyi= —NmH- ar kT In Ay(u,2), 
(14) 
N! 
ONG, (2) = se es tg iN NalN 
n (we ) ven N,!Ny! 
For large N we can as usual replace the sum in (14) by its maximum 
term. The logarithm of the general term is asymptotically equal to 


qN N. 
NinN-N,nN,-N,mN, + Nine + 4 7H ee eee 1) 


differentiating with respect to NV, (NV.=N —,), we obtain for the maximum 
term the equation 
(NV, wee N,) 
IN 


—InN,+inN,4+np+q it =) oteen ee, ise EO) 
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Writing N,=N(1—5)/2, N.=N(1+8)/2 we obtain 
148 _2mdl , 2a 
1-5 kT ais 
which is exactly equivalent to (4). Thus when H =0, for T’>7'¢ (15) has 
only one stationary value which corresponds to 6 = 0, or random orientation 
and mixing of spins. But when 7’< 7’, (15) has three stationary values. 
Differentiating the left-hand side of (16) we obtain for the second derivative 
of a typical term 
1 ss) Agee 1 1 =| 
Se te 148-2128 0s eee 
Hence when 7'> 7’, the stationary value at 6=0 is a maximum; but when 
T<T. the stationary value at 6=0 is a minimum, and the non-zero 
values given by (17) are maxima. 

We can now replace A,(y,z) by its maximum term given by (17) and 
(15), and hence derive the values of any thermodynamic functions in this 
approximation. Thus when H=0 the internal energy is given 
a N os i 7_ — qNI® 
Hence we derive the inne ance for the internal energy : 
Ni = 


(17) 


Ey= J + 2q —-— (19) 


[l1—4exp(-—27',/T)], | 


T'.—T small fig icc isha Pert _9 9/Te_ : (20) 
and positive, * ~ 2 fi Ben . 
DST ne lixg 20; 

For the specific heat we have similarly 


Tel ey li 


RP lars y= Nk ) exp (—27'./T), 


h 
T.—T small a ate T'¢ mere ee 
and positive, One: Sei E Coe (a > 1) |. 
ARES aH). 


It was noted by Bell (1952) that inconsistent results arise if the partition 
function is constructed using the mean field of the previous section and 
treating it in the same manner as an external field; the simplest method 
of obtaining a consistent partition function is the use of a combinatorial 
factor as in the present section. 


4.3. First-order Approximations 
Both of the approaches of the previous sections have their natural 
generalizations to an improved approximation. The mean field approxi- 
mation leads naturally to the Bethe approximation which treats more accur- 
ately the interaction of a given spin with its nearest neighbours. Im- 
provement of the approximation of random mixing leads to the quasi- 
chemical approximation. For the Ising model the results of the improved 
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approximations are mathematically identical; however, we shall again 
give a brief account of both approaches, observing that the extension of 
the approximations to higher order, and to the Heisenberg model, no 
longer produces identical results. 


4.3.1. The Bethe approximation 


The Bethe approximation takes exact account of the interaction of a 
given spin with its nearest neighbours, but uses a ‘mean field’ to take 
account of the interactions of these neighbours with other spins of the 
lattice (fig. 2). The energy of interaction of an elementary cluster con- 
sisting of the central spin and its q nearest neighbours is then 


q q 
Hy= — > IJ0o,0;,-mHoy,—mH, > 65, teen eae (22) 
jel jut 
Fig. 2 
H, 


i 


A cluster consisting of a central spin and its nearest neighbours for the Bethe 
approximation. The magnetic field acting on the central spin is the 
external field H; an additional field acts on the outer spins arising 
from the coupling with the remaining spins of the lattice. 


where as before the magnetic field is taken in direction 2. The partition 
function for this elementary cluster is 


Z, s(n, gd 4 sO Fy to 242+ + 8Oou nauess gree...) 
+ pi ( aq? z—a2 4. 8Q), p22 g-a2+1 4 sO, ON te Baer ei) se, ) 
= po? ( Wee ge 4 ae g2)a 4 12 ( poy? ge eke Blade SS a (23) 
where 


w=exp (—2mH/kT), py=exp (—2mH,/kT), z=exp (— 2J/kT). 
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The mean field H, is determined by the consistency condition that the 
probability of a spin pointing in directions 1 and 2 is the same for the central 
spin or any of the outer spins. Thus 


On ZG = OY in 7 ee os 
Yon te gaia 


Applying (24) to (23) we find after a little simplification 


a = (ee) 1 eee os 
be I+ py2 
Fig. 3 
ie) 
0.8 
0.6 
tanh 9/(q-1) 
tanh 
0.4 
0-2 
1.0 2:0 n 3-0 4.0 5-0 


Condition for the existence of a spontaneous magnetization in the Bethe 
approximation. When tanh (J/k7')<1/(q—1) there is no non zero 
solution for 7 ; when tanh (J/k7)>1/(q¢—1) there is a non-zero solution 


(n = mH, /kT’). 


In the absence of a magnetic field (= 1), H,=0 is always a solution. To 
find if a non-zero solution exists for H, we readily transform (25) into 

tanh [aeH,/(q—1)kT') 

tanh (sH,/kT') ) 

We have sketched the form of tanh (7/q—1)/tanh in fig. 3; it starts with 


the value 1/(q—1) at 7=0 and increases steadily to 1 at y= 00. Hence at 
any temperature 7' satisfying 


tanh (J/kP)<1j/(q—1)) sos oe ws ee LT 


=tanh(J/kT). . . . . (26) 
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there is no non-zero solution for H,; for all other temperatures the non-zero 
value of H, is given by the intersection of the curve in fig. 3 with the line 
y=tanh (J/k7’). The Curie temperature is thus given by 


2 

tanh (J/kT'c)=1/(q¢—1), i elgd = — -. (28) 
For a linear chain q=2 and the Curie point is at zero temperature, as it 
should be. As q increases k7'¢/qJ increases steadily reaching the limiting 
value of 1 as q+ oo. In fact if in (26) we allow g to become large and J 
to become small, q.J remaining finite, we recover the mean field eqn. (5). 
From the partition function for an elementary cluster, (23), we find that 

the mean fraction «, of spins pointing in direction 1 is given by 


Oy => = < (29) 
Z . 1+ Laas ae 2142 
Using (6) and (25) we find for the magnetization 
Spee 
fice vou eaiasaal ts RD ee, erg et ( SO) 


As H +0 when H, 40 we obtain a spontaneous magnetization. Near 7'=0 
we can derive from (25) an expansion for p, in powers of z starting with 
f,=221; from this we can derive a corresponding expansion for Io. 
Near 7'=T',, it is better to use (26) and expand for small H,; we find that 


if 
2 T 1/2 
Tye me} xo Lin alta—2)1(FE -1)} 


Hence the spontaneous magnetization curve is steepest for small g, and 
becomes less steep as g increases ultimately tending to the limiting form (9) 
given by the mean field approximation. 

The susceptibility can be obtained by differentiating J in (30) with respect 
to H, using (25) to obtain dy,/dH. For the susceptibility in zero field above 
the Curie point we put H =H,=0 after differentiation, and a substantial 
simplification occurs. We find that 


(mH [kT ~ 3(q —1) (In q/(q—2)] (7 =) 
(31). 


oa 2 G 
a formula first derived by Firgau (1941). 

The internal energy per spin can be obtained by considering the ‘bond’ 
energy of an elementary cluster, i.e. the contribution from the interaction 
of nearest neighbour pairs. If we then sum over all spins of the lattice 
each bond will be counted twice. Hence the total ‘bond’ energy of the 


lattice is given by 


ewe _ NgJ a as 
od 5, i a) = oo aero Son SRE ES ( ) 
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after simplification by means of (25). In the absence of an external mag- 
netic field this is the only contribution to the energy, and we shall confine 
our attention to this case. For 7'>7'c, H,=0 and we find 


NqJ (1-2) 
aie Se Bee ’ 34 
By 2 (1+2) See 
The value of the energy at the Curie point is thus 
E.=—NqJ/2(q—1),, (q@>=2). - e -eed 
and for temperatures just above 7’¢ 
2(q—2 i} oh 2 
By=Be+ NILE tinghtg— 201 (F -1)+0(7 -1)", 
4(q—1)° Tc Te (36) 
= Nk L2—*) fin g/(q—2) 2 +0 7~-}) 
y= Wk ETE tm gi(g—2)P +0 (ge -1). 


Below the Curie point we must calculate jx, from (25) and substitute in (33). 
For temperatures just below 7'c, expanding as in (31), we find that 


¢ ¢ 2 
By=Be— ws TA (In q ) (= -3) +0 -1) 


4(q—1)? q—-2/\T (37) 
_ yp F(4= 2)(89— 2) q ¥ Te 4 
c= wk") (in 4) niles 


The formulae which we have derived here using the Bethe approximation 
are most usefully applicable to loose packed lattices in which there is no 
interaction between spins of the outer shell in fig. 2. For lattices such as 
the triangular and face-centred cubic these interactions should be taken 
into account in the partition function Z,, and if they are ignored the 
approximation will be less accurate. The shapes of the thermodynamic 
functions in this approximation depend only on the coordination number q, 
and we shall discuss their detailed behaviour subsequently in § 4.7. 


4.3.2. The quasi-chemical approximation 


In 1935 Guggenheim, using the analogy of the equilibrium theory of 
chemical reactions, suggested an improvement of the approximation of 
random mixing in the theory of solutions. From (13) and analogous equa- 
tions for (N,,) and ¢N,.,) we find that in this latter approximation 


eee a) 
(Ni) (N22) 
However, if we consider the chemical reaction between molecules 
AA+BB=2 AB, 
ignoring vibrational and rotational effects, the equilbrium state is given 
by (Rushbrooke 1949, Chapter 11) 


(Ny)? 4lexp (= €o/k7’) 


(Ni) (Noe) - [exp (— €,,/k7’) [exp (— €g9/kT’) | 


(38) 


(39) 
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Guggenheim therefore suggested that (39) should be applied to the theory of 
solutions as an improvement of (38), the right-hand side of (39) being in fact 
equal to 427; he also proposed the term quasi-chemical to characterize the 
approximation. 


From (39), using the fundamental relations in eqn. (2) of §2, we easily 
deduce that 


Nye) _ — 2+ v/[24(1 + doar) — 42% or,09] 


qN 2(1 —2?) eo) 


Now the internal energy of a solution (fixed N,, N,) is given, apart from a 
zero energy term, by (N,,) « (§ 2, eqn. (3)); also 


F) a 
Ey=—kI?—,(InZy)=e2—(InZy). . . . . (41) 


Hence using (40) and (41) we can integrate to find InZ,; the constant of 
integration is fixed by considering the value at z=1 (7’= oo) where random 
mixing occurs. We thus find that (except for zero energy terms) 


B+a,— B-ata, |} 
InZy= —3qN SS se 


aie (a, In a, + a, In a2) 
B= v/[Aeryorge + (044 — X%)”].- 


(Details of the calculations are given in Rushbrooke (1949), Chapter 18, 
§ 6.) 

To determine the thermodynamic behaviour at different temperatures 
we must plot the free energy as a function of the concentration «, of one of 
the constituents. The results obtained are equivalent to those of the 
Bethe approximation. Above a temperature 7’, given by (28) the free 
energy minimum corresponds to «,=4; but below 7’, there are two values 
of «,, (+8) and ($—8) (840), giving rise to minima in the free energy. 
Below 7’, the stationary value at «,=4 becomes a maximum and does not 
correspond to astable state. The results can then readily be re-interpreted 
for the Ising model. 

The partition function for the Ising model for H = 0 at temperatures above 
the Curie point is readily obtained from (42) by putting o,=a,=$; 


sti9 
inZy=4qN | ——*in2+In(1 +2) | oes. (43) 


From this we can write down the critical value of the entropy 


a 7 Te ee gig 2) fn 2 
Sem 2 (Pin Zy)e=Wk| n2+ 4in(1 *) fain (1 )] 


(44) 
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Subsequently, Fowler and Guggenheim (1949) (see Guggenheim (1952)) 
and Li (1949) developed independently an alternative approach to the 
quasi-chemical approximation by deriving an expression for the combi- 
material factor g (N ; N,, Nj»), 

e (tq)! Ning |e 

Ns No Nwl= Fw wT | * + 
This provides a simpler method of calculating thermodynamic functions on 
the quasi-chemical approximation. The formula (45) is only a valid 
approximation in the neighbourhood of random mixing, i.e. at high tem- 
peratures (Prigogine, Sarolea, and van Hove 1952) ; but at low temperatures 
it gives rise to minima in the free energy corresponding to phase separation, 
and these provide the appropriate approximation in this region. 


4.4. Higher Order Approximations 


There have been many attempts at improving the Bethe and quasi- 
chemical approximations. The most obvious method of extending the 
Bethe approximation is to take rigorous account of an additional shell of 
atoms, and to take account of the interaction of this outer shell with the 
remaining atoms of the lattice by the use of a mean field. In principle 
the approximation could then be further extended by steadily increasing 
the number of atoms which are treated rigorously ; in practice the calcu- 
lations rapidly become cumbersome, and the approximation has not been 
pursued beyond the second stage. 

Using the combinatorial approach described at the end of the last section 
the quasi-chemical approximation was similarly extended by a number of 
writers (Yang 1945, Li 1949, Hill 1950, Guggenheim and McGlashan 
1951). However, comparison with the results of known exact solutions 
in two dimensions (§3) indicates that the improvement proceeds rather 
slowly. 

We have therefore preferred to present the approach of Kikuchi (1951) 
since this leads most simply to the maximum improvement; many other 
approximations will then be seen to arise as particular examples 
of this approach. We shall briefly survey other ‘closed form’ approxi- 
mations in § 4.4.3, and analyse in more detail the nature of the various 
approximations in § 4.6, 


4.4.1. The method of Kikuchi. Combinatorial factor 

Kikuchi (1951) considered the process of constructing an ensemble of L 
assemblies, one unit being added to each assembly at the same time. The 
ratio of positive and negative spins is kept fixed in the ensemble, as are 
certain other ratios such as those of different types of bond, or of configura- 
tions of a square, and a direct attempt is made to calculate the number of 
possible complexions of the ensemble subject to these restrictions, and 
hence the entropy per added unit of the assembly. The method will be 
understood more precisely on reference to table 1: 
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Table 1 
Fraction of Configuration ae of B, 
Configuration} Occurrences ccurrences Gas 
aie! oa Y1 1 
5_§ Toa abe 
Cans Y3 l 


contguracon| Brvation of | 5, || contgumation | Brsetion of | 
(Cran ae @ 
| | Ail 1 | Ww, il 
O—-& o-® 
eee ® ‘S) 
| | ete a9 4 | | Wy 2 
a Ore O-e 
aan ® 
bere | rete 25 4 | We 1 
iy mae, 
-— ® © 
| |. ete ZA 2 | | Wa 2 
re Orne Mao ae, 
Ce. S) 
bo Ble te Be 4 Ws l 
Sa) e—-O 
ae) © 
| | 6 1 ee We 1 
Sane, Crane 


(¢) (d) 


where various elementary configurations are classified together with their 
fractions of occurrence. 

The simplest configuration of a single spin is shown in (a); clearly the 
fractions satisfy the relation 


L,+%,=1, . . ° . ° . . . (46) 
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so that one degree of freedom arises from this configuration. The next 
configuration is a bond shown in (b); the third column records f; the 
number of equivalent configurations of a given type. Now the following 
relations hold 


Yt Y2=*p \ OS ne anes 
Yor Y3= 2, 


so that this configuration introduces one additional degree of freedom. 
The next configuration which we shall consider is a square shown in (¢); 
the connecting relations are now 


2j+2%,+23, =Yy, 
Za + 2% +24 +25 = Yo, eee | ns Sh 
2% + 225 +2, =Ys3; 


so that this configuration introduces three additional degrees of freedom. 
If we now confine our attention to (a) and add a new unit to a lattice 
keeping the «, constant, the number of possible complexions is 


L! D} 


(a,L)!(@b)! Xz’ Sa 


2 
where X,=[] (x, L)!. 
i=1 
The entropy per unit is obtained by taking the logarithm of (49) for large L, 
and dividing by L, and is hence 


S(2,)=k( —2,Inz, —2,Inz,). 


The internal energy must be written as a function of the ;, and hence 
correspondingly the free energy; minimizing this with respect to the free 
variable, we obtain the equilibrium distribution of the 2,;, and hence all 
thermodynamic properties can be deduced. This corresponds exactly to 
the mean field approximation. The next stage introduces y; and configura- 
tions (b) and attempts to calculate S(x;, y;); minimization of the free energy 
with respect to the free variables gives the Bethe approximation. The 
third stage introduces z; and configurations (c) and attempts to calculate 
S(x;, Y;, %); minimization of the free energy now leads to an improved 
approximation. By taking larger elementary figures the method can clearly 
be extended. 

We shall first consider the calculation of S(«;, y;) for a linear chain as 
illustrated in fig. 4. At any stage in the construction of the ensemble 
v,L of the last units added are (+) spins, and x, are (—) spins. Among 
the x, units we can select y,/ at random and attach (+) spins to them, 
giving y,L(+ +) bonds; tothe remaining y,/ of the x,L units we attach (— ) 
bonds giving yl (+ —) bonds. The combinatorial factor associated with 
this procedure will be (7,L)!/(y,L)! (y,L)!... We now attach (+) and (—) 
spins in an analogous manner to the 2,1 (—) spins with combinatorial 
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factor (x,L)!/(y,L)!(y3L)!. Hence the total combinatorial factor arising 
from the addition of a new unit is 
(selene) Xz 


(CHEV CM AMLIOMD ip ke au on peammciisose 


3 
where Y,=[[ (y,L)!4. The corresponding entropy is 
i=1 


2; 3 
See.a=k| Sana > Binlnas eee OL) 
4=1 tJ 


The relations (47) will ensure that of the new units added x,L are (+) spins 
and Xb (—) spins. For the one-dimensional chain there is no approxi- 
mation in the procedure and formula (51) is exact. 


Fig. 4 
NEW BOND 
| ADDED 
(Grane Orr OF Ran ae eae Oa On Sh aU AT TICE 
0 O Sp aoa es O-——@® 2ND LATTICE 
lo O Oe O-——@® LTH LATTICE 


Building up an ensemble of linear lattices. The last bond is added in such a 
manner as to maintain a given ratio of + to — spins in the ensemble. 


We now attempt to apply the same procedure to the two-dimensional 
quadratic lattice as illustrated in fig. 5(a). We wish to find the combina- 
torial factor arising from the addition of a unit A, when units B and C have 
the distribution x, and the bonds AB and AC have the distribution y; in the 
ensemble. Let us first add bonds BA’, CA” (fig. 5 (0)) as in the one-dimen- 
sional case, the combinatorial factor for this procedure being then (X ,/ Y ,)?. 
Both A’ and A” will have the distribution x, in the ensemble; however, to 
add the unit A correctly we must examine the detailed distribution of (+ ) 
and (—) spins in the ensemble for every pair of configurations A’ and A” and 
only count those for which these detailed distributions are identical. Hence 
the expression (X ,/ Y ,)? must be multiplied by a correcting factor to obtain 


260 C. Domb on 


the appropriate combinatorial value. This factor cannot be assigned 
exactly in two dimensions; but Kikuchi suggested that a reasonable 
approximation could be obtained by assuming that A’ and A” were random 
with respect to each other, in which case the factor would X,/L! (a single 
distribution of the unit A among members of the ensemble is required 
instead of two independent distributions A’, A”). We thus finally derive 
as an approximation to the combinatorial factor for the quadratic lattice 


the formula SP eh ce aN Re ee ees 


Fig. 5 


(b) (c) 

(a) An intermediate Stage of constructing a two-dimensional quadratic lattice. 
(b) Bonds BA , CA” are added independently each having the distribution 
y;- (¢) A and A” are added independently so that square A’ BDC 
has distribution z,, and angle A” CE has distribution w;. 


More generally for a lattice with coordination member q it can be shown that 
the corresponding combinatorial factor is 


Kt Y OF DG aY es ae ee eS 
and the entropy 


' 2 3 
S(«;, ¥:)=(q-1) > wIna,—¢/2 > BG, eee (54y 
i=1 t=1 
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To improve the approximation we build up the lattice using a square as 
the basic elementary figure. We introduce the angle configurations in 
table 1 (d) to help in the construction of the lattice, but they do not enter 
into the final result. The unit A is to be added so that the square ABDC 
has distribution.z,,.and the angles BDC, ACE have distribution w; in the 
ensemble. As before we construct the configurations BDCA’, ECA” with 
appropriate distributions (fig. 5(c)) having combinatorial factors 

6 6 
W ,/Z,, Y,|/W, respectively (4, TL @£)!", W,= TLL) Dy 
i=1 i=1 

To obtain the correct combinatorial factor we must now pair together 
complexions for which the detailed distributions in the bonds CA’, CA” are 
identical. Again we cannot correct exactly, but we assume as an approxi- 
mation that the bonds CA’, CA” are uncorrelated ; we must therefore divide 
the total number of complexions obtained by independent pairing by the 
number of possible complexions of a bond CA’, C having given distribution 
x;, and we thus obtain for the combinatorial factor 

Wie ae eee 

— eS aH YY B(X,Z_). ee a pe AGED 
The entropy in this approximation is then given by 


2 (3 6 
S (x, Ys “)=k| - >» x, In x;+ 2 > Byy,ny,— » vein |. (56) 
(=1 jai = 


Kikuchi (1951) applied a similar method to the simple cubic lattice, 
obtaining for the combinatorial factor Y,°L!/(Z,>X,’), and for the 
entropy ; 

2 3 
S (aj, Yj, %) =k | —7 Y wInaj+9 > Bylny,—3 2 yeiln | - (57) 
i=1 i=1 i= 
Similarly he obtained improved approximations for the triangular and 
face-centred cubic lattices. 

However, the above argument is not easy to follow in the more complex 
cases, and in‘a subsequent paper, Kurata ef al. (1953) suggested an alterna- 
tive development of the approximation which enables the appropriate 
combinatorial factor to be written down more readily. This is briefly 
illustrated in fig. 6; the elementary figures from which the lattice is 
constructed are all drawn in separately, and the factor necessary to select 
only configurations in which all spins at a given lattice point are identical 
is estimated. It is found that, for example, if the lattice consists of £N 
squares, sNV bonds and rN points the combinatorial factor is 

(Liste Y M-8/X  a-28+0 Zt, eres ae (58) 
and this will be found to include (55) and (57) ((=1, s=2, r=1 and t=3, 
S=2o,7 = 1), 

A ca inal analysis of the Kikuchi method together with a technique 
for determining the combinatorial factor in the most general case is provided 
in a comprehensive series of papers by Hijmans and de Boer (1955, 1956). 


P.M.S. T 
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4.4.2. Evaluation of thermodynamic properties 
The internal energy per spin of any configuration is given by 
E! = —mH (ay — 2) — 394) (Yy + Ys — 2Y2) =H (wy — Xq) — 39d + 2S Yo 
e e e . . . . (59) 


and the equilibrium values of a, y;, 2; are determined by minimizing the 
free energy Ff! =(H#!—7'S) with respect to the free parameters. We have 
noted that relations (46) (47) and (48) leave five degrees of freedom. It is 
useful to choose free parameters which have a direct physical significance ; 


thus 
Sh Pie \ pe Sees 
£5 = %— 2% 
Fig. 6 


From the above assembly only configurations in which all spins at a lattice 
point are identical must be selected. 


have the significance of long-range order parameters and will be zero for a 
disordered state in the absence of a magnetic field. For the other free 
parameters we choose Ys, 2, and z,; the remaining dependent variables can 
then readily be expressed as linear combinations of these e.g. 


2z,=14 €,—2&,—4y,4 22. ee er (61) 


Taking as an example the second approximation for the simple cubic 
lattice, (57), we obtain: 


OF" /0E,=0, (a4/2y)" (Ys/Y1)° (2/26) =p, (i) 

OF" [y= 0,[Yo"/(YrYs) ¥ [2126] (2225)? =2%, (ii) 

oF! /dz,=0, re = 225, (iii) el oa) 
OF! /dz,=0, gPo* =2,72,2,, (iv) 

OF'/0g,=0, (29/%5)" =21/Z, (v). 


These equations can in general be reduced to a single algebraic equation in 
one variable, and the thermodynamic functions can be derived from (57) 
and (59) by numerical computation, 

In the absence of a magnetic field (= 1) we look for critical phenomena 
corresponding to the appearance of long-range order. The disordered state 
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at sufficiently high temperatures will be characterized by G6, = 0, and 
eqns. (i) and (v) above become identities. The remaining equations 
simplify considerably, and an explicit solution is possible in terms of a 
single parameter ¢ as follows: 


2 = [(3¢ + 1)/($+3) P/¢, 

44 = [2(6? + 66+ 1)] 7, 

23 = P24; 

Yo= (36 + 1)zq. 
(For details of the calculations see Kikuchi (1951).) 

The ordered phase is more complicated since it depends on five parameters, 

and the detailed properties must be determined by numerical computation. 
However the critical point (and hence critical values) can be determined by 
examining eqns. (62) for the onset of non-zero solutions for ¢, and €. By 
analogy with the conclusions of the Bethe approximation (31) we may 
assume that just below the Curie point the values of d5€, and dé, will be an 
order of magnitude larger than the changes in any of the other parameters. 
Hence if we consider a small change in equations (i) and (v) taking account 
only of 6€, and 6€,, and eliminate 6é,, 6€, we shall obtain an additional 
relation to (63) from which the Curie point can be determined. We find 
that (putting £,=&,=0 after differentiating) 


(63) 


Bie 22. oy 
Of. 2 a = cal el ©2) | &4 
From (63) and (64) we deduce that 
odc=3, kT ./J =4-6097. tee eee. tote O0) 
The specific heats on the high and low temperature sides of 7’; are given by 
c, /k=0-389, \ (66) 
c_/k= 2-902. 


We shall defer until §4.6, §4.7 a detailed discussion of the physical 
conclusions arising from the application of the Kikuchi method to various 
lattices. 

It is perhaps worth observing that by giving some parameters fixed 
values corresponding to a high temperature and varying the others, some 
estimate can be obtained from this approximation of non-equilibrium effects 
such as ‘frozen in’ order. 


4.4.3. Other approximations 

The first approximation which registered a marked improvement over 
the Bethe approximation was the variation method introduced by Kramers 
and Wannier in 1941. We have seen in § 3.2.2. that the partition function 
for the Ising model is determined by the largest eigenvalue of a symmetric 
matrix. By assuming an eigenvector § dependent on a number of para- 
meters and minimizing §’V&/(§’§) with respect to these parameters a 


T2 
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useful estimate can be obtained of the largest eigenvalue. Kramers and 
Wannier applied this method to the quadratic lattice and obtained results 
which were exactly reproduced by Kikuchi (1951) from eqn. (56). ‘The 
method was extended to the simple cubic lattice by ter Haar and Martin 
(1950), but the detailed calculations for three-dimensional models are 
considerably more complex. The accuracy of the results in this case seems 
comparable with that of Kikuchi. (The technique of ee a various 
approximations will be discussed in § 4.6.) 

In their 1953 paper Kurata et al. considered the possibility of using the 
Kikuchi method with different types of basic elementary figure. They 
pointed out that the choice of a larger basic figure could result in little or no 
improvement of a preceding approximation. In fact they suggested that 
in order to improve an approximation it is necessary to adopt a basic figure 


Fig. 7 


rae 
er 


(d) 


Examples of open and closed clusters. (a) Represents an elementary figure, 
(b) are open and, (c) are closed with respect to (a). 


which is ‘closed’ with respect to the basic figure in the preceding approxi- 
mation. The terms ‘open’ and ‘closed’ are defined as follows. If the 
fraction of occurrences of a cluster B (table 1) can be (cannot be) derived 
rigorously from those of its sub-cluster A, B is said to be open (closed) with 
respect to A. For example, in fig. 7 (b) are open and (c) closed with respect 
to the corresponding (a). Subsequently Hijmans and de Boer (1955) showed 
that this statement is not completely valid, and that a more careful analysis 
is required; it can nevertheless serve as a useful guide. The approxima- 
tions of Yang (1945), Li (1949), Hill (1950), Guggenheim and McGlashan 
(1951) and Fosdick and James (1953) correspond to the choice of a particular 
elementary figure, and are intermediate between the first and second order 
approximations of Kikuchi. 

Prigogine, et al, (1952) suggested an ingenious improvement of the quasi- 
chemical combinatorial factor (45) based on a comparison with the known 
exact factor for the Onsager solution. However, they did not undertake 
any detailed thermodynamic calculations other than the evaluation of the 
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Curie point (which resulted in a substantial improvement) since such cal- 
culations would seem to be quite involved. Itis difficult to assess the precise 
status of this approximation. 

It is similarly difficult to draw definite conclusions regarding the approxi- 
mation of Cowley (1950). Following Zernike (1940), Cowley attempted 
directly to calculate the order between different pairs of lattice sites, but 
he assumed that the order parameters corresponding to pairs of sites at 
different distances are independent. The resulting specific heat singularity 
is infinite of order (7’.—7')-1?, but the estimate of the Curie temperature 
is no better than the mean-field approximation. The method of Yvon 
to be described in § 4.5.2. seems to provide an exact description of this ap- 
proach with no arbitrary assumptions. 

The ‘constant coupling’ approximation was introduced by Kasteleijn 
and van Kranendonk in 1956. Its most important achievement was the 
provision of a ‘quasi-chemical’ approximation for the Heisenberg model. 
However, the authors showed that the same technique could be applied to 
the Ising model. 

Recently Levitas and Lax (1958) introduced a technique which is a com- 
bination of the ‘spherical model’ and the Bethe approximation (the 
‘spherical model’ replaces the Ising interaction by an alternative inter- 
action for which the partition function can be calculated exactly; see 
Berlin and Kac 1952). In two dimensions the results obtained are 
appreciably wrong; but for the simple cubic lattice in three dimensions a 
remarkably accurate estimate is obtained for the critical temperature. 
However, the precise status of this approximation is again difficult to 
establish. 

4.5. Methods Based on Series Expansions 

We have seen in §3.6 how exact expansions for the partition function 
can be formulated at high and low temperatures. The derivation of 
substantial numbers of terms of such expansions requires a specialized 
technique which is described in detail in §5. We shall quote the results 
of such derivations, and discuss how they can best be used to provide 
information regarding the physical properties of the model. 

Two alternative procedures have beenemployed. A particular selection 
of terms which can readily be evaluated has been summed to provide 
closed form approximations similar to those discussed previously ; however, 
the nature of the approximation and of the terms ignored is known precisely 
at each stage. This is the basis of the methods described in §§ 4.5.1, 4.5.2. 
Alternatively, direct use can be made of the known terms of a series expan- 
sion to assess its asymptotic behaviour, and hence estimate its sum. This 
approach is described in §§ 4.5.3, 4.5.4. 

In the first method extensive use is made of the Mayer theory of cluster 
integral expansions (Mayer and Mayer 1940, Chs. 13, 14). The theory 
was first applied to the Ising problem by Fuchs (1942), but a substantial 
amplification was made by Rushbrooke and Scoins (1955). A somewhat 
analogous approach to the problem was initiated independently by Yvon 
in 1945. 
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We shall quote here for convenience the main forumlae needed from the 
Mayer theory. This theory applies to a gas in which the interaction energy, 
V(7,;) between any pair of molecules 7, j is a function of their distance apart, 
Nis and can be taken as zero beyond a certain range. Introducing the 
functions 


Sa=T(rg) =exp [—V (ty) [RT] 1, ee 
the grand partition function can be spa in the form 
2 m= FAG) 3 ee a 
V t>1 


where b, are cluster integrals defined by 


Voy= 73 -.. | fi foes Seka 


Each cluster integral corr stants to a different connected graph among the 
t labelled points, and the sum must be taken over all possible graphs. A is 
the absolute activity, and is determined in terms of the density p by 
0 eine ; . 
p= a7 ne) = ¥ aby hie be ge ee 
It is one of the most remarkable properties of the Mayer theory that the 
elimination of A between (68) and (70) leads to the formula 


pme=e(1— 3 Ao Re"), St ee 


ul 
where the f,, are ‘irreducible’ cluster integrals. These are defined similarly 
to (69) by 


1 
Bam ay | oe | false Satta dren . . (72) 


e 


but the sum in (72) is taken over graphs among (w+ 1) labelled points which 
are multiply connected (when w> 1), i.e. in which any two points of the 
graph are connected by at least two paths having no common point other 
than their ends. The elimination of all but multiply connected graphs 
is a substantial simplification. 

We shall first describe the approach of Fuchs and Rushbrooke and Scoins 
because formulae are derived for the regular assembly problem which 
directly parallel eqns. (68) to (72). 


4.5.1. The method of Rushbrooke and Scoins 


In order to simplify the analogy between the Mayer theory and regular 
assemblies it is better to consider the problem in which 1 units are molecules 
and 2 units are vacancies, so that in the notation of §2.1 €9=€9.)=0, 


€=-—e,,. From the general discussion in § 2 the function 
Ay(s2)= D> g(N3Ny,Ny)pMe% 2... (78) 
My, Nis 


can be regarded as a configurational partition function from which all 
thermodynamic properties of the Ising model and other corresponding 
physical assemblies can be derived. 
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To introduce the Mayer formalism for systems restricted to certain lattice 
sites, we observe that interaction energies arise only for certain discrete 
values of the distance r,;,. In fact, V(7;;) now has the three possible values, 
+ 0 (when 7;;=0), €,, (when i and j occupy neighbouring sites), and 0 
(otherwise). Hence instead of (67) we shall now expect 


—1 if the two systems are on the same lattice 
point, 
fig= 4 27-1 ifthe two systems are on adjacent lattice . (74) 
points, 
t 0 otherwise, 


where z= exp [e,,/2k7'] as in § 2 (4). 
Making the substitution 


Se Le) emer ys eee ees ee LO) 
in (73), and using the relation [§ 2 (2)]2N,,+N,.=qN,, we obtain 
Wie) a> GIN, Nio)(uet)*2(1 +f)... (78) 


Is 


To complete the analogy with the Mayer theory we must now label the V, 
systems (thus introducing a factor 1/N,!) and rewrite (76) in the form 


a)N, 
Ay (im2)= > sq +i) ey ah 
N, ae t>j 
where the /;; are given by (74), S denotes summation over discrete values 
instead of integration, and the product is taken over all pairs of molecules. 
The Mayer formalism can be applied exactly to (77) (Mayer and Mayer 
1940, sections 12b and 13 a), and we obtain formulae equivalent to (68) and 
(70) 

Tee) > OAc) ee (A= pet). 2. (78) 

t>1 


“ SSE i ae ne) 


tS1 
where the 65,(z) are now cluster sums defined by 


Nb(2)= 5 DShabyn----Sw cement, go (80) 


the connected graphs being defined as in (69). Comparing (78) with §3 
(146) or §5(115) we see that 
2!2b(z) =g, (2). SA Pcneg O (81) 
As a simple application of (80) we shall evaluate the first three cluster 
sums corresponding to the graphs -, —,/\ /\. For each type of graph the 
sum (80) must be evaluated for all possible configurations of the systems on 
the lattice, remembering to take account of contributions which occur when 
two or more systems occupy the same lattice point. 
(i) - clearly gives 6,(z)=1. 
(ii) , ,» When 1 and 2 are nearest neighbour lattice points the graph 
gives a contribution f. There are Nq/2 such links in the lattice, and taking 
account of interchange of 1 and 2 we obtain a contribution to b,(z) of gf/2. 
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When 1 and 2 are the same lattice point the contribution is —}. Hence 
bo(z) = (af — 1)/2. 


1 . . 
(iii) ,/\3 When all the vertices are nearest neighbours the contribution is 
pf? where Np, is the number of triangles in the lattice (§ 5.1.1); when two 
of the vertices coincide we obtain a contribution —qf?/2; when all three 


1 
coincide the contribution is —}. ,/\, When 1-2 and 2-3 correspond to 
nearest neighbour links the contribution is 


g(q—1) os, 
cer ree be 


when 2 and 3 coincide the contribution is (¢/2)f?; when 1 and 2 coincide 
the contribution is —qgf; when all three coincide the contribution is $. 
Hence 

b3(z) = psf? + a —") pp—af +t. 

The problem of deriving higher terms depends on the enumeration of 
different types of configuration on lattices, and is discussed in § 5.4.1 (see 
also Katsura 1958 and Temperley 1959). 

Continuing the analogy with the Mayer theory Rushbrooke and Scoins 
used (78) and (79) to eliminate A and hence express In A as a power series 
in a, analogous to (71), 


u 
InA =a,(1- >} —— gay oa ae 
nA (a4, 2) a ( 2 x3 By oy ) (82) 
B,, are now irreducible cluster sums 
1 
Be = a! 2" Shs cee Siw . . . . . . (83) 


taken over multiply connected paths among the (w+1) points. £8, is 
evaluated by keeping system | fixed and allowing systems 2, 3....(w+1) to 
take up all possible positions on the lattice. Rushbrooke and Scoins then 
proved that, despite the multiple occupation of lattice points, non-zero 
contributions are obtained only from configurations in which the lattice 
points are multiply connected. £,(z) and £,(z) are thus derived from the 
graphs —, /\ respectively, and proceeding as for b,(z) we find that 
B,(z)=9f—1 \ 
Ba(z) = 3psf*— 3qf*/2— 1/2. Sic 
It will be noted that each £,,(z) can be divided into various contributory 
terms; pure numbers arise from configurations representing multiple 
occupation of single lattice points; terms linear in q arise from con- 
figurations representing occupation (single or multiple) of two nearest- 
neighbour lattice points ; terms linear in p, arise from configurations repre- 
senting occupation (single or multiple) of three nearest-neighbour lattice 
points forming a triangle; and soon. Thus we may write 


By=BY+ £ Bu + PoB,. + DaBy® + Psa, + pads : % (85) 
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where Nps, Np,, Nps, represent the numbers of various types of 
multiply -connected configuration on the lattice (§ 5.1.1). Hence 
substituting in (82) and summing each particular class, we may write 


In A (a, 2) = Hi, (a,2)+ 2H, (x; 2) + pglf (ay, 2) 


Se Dyld (Oy, 2) + Datta (Asse) ees sce 8 (80) 

To determine H, (a,, 2) we may consider a pseudo-lattice of N isolated 

systems (fig. 8a) for which g=0, p,;=0, py=0.... but for which all the 

conditions of the Mayer theory apply. For such a _pseudo-lattice 
A(u,z)=1+p; hence a,=p/(1+p), w=0,/(1—o,), and 


on, 2) — (t= a). oe a ke ee (8a) 
Using (82) and (85) we deduce that 
Ge Ocae a irene) sie ubyuaer non wrt a {BB} 
Fig. 8 
OMe. Oo—o o-0 O—Oo 


OF ©) ae On) Oa) 


Oe O at Oa OS OO 


eA 
eS 


(c) 


Examples of pseudo-lattices to which the Mayer theory is applicable. (a) 
Isolated spins. (b) Isolated pairs of spins. (c) Isolated triangles of spins. 


Similarly to determine H, («,,z) we consider a pseudo-lattice made up of 
N/2 isolated pairs of points (fig. 8b) for which y= 1, p3=9, p4=9. 
For this pseudo-lattice 
In A (pu, 2)=3 (1 + 2u2+ ?), 
ee p(w +2) toate as ew S9) 
ple 
and we deduce that 
Hi, (04,2) =In [(1 —o,) +427 (20,-14+2)]. - - 22+ - (90) 
(C= [1+ 40, (1 — 04) f ]) 


Again by expanding (90) as a power series in «, we can calculate 6. 


270 C. Domb on 


We can determine H, («,, 2) for any simple /-sided polygon by using a pseudo- 
lattice made up of N/l isolated polygons. For such a pseudo-lattice 
q=2, p,=1/l and all other p,, are zero; also for this assembly (§ 3.2.1) 


InA (p,2)= rin [A,!-+ As!], 

7 (91) 
Ava ey [ee ae js 
and hence H,(a,,2) and 8, can be determined. This powerful technique 
of rearranging the terms of (82) in the form (86) and using finite clusters to 
find the initial terms of the latter series seems first to have been used by 
Yvon (1945). 

Rushbrooke and Scoins pointed out that stopping the series (86) after the 
first term corresponds to the approximation of random mixing ; inclusion of 
the second term corresponds to the quasi-chemical approximation; and 
the addition of further terms gives rise to improved approximations. They 
calculated the critical points for various lattices using such improved ap- 
proximations, and found the results in general better than the approxi- 
mations of Yang—Li-Hill and Guggenheim and McGlashan (§ 4.4.3) but 
not as good as the approximation of Kikuchi (§§ 4.4.1, 4.4.2). In fact the 
method as outlined so far can only take account of homogeneous con- 
figurations (§ 5.2.9) in which all vertices have the same coordination number ; 
a great improvement would probably be registered by calculating terms in 
(86) corresponding to non-homogeneous configurations (e.g. H;, (a, 2)), 
but this would first require a generalization of the theory to apply to non- 
homogeneous lattices. 


4.5.2. The method of Yvont 


It is unfortunate that Yvon’s papers (1945, 1948) describing his method 
were published in a rather obscure journal. His approach is particularly 
comprehensive, and has not received the attention which it deserves. The 
extension and application of the method to a variety of physical problems is 
due to Fournet (1952, 1953, 1957). 

Yvon starts from his own development of the cluster integral theory 
enunciated in a classical memoir on density fluctuations (1937). Con- 
sidering an assembly of NV interacting molecules in an external field U(r), 
and using the notation of the previous section for f,,, the configurational 
partition function is given by 


N 


Baym [os TIL 0 +h) Ear, Aa Res 


ay J ist bi=3 
<9 
where f=exp [— U(r,)kT]. 


pt na a a eR ae ae 

} For the cluster integral development in this section it is convenient to use 
the numbers 1, 2 to refer to different systems of the assembly, and not as 
previously to different states of individual systems. 
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For the density at a given point 1 we have similarly 


N 
ely= |. (TIL +idELen (98) 


Jt 4 j= 
(i<j) 


The Product TL 1+f,;) willbe denoted by 4,, y. This can be expanded 
‘Jat 
<j 
in the usual Tanner in terms of clusters of various order; a grouping can 
then be undertaken with special reference to the par agiks system 1. The 
first group consists of all clusters not involving system 1, and these sum to 
A,3 yx. The second group consists of all clusters in Bren system | is 


part of a 2-cluster (e.g. fis fsa fas fos OF S13 feats fs); the terms involving f,, 
sum to f> A434; ___ y and there will be (NV — 1) terms of the same kind which 


we denote by }f,,434; __y- The third group consists of all clusters in 
2) 


. ( . 
which system 1 is part of a 3-cluster; as in the Mayer development it is 
convenient to label the systems other than 1 and divide by (t—1)! for a 
t-cluster; this group may then be written 


2 (aSiefis + Sietes + BhieSisSes) A 
the sum consisting of (V —1)(N—2)terms. Thus we may write 


SEG ile alta (Sfrates Hh iates +o Sitistes) 


Age ont: . . . . . (94) 
Substituting (94) in (93), ah using (92), we find that 


piZy =), [Zy_1 aa Zn-2( (1 a 1/N) | fosis dr, 


+ Zy_,(1—1/N) (1—2/N) | [i Tiehethats tihahate) 


+JoJs dr. dt3+ aersie t's (95) 


where ae | | sian: 


In the limiting case of large V, assuming Zy ~ Z”, and neglecting terms of 
order 1/N, (95) can be written 


Pah [ 1+ [fs q2dt.+ | [ (kfiefist+hiefes + b.fiefis Ses) 12 Is 472473 + 
(98) 
where q,=41(Z: 

Considering as an example the equilibrium distribution of a fluid under 
gravity Yvon concluded that the series (96) was very slowly convergent. 
However, treating the series as an integral equation for the determination of 
an unknown function q, (2, y, z), the solution could be obtained by inversion 


as follows :— 


ts=rs| 1 | fevedra | [ (—Malie+ Lfahaledeapedrades— «| 
; bpiemceeat bee 162497) 
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Finally, taking logarithms in (97) a series is obtained, as in the Mayer theory, 
in which only multiply connected clusters enter : 


Ing, =Inp,— | fispedre I] hiefisfesP2P3d72473+ .... - (98) 

=O ORE 
Here C=In[N/Z ({f,dz;)] is a function of temperature and volume alone, 
and independent of position. In the series (98) the general term is the sum 
of all multiply connected clusters of order ¢ divided by (¢— 1)! remembering 


that the system 1 is specialized and that the remaining (t—1) can be 
permuted. For example the fourth term can be written symbolically 


hea ] 3 3 3 
pf eS) oi bx 
2 4 2 4 Pa 4 2 4 


<5 intense (99) 

Considering again the example of a fluid under gravity Yvon concluded that 
the series (98) is rapidly convergent. 

The same procedure is now applied to the correlation p,, between the 
systems 1 and 2 using the equation 


Pel v=. JIL, (1+f;) Il dr. ° evpemen 00h 


(t<j) 
In all expansions based on (100) systems 1 and 2 are given special attention, 
and instead of (94) the appropriate expansion is 


Arzs....w=Asa....v + 2 (hia tfos + Sisfos) Aas eg N 
+ & holst Siasfort +++.) Ase Peed Nos a OLD 


Substituting in (100) we obtain, in the limit of large NV, 


Py2 = 192 exp [— V(ry2)/kT] {1 te J (fis +Ses +hisfes) Ys ITs 


ca [i TP f Pho ceeee oeaean deme a _. (102) 


Using (97) to eliminate the q;, we finally derive as an analogous equation to 
(97) 


P42 = P1P2eXp [— V(ry9)/kT'] {1 a0 |e 3dr, 
+f Sisfeafsat+Sishiafoafsa+ «+++ )p3psdtgdt4 + si . (108) 


Here the general term is the sum of all multiply-connected clusters of order 
t (with the completion of the bond 12) divided by (t—2)!, remembering 
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that systems 1 and 2 are specialized and that the remaining (t—2) can be 
permuted. Thus the terms involving 4— clusters can be written 


l 3 | 3 | 3 | 3 I ) 
; O 
: 1+ A+ Nite fy] +e! 
Dena 2a A aeeiia yey) 24 
(104) 


Equations (98) and (102) are the master equations in Yvon’s theory. 
The application to crystal lattices proceeds as in the previous section, the 
integrals being replaced by sums. We may thus write instead of (98) and 
(102) for two neighbouring lattice sites a and by : 

U 


Inag— Dfaiui-t Laifasfiy +++ =C— FB. . (105) 
v rey) ; 


Ogh = %_Xp (1 Jap) [ 1+ D Sako Ou + mS (Said os Fig 


t+falaglojSig + + wos JOO; + 73 nae Bee ois (106) 
Here f;; is defined by (74), and the sums f,,;, f,;, fi; ....fori,j....4 must be 
taken over all pairs of lattice points which coincide or are nearest neighbours ; 
the detailed enumeration of the terms in (105) follows (99), and in (106) 
follows (104). 

The terms in (105) may now be rearranged as in the previous section into 
groups corresponding to occupation (single or multiple) of a single lattice 
point, occupation (single or multiple) of two nearest-neighbour lattice 
points, occupation (single or multiple) of three nearest-neighbour lattice 
points forming a triangle, and so on. By analogy with (86) we can write 
for the left-hand side of (105). 

J 1(aq > Z) a5 p2 J (oa, Xp; 2) or 2 J; (aq; Xp, Xe 5 z) 
, C 
ee POO apie ge 2) te ele g + Let (107) 
b,e,d 
the first sum being taken over the q nearest neighbours 6 of a, the second 
sum over the other vertices b, c of the 3p, triangles passing through a, and 
so on. The terms J,(a,;2), J2(%g:%p3 2), Ja(%q Xp, M3 2%)... CaN again 
be found explicitly by applying (105) to pseudo-lattices (fig. 8); the 
constant C can be eliminated by considering explicitly the first few terms 
of the expansion (105). We thus readily find from the pseudo-lattice in 


fig. 8a that 
Broo) sees ee LiL or (haa) Joe ae be eb oles) LOS) 


+ We drop the suffix 1 from a to avoid confusion. a, is thus the mean 
fraction of systems on lattice point a in state 1 (a, <}); «4, 1s the mean fraction 
of pairs of systems both in state 1 on lattice points a and 6. 
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To preserve complete generality for the calculation of J, (a,, a, ; 2) we must 
take the external fields U, and U, as different. We then have, for the 


pseudo-lattice in fig. 8b 
og =e(ve2+1)/(l+_t+vt+pre), p=exp(—U,/kT), 
Op = v(me2 + 1)/(l+yt+vt+pve), v =exp(—U,/kT), . (109) 
z=exp (¢€,,/2k7). 
If py is now eliminated, the logarithm taken of y, and the result compared 
with (105), we can deduce that 


Jo (Og, Xp 3 z)= 
(o%4 + Xp — L)2? + (oq — an) + [aa + op — 1)?24 + 2 (oq + Oty = Og — 0147)2? + ecu : 


In 


9 
ay 


The complication in (110) arises since the equation can deal both with 
ferromagnetism and antiferromagnetism; for the ferromagnetic case in 
which «,=«, a considerable simplification results. In fact in the latter 
case equation (105) can be identified with the well known formula in the 
Mayer theory expressing the absolute activity in terms of the cluster 
integrals, and may be written (Rushbrooke and Scoins 1955, eqn. (27)) 
Ing =Ina~qine—)} Biot) 2 eee 
uS1 
Yvon went on to deal with higher-order terms such as triangles and 
squares, but only a maximum of two independent sub-lattices need usually 
be considered, giving two variables a, and a,. Thus J, (ag, a), &, &}; 2) 
required in the next term for loose-packed lattices is obtained by eliminating 
pand v from the equations 


2 Sy 2p? + 2-4(2Qu2v + wv?) + 22-7 uy + wp? + 
 egleen a Er WET SO" VER WRENS PRC CRISEET TST ear RO ee 
* g Sup? + 22-4(u?2y + pv) + 42 uy + p2+ 2+ 2(u+v)4+1’ 
% 2S u?y? + 2A*(u2y + Quy?) + 22 + v2 + v 
 z-8y Py? + 22-4(2y + py?) + 422 + pw? + 2+ Aut v) +1’ 


(112) 


Xp 


using (105) for ~, and «,. Successive approximations correspond to stopping 
at the end of a finite number of terms of (107); the first two terms again 
correspond to the Bethe approximation. 

Equations (105) and (107) are used by Yvon to determine the long-range 
order. ‘Thus for the Ising model in the absence of a magnetic field we must 
look for a solution with a, and a, equal to 4(1+R) where R40; we know 
that a symmetric solution having the same partition function then exists 
with «, and a, equal to $(1—R). Hence if we denote the sum of the 
first ¢ terms of (107) with a,=a,=a...=a by J (a; z), the long-range 
order is given by the equation 


JORIL+R);z]=JHRI-R);2) 2 2... (118) 


(110) 
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Thus for the first two terms giving the Bethe approximation we find from 
(108) and (110) that 


14+ R\@ig-1 — —-14[14+(R-2- 12-22 
ik 1+ [14+ (R122 pe 


Once & has been determined from (112), C can be calculated in (105), and 
hence « determined in the presence of a field. 

For the ordered phase of an antiferromagnetic, and for alternating order 
generally, the two sub-lattices must be treated separately (cf. § 3.6.4), and 
the two functions J!(a,,«,; 2) and J" Vop, 3 2) must be used. In the 
absence of a magnetic field, and when the sub-lattices are equivalent, the 
resulting equation is the same as (113), asexpected. When the sub-lattices 
are not equivalent a different result is obtained (see e.g. Fournet 1952b). 

To determine the short-range order and internal energy, eqn. (106) must 
be used. This can be transformed like (107) into the form 


(114) 


Inog, = Kola gs Xp 3 z)+ ¥ K3(a,, On, Oe 2)+ > Kila, ps ey gs 2) as 
c c,d 


(115) 
where the first term corresponds to occupation (single or multiple) of a pair 
of nearest-neighbour lattice points, the second to a triangle, ... . as before. 


The terms of the right-hand side of (115) can again be calculated by the use 
of pseudo-lattices. Thus Yvon gives, for K,(a,, «, ; 2) the expression 


in {- + (2-8 — 1) (aq + op) — [(&q + op — 1)? + 22-F8( ag + 44 — 1g? — wy”) 274g, — Tt 


2(z-#— 1) 
(116) 
and K 4(0q, %p, %q> %p; 2) is obtained from 
Ko(aq; Xp 5 z) TA alorg, Xp > z) 
282 y? + 2-4(uv? + pe? v) +2 uy 
eB y2y2 4+ Qz-48( uv? + pv) + 44 + pw? +? + 2(u tv) +1? 
(117) 


pand v being derived from (105) and (107). oe ; 

As might perhaps be expected the second approximation using the 
method of Yvon yields results identical with those of Rushbrooke and 
Scoins. However, in Yvon’s method higher approximations can be pursued 
consistently for all multiply-connected configurations. The Curie tem - 
perature is that for which R, given by (113), becomes zero ; hence it satisfies 
the relation 


 (JOYL+R); 2D }n-0=0. pemee eye 118) 


The latter equation can be quite simply expressed in terms of the partition 
functions of the finite clusters which constitute Ja; z), and a series of 
improved estimates of the Curie temperature can be derived (Fournet 
1957). These will be discussed in §§ 4.6, 4.7.1. 
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Although the Yvon technique may seem a little complicated it should be 
emphasized that once the basic formulae have been established they can 
be applied equally to direct or alternating order; that the method is exact 
in principle and a definite procedure is described for obtaining successive 
terms; and that further neighbour interactions can be taken into account 
with little additional effort (see e.g. Fournet 1952a; the latter aspect will 
be referred to in more detail in a subsequent publication). 


4.5.3. High temperature series expansions 


It is a well-known property of power series, all of whose coefficients are 
positive, that the first singularity lies on the real axis, and is hence deter- 
mined by the radius of convergence of the series. | Expansions of the parti- 
tion function (for the Ising model) in zero field, or of the corresponding 
initial susceptibility in inverse powers of the temperature have all their 
coefficients positive, and can therefore be used to estimate the Curie tem- 
perature. The method will be effective if a sufficient number of terms is 
available for a reasonable estimate of the asymptotic behaviour of the 
coefficients to be possible. 

The technique of deriving such expressions is discussed in detail in § 5. 
For the partition function in zero field for three-dimensional lattices we 
find (Domb and Sykes 1957 a) : 


Simple cubic: 

Z = 2(cosh K)3[1 + 3wt + 22w® + 192w8 + 2046w! + 24853wl" + ...] 
Body-centred cubic: 

Z=2(cosh K)4[1 + 12w* + 148w® + 2568w8 + 53944wl + ...] 
Face-centred cubic : 

Z = 2(cosh K)*[1 + 8w? + 33w! + 168w* + 962w® + 5928w? 

+ 38907w® + 268056w* + ...] 
K=J/kT, w=tanhk. 

The first series was obtained by Wakefield (1951) who used it to estimate the 
Curie point of the simple cubic lattice. Wakefield derived from this series 
a corresponding expansion for the logarithm of the partition function, and 
analysed the coefficient c, of w?” in this expansion. Assuming it to be of 
the form n?/w," the best fit was obtained by taking g= —3 (only integral 
values were considered) ; using this result to approximate to the asymptotic 
form of c,, Wakefield estimated the thermodynamic properties of the 
model at temperatures above the Curie point. In fact there were errors 
in the last two coefficients of Wakefield’s series which have subsequently 
been corrected}; but his estimates of the critical point (k7'.=4-50J) and 
thermodynamic functions seem quite reliable. 

Domb and Sykes (1957b) have preferred to use the susceptibility for 
locating the Curie point. This series contains both odd and even terms, and 


(119) 


i ee 
+ The small error in the final term of the series given by Domb and Sykes 
(1957 a) was pointed out by Rushbrooke and Eve (1959). : 
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the coefficients seem to settle down more rapidly to smooth behaviour 
(probably because a greater variety of configurations are needed for its 
derivation). If we write for the initial susceptibility, 


m2 © 
Xo= TD MK”, eR Oh ae eee 120) 


n=0 


Table 2. Coefficients in the High Temperature Susceptibility Expansion 


Xo=2[kT > a, K™ 


n=0 


Body- Face- 
centred centred 
cubic cubic 


Simple 
quadratic 


Simple 


Lattice . 
cubic 


Triangular 


1-00000 1-00000 1-00000 1-00000 1-00000 
1-00000 1-50000 1-50000 2-00000 3-00000 
0-750000 1-87500 1-87500 3°50000 8-25000 
0-541667 2-12500 2-31250 6:08333 21-8750 
0°359375 2-28906 2:75781 10-1979 56-8906 
0-234896 2:39141 3°28203 17-0635 146-205 
0-149219 2:44661 3°84603 28-1373 372-766 
0-0938836 2:46368 4-50157 46-3590 945-047 
0-0581305 2-45098 5:22797 75-7855 2385-72 
0-0357232 6-06678 123-820 


the coefficients a, are given in table 2 for a number of two- and three- 
dimensional lattices. On physical grounds it is reasonable to assume that 


for sufficiently large n 
A, ~n9/K 0”, en ee west sp r( 121) 


where g is a measure of the curvature of the reciprocal susceptibility curve 
near the Curie point, g=0 corresponding to zero curvature. This 
~ assumption can be put in the equivalent form 


Cea incor a VaGatye yen es Ge (122) 


which is more convenient for analysis. Thus if a,,,,/a, is plotted against 
1/n the resulting curve should tend to a straight line for large n, and from 
the position and slope of this line 7’, and g can be estimated. 

In order to provide an unbiased comparison of the properties of various 
lattices it is convenient to introduce the factor g, and hence to estimate 
kT'./gJ. Domb and Sykes thus analyzed the behaviour of 


(ONC este Tels (net. (123) 


The successive values of w,, are plotted against 1/n in fig. 9. For any two 
values of n by taking the intercept of the straight line joining the corres- 
ponding values of u,, on the line 1/n=0 an estimate is obtained for kT'./qJ. 
It is found that if u, and w,,, are used, the loose-packed lattices show a 


P.M.S. U 
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marked distinction between odd and even terms. Instead, therefore, the 
mean of consecutive pairs of intercepts is taken, which corresponds to using 
U,_, and U,,,. The values of these means, or of 


L(+ lca — (8 1) gy be ee 


are reproduced in table 3. It will be seen that the estimates for the two- 
dimensional lattices based on this approach are very close to the true 
values. Estimates for successive n manifest a little irregular scatter, and 
the mean of the last two values is therefore taken as the final estimate. 


0.8 Pee Meare 
On tre 
qon. SS. 
0.6 linear chain - oad 
0.4 
0-2 
0.0 


1.0 0-8 0.6 0.4 0.2 LL 0.0 
n 


Estimation of the Curie point from high temperature susceptibility expansions. 
The values of wp(=dp/qa,_-1) rapidly approach linearity in 1/n, and 
hence the Curie point can be estimated quite accurately. 


Table 3. Values of $[(n + 1)u,,,—(m—1)up_4] 


Body- Face- 
Simple quadratic | Triangular | Simple cubic centred centred 
cubic cubic 
0-58333 0-63333 0-73333 0-80357 0-82576 
0-56792 0-60294 0-75676 0-80137 0-81714 
0-55072 0-60785 0-75014 0-78798 0-81585 
0-57884 0-60990 0-75360 0-79709 0-81580 
0-56802 0-60842 0-74757 0-79175 0-81616 
0-57095 0:60675 0-75328 0-79606 0-81632 
0-56334 — 0-75031 0-79278 — 
Mean of last two: 
0-56714 0-60759 0-75180 0-79442 0-81624 
Exact value: 
0-56730 060682 — — = 


Cooperative Phenomena in Crystals 279 


Once the critical point has been determined, an asymptotic approxi- 
mation to the coefficients in (119) can more readily be established, and the 
thermodynamic functions calculated in a manner similar to Wakefield. 


4.5.4. Low-temperature expansions. Metastable approximations 


The low-temperature expansion for the simple quadratic lattice in the 
absence of a magnetic field has all positive coefficients, 


A=1+24+ 228+ 528+ 142104 447124 1522144 ..... . (125) 


and the methods of the previous section could be used to estimate the 
Curie point. However, this seems to be peculiar to the simple quadratic 
lattice, and for most three-dimensional lattices the coefficients are not 
consistent in sign. 

For the simple cubic lattice they alternate regularly 


A=1 +28 + 3210 — 3212 4 15214 30216 + 101218 — 2612794 .... . (126) 


This series has a radius of convergence of about z=0-57, whereas from the 
high temperature expansion we have seen that the critical point is at about 
z=0-64. Thus the series does not converge right up to the Curie point. 
Wakefield (1951) endeavoured to subtract a term corresponding to the 
spurious singularity z = 0-57 which gives rise to the alternation in signs ; since 
this is the dominant term, and the assessment of the behaviour near the 
Curie point depends on the remainder, the procedure is not very satisfactory. 
The series for the face-centred cubic lattice is 


A=14+ 2124 6222 — 6224 + 8730 4 49732 — 1142344 66236 + 2.4738 4 123240 
4+ 134284... Paar e res Bt - (127) 


This is much more irregular, and it is difficult to see how any sort of extra- 
polation could be made. 

Domb (1952) therefore suggested that the appropriate series to use in this 
region is the low density expansion (§ 3 (143)) 


A=1+ pf,(z) + pfol(z) +... wh(z)+.... Petes oe Les) 
The term in ! in this series corresponds to configurations involving ¢ 
overturned spins. If this series is stopped after the term in p! a physical 
interpretation can be given to the resulting partition function A, Re- 
membering that the partition function for the N spins of the assembly is 
obtained by raising the series (128) to the Nth power, we see that replacing 
A by A, in (128) corresponds to taking clusters of order ¢ or less into account 
correctly, but ignoring all terms which involve clusters of order greater 
thant. Hence this represents an assembly in which clusters are allowed to 
grow only up to a size ¢, and can therefore be interpreted as a metastable 
assembly. This interpretation is analogous to that of Mayer in his theory 
of gas condensation (Mayer and Mayer 1940, p. 312). Because of this 
physical significance of the partial sums A,, we should expect the series ( 128) 
to converge for all u,z<1, although the convergence will only be rapid for 
sufficiently small z or ». We have already noted (§ 3.7) that this property 


U2 
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has been established rigorously by Yang and Lee. It is significant in this 
connection that the series (128) can be used to derive high temperature 
expansion (§ 5.4.3). 

To give an illustration of this interpretation fig. 10 shows the free energy 
of the partial sums as a function of concentration for the body-centred 
cubic lattice taking¢=4and6. It will be seen that these free energy curves 
have the shape characteristic of an assembly with a metastable region, and 
as ¢ increases a steady series of approximations is obtained to the stable 
phase (it should be emphasized that the approximations are applied only 
in the region to the left of the first minimum or to the right of the second 
minimum in fig. 10, ef. §3, fig. 42 and § 3.7.1). 


Fig. 10 


r=4 r=6 


free energy 


stable estimate 


0 0-1 0:2 0:3 0-4 0-5 0-6 0:7 0:8 0-9 1-0 
concentration 


Metastable approximations for the b.c.c. lattice. The approximations are 
applied for concentrations below that corresponding to the minimum 
in free energy. The approximations correspond to a metastable phase. 


The detailed application of the method is discussed by Domb and Sykes 
(1956), where tables are given of polynomials f,(z) for a variety of lattices. 
A partial sum <A, is referred to as a metastable approximation of order f, 
and if an approximation is required to any physical quantity in the stable 
assembly, the corresponding quantity is determined for metastable assemb- 
lies of order ¢, and extrapolated to t= 00. For example, by determining 
the maximum specific heat of a metastable assembly for a given ¢, and 
analysing its variation as a function of t, it is possible to conjecture whether 
the specific heat is infinite in the stable assembly. The method is then 
somewhat analogous to that used by Kramers and Wannier (1941) for the 
simple quadratic lattice before Onsager’s exact work (§ 3.5.1). The known 
exact solutions can serve as a valuable guide to the interpretation of results, 
and to the estimation of the order of magnitude of errorsin extrapolation. A 
useful feature of the method is that quantities can be compared for various 
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lattices at corresponding stages in the approximation, and hence any 
definite trends can readily be identified. 

As a typical example of the application of the method, the specific heat 
maxima of metastable assemblies with different values of t are reproduced in 
Table IV for various lattices. The results are plotted in fig. 11 against 
log,9#, and it is seen that they very nearly satisfy a linear relation. 

As a series of values approximating to the critical point, Domb and 
Sykes used the temperatures of maximum specific heat of the metastable 
assemblies. The results are shown graphically in fig. 12. In order to 
compare the properties of different lattices on the same scale the estimate 
of k7’./qJ is plotted as a function of 1/t. The variation in the estimates for 
different values of ¢ decreases as the coordination number of the lattice 


Fig. 1] 


logio t 


Maximum specific heats of metastable approximations. simple 
quadratic, —~- --— BIMPIS = OUDIC, ors salh ss 5 body-centred cubic. 
The variation of the maximum specific heat with log ¢ is approximately 
linear. 


increases (a similar behaviour occurs in the high temperature estimates, 
fig. 9). The final estimates of the critical point obtained in this manner 
differed slightly (~ 1%) from the high temperature estimates. Recently 
Burley (1959), using the alternative possibility of locating the critical tem- 
perature as the point of disappearance of the spontaneous magnetization, 
secured very good agreement with the high temperature estimates. 

Having located the Curie point, estimates of critical values of thermo- 
dynamic functions can be obtained by extrapolating the corresponding 
metastable approximations. Thusif H,is the metastable approximation of 
order ¢ for the internal energy, and AH,=H,—H,,, a useful asymptotic 
formula for AZ, is A/t”; any successive pair of terms determines a value of h, 
and if the values for different pairs remain fairly constant the extrapolation 
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is a reasonable one. In this manner Domb and Sykes derived estimates 
of the critical values of energy and entropy, and these will be discussed in 
§ 4.7.2, 

The method of metastable approximations depends on the provision of a 
smooth and steadily behaved series of approximations rather than ona good 
single closed-form approximation. It is exact in principle, and makes 
maximum use of available configurational data for the lattices considered. 


Fig. 12 


kT ,/qJ 


1:0 0:8 0-6 0-4 0-2 0 


Estimation of the Curie point by metastable approximations. Deviations 
from the limiting value are less marked for three-dimensional lattices. 


Table 4. Specific heat maxima of metastable assemblies (C,,,,,,,./k) 


Simple 
quadratic 
lattice 


Body- Face- 
centred centred 
cubic lattice | cubic lattice 


Triangular Simple 
lattice cubic lattice 


0-4392 0-4392 0-4392 0-4392 0-4392 
0-6722 0-6634 0-6634 0-6582 0-6526 
0-8168 0-8084 0-7989 0-7885 0-7819 
0-9242 0-9137 0-8948 0-8803 0-8722 
1-0052 0-9953 0-9671 0-9499 0-9409 
1:0721 1-0600 1-0254 1-0059 

1-1269 — 1-0740 ann 

11745 —- 1-1156 — 


it 
2 
3 
4 
5 
6 
7 
8 


4.6. Comparison of Approximate Methods 


In their 1941 paper, Kramers and Wannier undertook a detailed assess- 
ment of the various approximate methods then available by comparing 
series expansions of the partition function of the closed-form approxi- 
mations at high and low temperatures with the corresponding exact 
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expansions. They were concerned with the simple quadratic lattice and 
some examples of their results are as follows. 
High temperatures : 


4a al T- 1009 
Exact Bao] eK KY Ke OO Kee], 
Variation Z=2| 1 sedrRley qe alll 7g Wop a A 
3 45 315, ‘ 
Bethe Z=2(1+K?+}4 K++ ....]. (129) 
Low temperatures : 
Exact Z=z27[1 +244 226+ 5284142104 |. .], 
Variation Z=z[1 +244 226 + 5284 142104... .], 
Bethe Z=2[1 +244 2284 4284 2... 


Hence they were able to show that their own variation method represented 
a substantial improvement over the Bethe approximation. 

Making a similar comparison for the high-temperature expansion of the 
simple cubic lattice following Kurata et al. (1953) we may write : 


Exact Z = 2(cosh K)3[1 + 3w4 + 22w® + 192w8 + ...], 


Kikuchi 3rd approximation 
(cubeaselementary figure) Z=2(cosh K)3[1+ 3w*+ 22wt+189w8 + ...], 


ter Haar and Martin 

(Variation) Z=2(cosh K)3[1 + 3wt + 22w' + 173w8 + ...], 
Kikuchi 2nd approximation 

(square as elementary figure) Z = 2(cosh K)?[1 + 3w* + 18w® + 147w8 + ...], 


Bethe Z=2(coshK)?, (w=tanh K). 
(130) 


A better understanding of the nature of the various approximations is 
provided in a study, suggested by the above authors, of the ‘configuration 
counting’ associated with each approximation; for this purpose it is con- 
venient to use the terminology of §5. The Bethe approximation takes no 
account of closed configurations, and would berigorously correct for a lattice 
of coordination number qin which there are no closed circuits, which may be 
called a ‘Bethe Lattice’ (fig. 13a). Series expansions in terms of lattice 
constants derived in § 5.2.10 reduce to the Bethe approximation if all the 
Pin are put equal to zero; thus all lattice constants for a Bethe lattice 
_ depend only on q. 

In the next approximation both q and p, are taken into account, and the 
remaining lattice constants taken to be zero. This approximation would 
be rigorously correct for a lattice in which there are no closed circuits larger 
than triangles (fig: 13b). It corresponds to the second approximation in 
the development of Rushbrooke and Scoins or Yvon. To find whether a 
given configuration is counted correctly or not in this approximation, all 
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cutting points must first be removed by the methods of §§ 5.2.4, 5.2.6 (but 
not cutting lines)+; the counting is correct if the resulting expression 
depends only on p, andq. For example from § 5 (9), 


9 3 
Poo= 5P3— Psa 5 Ps liek) 


Fig. 13 


(a) © 
(d) 
Pseudo-lattices for which certain approximations are exact. (a) A ‘ Bethe 
Lattice ’ containing no closed circuits. (b) Pseudo-lattices whose only 


closed circuits are triangles for which the 2nd approximation of 
Rushbrooke and Scoins is exact. 


Therefore p,, will not be correctly counted because of the p;, term. By 
carrying out the reduction of simply-connected configurations, and 
putting all lattice constants other than p, and qg equal to zero, the series 
expansions in § 5.2.10 can be applied to this approximation. The generali- 
zation of configuration counting to higher-order approximations of this 
kind is immediate. 


} The Mayer development on which these approximations depend removes 
all simply-connected configurations, but not configurations which are reducible 
in the sense of § 5. 
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The Kikuchi approximation counts correctly all configurations which are 
open with respect to the basic figure (§4.4.3). To analyse this approxi- 
mation all reducible configurations of §5 must be expressed in terms of 
irreducible lattice constants; if these expressions involve only configura - 
tions contained in the basic figure, the counting will be correct. Thus from 
§5(11) both p;, and pg, will be correctly counted by the Kikuchi approxi- 
mation using a triangle as basic figure, and hence this approximation is 
better than the corresponding approximation of the previous paragraph. 
Similarly for a loose packed lattice with a square as basic figure p,, and p,, 
will be correctly counted. However, when we look at simple polygons, 
the individual space types must be examined. Thus considering hexagons 
on the simple cubic lattice, the first two space types in §5, fig. 6, will be 
counted, but not the last one; hence the coefficient 18 instead of 22 of w® in 
(130) for this approximation. When a cube is taken as basic figure all 
hexagons are correctly counted; but when we examine octagons in §5, 
fig. 7, all space types except (i) are correctly counted, and the coefficient of 
w® is 189 instead of 192. 

The approximations of §$4.5.3, 4.5.4 make direct use of exact series 
expansions as far as they are available, and depend on the smooth con- 
vergence of partialsums. For the Ising model of a ferromagnetic they seem 
particularly suitable; for the antiferromagnetic in a magnetic field (§ 3.7), 
where terms alternate in sign, they cannot be readily applied. It was noted 
by Domb and Sykes (1957 b) that high and low-temperature estimates of the 
Curie point differed slightly. The properties of the disordered state are 
simpler than those of the ordered state, and for most approximations the 
Curie point can be more readily located from the disordered state. It is 
likely that this is also true of series expansions, and thus high temperature 
estimates are probably more reliable (this is also the conclusion of Burley 
(1959)). This would mean that the curve of metastable estimates (fig. 2) 
passes through a minimum and turns up slightly in the limit. Dr. M. F. Sykes 
has suggested (privately) that the use of the series for In A instead of A 
in (128), corresponding more precisely to the Mayer theory, might lead to a 
more uniform approach to the critical point. 

The reliability of the closed form approximations would be substantially 
increased if they could be made to supply a steadily behaved series of 
values of any thermodynamic quantity to be estimated. The initial im- 
provement achieved by the Kikuchi approximation is striking, but it is 
difficult to proceed to high approximations. The method of Rushbrooke 
and Scoins as formulated at present can only make use of homogeneous 
configurations, and hence it is not possible to provide systematic improve- 
ments. The method of Yvon can take all configurations into account, and 
the formulation of a series of successive approximations to the critical 
point was undertaken by Fournet (1957). However, only three terms of 
this series were calculated, and although the results of extrapolation seem 
to provide better estimates than any other single closed form approxima- 
tion, we do not think that they are as reliable as the results of direct series 
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expansions in §§ 4.5.3, 4.5.4, which make use of several more terms. The 
results of Fournet might well be extended with the aid of the configurational 


data of § 5. 


4.7. Physical Properties of the Three-dimensional Ising Model 


Before discussing in detail the characteristic features of the three-dimen- 
sional Ising model we shall briefly summarize the main conclusions to which 
we have been led. The specific heat of the three-dimensional model seems 
to remain infinite at the Curie point, but its detailed shape is appreciably 
different from the Onsager curve in two dimensions. The ‘ tail’ is far 
smaller, showing that most of entropy change of the assembly takes place 
in the temperature region below the Curie point ; this is in far better agree- 
ment with experimental observations of specific heat anomalies. 

The dimension of the lattice is of primary significance in determining the 
properties of a model, and the coordination number plays only a small part. 
The differences between lattices of the same dimension are thus only of 
secondary importance. As an example of this, the two-dimensional tri- 
angular lattice and the three-dimensional simple cubic lattice both have 
coordination number 6; but the properties of the former are very similar to 
the honeycomb lattice of coordination number 3, and those of the latter to 
the face-centred cubic lattice of coordination number 12. 

Simple closed form approximations (such as the mean field and Bethe 
approximations) are closer to the true solution in three dimensions than in 
two dimensions. Hence the inverse susceptibility above the Curie point 
is more nearly linear in three dimensions than in two dimensions, and the 
spontaneous magnetization curve at the Curie point is less steep. 


4.7.1. Position of the Curie point 


It is interesting to compare the values of k7'./q¢J for various lattices as 
estimated by different approximations (table 5); this quantity is con- 
siderably more sensitive than z_ [=exp—(J/k7'<)] for larger coordination 
numbers. For the mean field approximation it is 1 for alllattices. For the 
Bethe approximation, (28), it depends on q only, being zero for g=2 and 
increasing steadily to the mean field value as g+oo. The Kikuchi 
approximation takes into account the first few closed configurations, and a 
marked separation becomes manifest between the triangular and simple 
cubic lattices. 

The first four approximations which are of closed form all overestimate the 
Curie temperature. The remainder depend on some form of extrapolation, 
and may either overestimate or underestimate it. Fournet’s results, 
based on the method of Yvon, use the first three stages of the approximation, 
which are fitted for t=1, 2, 3 to a second degree polynomial in 1/t; this is 
then extrapolated to 1/t=0. The value for the simple cubic lattice seems 
appreciably better than those for the body-centred cubic and face-centred 
cubic lattices. ‘The estimates of the metastable approximations and high- 
temperature series expansions are quite close to one another, and for the 
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reasons stated in the previous section we have adopted the latter as final 
estimates. Polynomial fitting in 1/f was not considered to be the most 
suitable method of extrapolation since it pays too much attention to slight 
irregularities. Instead some suitable asymptotic formula was postulated, 
the parameters of which were fitted from the results of successive approxi- 
mations (§§ 4.5.3, 4.5.4). 


Table 5. Approximations to the Curie Point 
(Recorded values are te =kT'¢/qJ) 


Lattice Simple ; Simpl DOr ain ee: 
approximation Sriileatis paengalar capes centred centred 
cubic cubic 
Mean field 1 1 1 1 1 
Bethe and quasi- 
chemical 0-721 0°822 0-822 0-869 0-914 
Rushbrooke and Scoins | 0-657 0-744 0-794 | 0-835 | 0-882 
Kikuchi (2nd 
approximation) 0-606 0-653 0-768 0-869 0-835 
Yvon—Fournet 
(extrapolation) 0-578 0-622 0-761 0-865 0-704 
Metastable 
approximations — — 0-748 0-784 0-805 
High temperature 
susceptibility series 0-567 0-608 0-752 0-794 0-816 
Exact 0-567 0-607 — == = 


4.7.2. Critical values of thermodynamic functions 


The estimates of critical values of energy and entropy provided by a 
given approximation are less accurate than corresponding estimates of the 
critical temperature, since a double process of approximation is involved. 
(Some typical results are compared in table 6.) A clear improvement can 
again be noted in passing from a two-dimensional to a three-dimensional 
model. 

The metastable estimates have been evaluated at the values of the critical 
point given by the high temperature expansions (table 5). Critical values 
of A and £ are determined by noting the size of successive approximations 
for the known values of ¢ and fitting the difference between the tth and 
(t—1)th terms to A/t*; any successive pair of terms determines a value of h 
and if the values for different pairs remain fairly constant it is reasonable to 
extrapolate to t= co. (Further details are given in Domb and Sykes 
(1956)). The results are in good agreement with independent estimates 
based on high temperature series expansions. 

We have collected in table 7 exact critical values for two-dimensional 
lattices and final estimates for three-dimensional lattices. Im order to 
compare the difference between the specific heat curves for the various 
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lattices, it is convenient to take 7’; as a unit of temperature, and consider C,, 
as a function of 7(=7'/7'-). Then 


(132) 


Table 6. Approximations to Critical Values 


Simple Simple Face-centred 
Lattice quadratic cubic cubic 
approximation 


— ol S elk —E old S elk —E Ac S elk 


Mean field 0 0-693 0 0-693 0 0-693 
Bethe or quasi-chemical | 0-667 | 0-580 | 0-600 | 0°633 | 0°545 | 0-668 
Kikuchi 

(2nd approximation) 1:123 | 0-428 | 0-857 | 0-587 “2 0-618 
Metastable 

approximations 1-405 | 0-311 | 1-015 | 0-553 0-589 
High temperature series 

extrapolations 1-406 | 0-311 | 1-071 | 0-541 0-597 
Exact 1-414 | 0-306 — = — 


Table 7. 


Lattice structure | q | kT ,/gJ(=t,)| S,/k |(S.—S,)/k| (EZ .—Eo)/kT .|\(Eo—E ,)/kT, 


linear chain é 0 0 0-693 0 oe) 

honeycomb 0-506 0-265 0-428 0-227 0-761 
simple quadratic 0-567 0-306 0-387 0-258 0-623 
triangular 0-607 0-330 0-363 0-275 0-549 


simple cubic 0-752 0-541 0-152 0-428 0-237 
body-centred 

cubic 0-794 0-586 0-107 0-460 0-169 | 
face-centred 

cubic é 0-816 0-597 0-096 0-463 0-150 


so that a tabulation of S, and S,,—S, for various lattices enables us to 
compare the magnitude of the specific heat curves below and above the 
Curie point. The sum of the two terms, S,,, is the same for all lattices, 
and is equal to kin 2(=0-6931k), For a linear chain the Curie point is at 
7’ =0, and all of the specific heat is above the Curie point. As the coordina- 
tion number increases, more and more of the specific heat curve moves below 
the Curie point, until for ¢= oo the mean field approximation shows us that 
all of the specific heat curve is below the Curie point. 
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Critical values of the entropy and the total entropy change in the 
assembly are particularly useful quantities for testing the validity of a 
given model; direct comparison with experiment is possible since these 
values do not depend on the interaction energy. 

The difference in passing from a two to a three-dimensional lattice is 
quite striking. For all the two-dimensional lattices more than 50°% of the 
entropy change occurs above the Curie point, but for the simple cubic 
lattice the ‘tail’ is very much reduced, even though this lattice has the 
same coordination number as the triangular lattice. For the body-centred 
cubic and face-centred cubic lattices the ‘ tail’ is even smaller, and this must 
lead to a specific heat curve in closer agreement with experimental curves. 

The difference between two- and three-dimensional assemblies is even 
more evident if we look at the values of 


1 
(B,— Bok P= 7 | Oar, | 
0 


ae (133) 
EH je? = 4 | Gian 
Bae 


The sum of the two terms, — H,/k7',, is not constant, but decreases from oo 
to a limiting value of 1/2 as goo. These terms represent directly the 
areas under the specific heat curve below and above the Curie point, and the 
decrease in the ‘tail’ is clearly demonstrated as the coordination number 
increases. 


4.7.3. The specific heat curve 


In the previous section we have made use of critical values to draw some 
general conclusions regarding the dependence of the specific heat curve on 
lattice structure. We shall now attempt a more detailed description of its 
behaviour. 

At sufficiently low temperatures the specific heat is given by the first 
term of the low-temperature expansion and we have 


Cj (Fir) exp(—2aFik), ig ch tera OUP 


kT 
Taking a reduced temperature, 7, as in the previous section, we may write 
this in the form 2\2 
C lnex(=) exp (—2/7t,). perry the (130) 
Tue 


Hence, for a given 7, C,, increases as t, increases, i.e. as we pass down the 


rows in table 7. 
Correspondingly, at sufficiently high temperatures the specific heat is 


given by awed 

Ci[k= 3 (am) kil SCIEN Bias tay eae ere @ LD) 
or, in reduced units, by 

CO, [k= 1/2q77t ” hee ee weit lot) 
Hence, for a given 7, C,, decreases as q or te increase, and this supports the 
conclusions of the previous section. 


290 C. Domb on 


We have seen in § 4.5.4 that metastable approximations indicate that the 
specific heat remains infinite in three dimensions. For the purpose of 
detailed comparison between two and three-dimensional lattices Domb 
and Sykes (1957 c) used high temperature series expansions for the triangu- 
lar and face-centred cubic lattices; odd and even terms occur in the 


Fig. 14 


Cy/k 


1 
ta 
(a) 
Comparison of specific heat curves given by various approximations. 
(a) Face-centred cubic lattice. ————— Based on series expansions. 
—----- Kikuchi approximation. ...... First-order Beth approximation. 


expansions for these close-packed lattices, and hence more terms are avail- 
able for analysis. The expansions in reduced form are 


Triangular : 
C,[k = 0-2263 x? (1 + 1-0986 « + 0°8298 x2 + 0-5525 «3 + 0-4021 4 
+ 0°3384 «> + 0-2971 «® + 0-2614«7+4 ....). 


Face-centred cubic: (138) 


C,, [k= 0-06257 x2(1 + 0:8169 « + 0:6778 x2 + 0-5680 x3 + 0-4818 x4 
+0-4190 «+ 0°3725 «6+ 0-3364«7+ 2...) («= 1/7). 
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In comparing these two expansions we may say roughly that the term 
outside the parentheses represents the magnitude of the ‘tail’ of the 
specific heat curve, and the term inside the parentheses represents the shape 
of the curve, particularly near the Curie point. It will be seen that the 
magnitude of the tail is much smaller for the face-centred cubic lattice. 


Fig. 14 (continued) 


(0) 
Comparison of specific heat curves given by various approximations. 
(b) (After Onsager) simple quadratic lattice. _——————_ Exact. 
------- Kramers and Wannier variation (= Kikuchi). ...... Bethe. 


But it will also be noted that the coefficients decrease less rapidly for the 
face-centred cubic lattice and the specific heat curve is therefore sharper 
(and is thus certainly infinite). By analysing the magnitude of the nth 
coefficient in the face-centred cubic expansion Domb and Sykes obtained 
indications of an asymptotic behaviour of the specific heat of the form 
C[k~(1—x)-1? (624). 
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Fig. 15 


(6) 


Estimates of the spontaneous magnetization for three-dimensional lattices. 


(a) Comparison of approximations for the simple cubic lattices, 


Metastable, - --—--—--— KikuGhines one cree Bethe. (b) Comparison 
of spontaneous magnetization curves for various lattices; t—————— 
simple cubic, — — — —- - body-centred cubic,......... face-centred 


cubic. (Two-dimensional simple quadratic for comparison — + — - =H 
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On the low temperature side the estimation of the specific heat is less 
reliable, and metastable approximations have been used. The final 
estimate of the specific heat curve for the face-centred cubic lattice is 
reproduced in fig. 14 (a) together with estimates of other approximations. 
A similar comparison for the simple quadratic lattice was given in the 
classic paper by Onsager (1944) and is shown in fig. 14(b). The improvement 


in the approximations in going from two to three dimensions will again be 
noted. 


4.7.4. The spontaneous magnetization 


The dependence of the spontaneous magnetization on lattice structure 
has recently been carefully investigated by Burley (1959). Near the Curie 
point the mean field and Bethe approximations both lead to the conclusion 
I,~A(1—T/T.)#* (see (9) and (31)). However, for two dimensional 
lattices the exact solution (§§3.5.1, 3.5.4) gives J)~(1—T7'/T.)"8, so that 
the approximations are seriously in error. Burley used the Kikuchi 
approximation and the metastable method to estimate the spontaneous 
magnetization for three-dimensional lattices. The difference between 
these two estimates was very small (fig. 15a). He found that the curve of 
spontaneous magnetization approaches the axis at 7’, less steeply in three 
dimensions than in two dimensions, the power of (1 — 7'/7'-) being estimated 
as between 1/2 and 1/4; also the standard approximations are more accurate 
in three than in two dimensions. These results are illustrated in figs. 15a, 
b. Thus we find again that the dimension of a model is the major factor in 
determining the shape of the spontaneous magnetization curve, the 
difference in crystal lattice leading only to a ‘fine structure’ variation. 

The conclusion regarding the steepness of the curve also results from 
an alternative method of using series expansions due to Park (1956) who 
attempted to approximate to the known terms of the expansion by a 
closed algebraic expression. Park’s estimate of the power of (1—7'/T,) 
for the spontaneous magnetization in three dimensions was initially 17/64, 
but was subsequently modified to 3/16 (private communication). 


4.7.5. The high-temperature susceptibility 


The initial susceptibility above the Curie point will be discussed in detail 
in a subsequent publication. We here quote the result that the inverse 
susceptibility is less curved in the neighbourhood of the Curie temperature 
for a three-dimensional model than for a two-dimensional model (and hence 
closer to the linear result of the mean field approximation (11)). In two 
dimensions y,~ A(1—T7'-/T)~74 near the Curie point (this has recently 
received rigorous justification by Fisher (1959)); in three dimensions the 
corresponding estimate is yy~ A(1 —7'¢/7’)~** (Domb and Sykes 19576). 

It has been pointed out by Rice (1954) that the infinity in specific heat 
corresponds only to an incompressible lattice ; if the lattice is compressible 
a first-order transition results (see also Domb (1956)). 


P.M.S. x 
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Notation for §4 
H, parameter in mean field approximation: =H) kT. 
§=(N,—N,)/N =a,— 4. 
H,=parameter in Bethe approximation; p,=exp(—2#H 1/kT). 
H ,, 4, interaction energy and partition function of an elementary cluster 
consisting of a central spin and its g nearest neighbours. 
Ly, Le} Yrs Yas Y3 3 21» Za» 23, Za» M5» 2p i Wis We, We, Wa, We, Wes 


parameters in Kikuchi approximation 
L number of assemblies in ensemble for Kikuchi approximation 


2 3 6 6 
Ay = Lee IL“) L_2,= Ihe) lL W,= II (2) ! 
iad fai in jun 


Cluster Sum Theory 
f=2-1, fj; defined by (74), b,(z) by (80), B,,(z) by (83) H,2(04, 2) by (86). 


Method of Yvon 


%q, %, mean fraction of systems on lattice points a, bin state 1. a, mean 
fraction of pairs of systems on lattice points a, b both in state 1. 


J 3(%q3 %)> J o(X%qs %p3 2), Sg (aq, Xp, %e3 Z).... defined by (107); J#(a; z) see 
(113). 


Ky(a%q; &p 3 Z), Ke(aq, Xp; 2) defined by (115). 


U., U, external field at lattice points a, 6; p=exp(—U,/kT), 
v=exp(—U,/kT). 


a, coefficient of K” is series for susceptibility (120), a, ~n9/K.”. 

Un = 2p |Q4y_1: 

A,=sum of terms of partition function up to and including term in pl. 
Tr=2/T ,: «=1/7: 
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§ 5. CONFIGURATIONAL PROBLEMS AND SERIES ExPaNsrIonsy{ 
5.1. General Introduction 


Perturbation expansions have been used widely in astronomy and 
physics to calculate the effect of small changes in problems for which exact 
solutions are available. However, the derivation of substantial numbers 
of terms of a perturbation expansion, and the use of such an expansion as a 
source of exact information have been more restricted. A significant 
example of the latter kind is the Mayer theory of the condensation of gases, 
and recent developments in quantum electrodynamics provide a striking 
additionalexample. In classifying the contributions of higher order terms 
one is naturally led to a consideration of the various types of linear graph 
which can be constructed from a given number of points and lines. Some 
general references on the theory of linear graphs are given at the end of 
the chapter. 

Fig. 1 


BaGligee 


ADRVAS 


(0) 
Classification of linear graphs as ‘ separated’ or ‘ connected’. (a) A separated 
graph of 12 lines; (6) A connected graph of 12 lines. 


The application of linear graphs to the Mayer theory has been particularly 
studied by Uhlenbeck and his collaborators. In this theory each 
topological type of graph leads to a particular cluster integral, and Riddell 
and Uhlenbeck (1953) solved the combinatorial problem of determining the 
number of graphs of each type. The problem of the evaluation of the 
resulting integrals for any reasonable physical potential remains extremely 
formidable. 

We shall classify the graphs which can be formed with p points and 1 
lines into groups. We first distinguish between 

(A) Separated graphs which consist of two or more graphs which are not 
connected together (fig. 1a) and 

(B) Connected graphs (fig. 16). 

We now sub-divide (b) into 
(i) ‘Stars ’—multiply-connected graphs (fig. 2a). 
(ii) ‘Trees ’—graphs which are not multiply-connected (fig. 2). 


} This chapter has been written in collaboration with Dr. M. F. Sykes. 
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By multiply-connected we mean that at least two independent paths 
which cross only at the end points connect every pair of points of the 
graph. The graphs of group (ii), which are not multiply-connected, must 
contain at least one articulation point or cutting point, which is a point 
where the graph could be cut into two or more separated graphs by cutting 
all lines going to this point. This group may be further sub-classified into 
(a) Cayley trees (Cayley 1889-98) which are connected graphs of p points 
and (p — 1) lines containing no closed circuits (fig. 3 a). 

(6) Husimi trees (Husimi 1950) which are built up similarly out of simple 
polygons so that no line lies on more than one closed circuit (fig. 3). 
A pure Husimi tree is one which contains only one type of polygon. 

(c) Star trees. The remaining graphs with articulation points in which 

the units are more complex than simple polygons (fig. 3c). 


Fig. 2 
(2) 
(0) 
Sub-classification of connected graphs. (a) Multiply-connected graphs or 


‘stars’; (b) Graphs which are not multiply-connected and therefore 
possess a cutting point (‘ trees ’). 


For the Mayer theory all the integrals which arise can be reduced to 
integrals of stars; some of these (e.g. simple polygons) can be evaluated 
easily, but in general the evaluation is not possible beyond three or four 
terms. The problems with which we have been concerned are confined to 
crystal lattices, and the cluster integrals are replaced by cluster sums. If 
interactions of all neighbours were taken into account, these sums would 
be almost as difficult to evaluate as the cluster integrals; the simplification 
in the lattice model arises largely since it is a reasonable physical approxi- 
mation to take account only of nearest-neighbour interactions. The basic 
problem is then to count the number of ways in which particular types of 
configuration can occur on a given lattice (using only nearest-neighbour 
links) and this combinatorial problem replaces the evaluation of cluster 
integrals. Also for such lattices many types of cluster configuration 
cannot occur, and hence the expansion can be followed for a substantial 
number of additional terms. 
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It is convenient to consider separately the topological problems arising 
from low-temperature and high-temperature expansions. For the Ising 
model the low-temperature (or low-density) expansions can be put into a 
form exactly analogous to the Mayer theory; the high-temperature series 
can be deduced from them, since the two expansions are not independent. 

However, it is important to develop direct methods of evaluating high- 
temperature series for application to more general forms of interaction in 
crystals, and even for the Ising model many more terms can be evaluated 
by these direct methods for particular problems. 


<< ND 
Phy Bp 
Wee 


(c) 
Sub-classification of trees. (a) Cayley trees which involve no closed circuits; 
(6) 1, Pure Husimi trees which involve only one type of polygon; (b) 
2, Mixed Husimi trees which involve lines and simple polygons; (c) 
More general trees. 


5.1.1. High-temperature configurations 


The form of the high-temperature expansions has been established for 
the Ising model in § 3.6.2. The topological problem is one of ‘line counts’, 
i.e. of classifying and enumerating all configurations of J lines which can 
occur on the lattice. The configurations may be sub-divided as in the 
previous section into connected and separated configurations; but whilst 
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in the Mayer theory the separated configurations are readily eliminated 
for crystal configurations they are of considerable significance and are 
difficult to enumerate. 

For the Ising model in the absence of a magnetic field the configurations 
of interest are closed configurations all of whose vertices are even (i.e. are 
the meeting point of an even number of lines) (fig. 4a). For the susceptibility 
in zero field the corresponding configurations are those with exactly 
two odd vertices (fig. 46). Other line configurations of particular interest 
are simple closed polygons », (fig. 4c) and simple chains c, (fig. 4d). 


Fig. 4 


(a) 


(d) 
Typical high-temperature configurations. (a) No-field configurations for the 
Ising model; (b) Susceptibility configurations for the Ising model; 
(c) Polygons; (d) Simple chains. 


We shall be mainly concerned with lattices in which all points and 
bonds are equivalent, which we shall call regular lattices. It will be 
assumed that the lattices are connected cyclically so that end effects are 
eliminated. The number of configurations of a given type which can be 
arranged on a lattice will be referred to as the lattice constant for the 
particular type. More specifically we shall define lattice constants per 
site so that if p, is the lattice constant for triangles, then the complete 
lattice contains Np, triangles. We shall describe as a closed configuration 
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a connected graph in which each vertex is the meet of at least two lines. 
Closed configurations of / lines other than simple polygons will be denoted 
by suffixes Pq, Py: Pe. Separated configurations will be denoted by 
brackets containing the connected portions, e.g. [ps, p,] for a separated 
triangle and quadrilateral, and [p5, p3, ps] for three separated triangles. 
All lattice constants for regular lattices can be expressed algebraically 
in terms of a selected number of ‘irreducible’ constants. However, the 
concept of irreducibility for lattice configurations differs from that used 
in the Mayer theory. For example, configurations consisting of two simple 
polygons with one side in common (¢.g. P5q, Peo, Appendix II, table A) 
are reducible in our sense but not in the Mayer theory; for this reason we 
have preferred the term ‘multiply connected’ for such configurations. 
Our task will be to evaluate the irreducible constants for different lattices, 
and to show how to express all other lattice constants in terms of them. 


Fig. 5 


A a 


Spin counts. Typical low-temperature configurations arising from 
four points (p=4). 


5.1.2. Low temperature configurations 


Series expansions at low temperatures are generally more difficult to 
derive or use than at high temperatures. For the Ising model the form 
of the expansions has been discussed in § 3.6.1 and the topological problem 
is one of ‘spin counts ’, i.e. of classifying and enumerating the configurations 
arising from a set of points of the lattice (fig. 5). We can define lattice 
constants per site as before, but we now use capital letters (e.g. P;,) to 
distinguish them from high-temperature lattice constants. We still find 
it convenient to classify low-temperature configurations by the number 
of lines which they contain. The two sets of lattice constants are related, 
but the relations are not always simple to derive. Thus P, is the set of 
all quadrilaterals none of whose opposite corners are nearest neighbours ; 
P, is contained within p,, but its complete determination will also involve 
Psa 2nd Peg. In general the line count of a given linear graph is therefore 
larger than the corresponding spin count. It will again be our task to 
show how to calculate low-temperature constants for various lattices. 
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5.2. High-temperature Lattice Constants 
5.2.1. The method of space types 


The first few lattice constants can readily be determined for all lattices 
by elementary geometry. Let us consider, for example, simple closed 
polygons on the simple cubic lattice. The first non-zero closed polygon 
is the quadrilateral, and it is clear that this lattice admits only of 
quadrilaterals that are plane squares. Through each site of the lattice 
there pass twelve squares, four in each coordinate plane. The number of 
distinct squares for a lattice of N points is thus 3N; hence for this lattice 
p,=3. When we come to examine the number of hexagons we find that 
there are three possible space types (fig. 6), and that these can be placed 
on the lattice in 12N, 6N, 4N ways respectively; hence p,= 22. 


Fig. 6 


Hexagons on the simple cubic lattice. (a) 12N, each pair of opposite edges 
determines two configurations; (b) 6N, there are two independent poly- 
gons in each coordinate plane; (c) 4N, each pair of opposite vertices 
determines one configuration. 


It is clear that great care and effort is necessary to obtain higher terms. 
The method of procedure was completely systematized by Wakefield 
(1952); to obtain any new term all the space types derived in previous 
terms must be expanded in all possible ways to include the extra bonds, 
and the new space types must be classified for symmetry. In fig. 7 we 
have further illustrated the method by listing the 11 types of octagon for 
the simple cubic lattice, and from this we deduce that p,= 207. However, 
the number of space types increases rapidly, and there are 73 types of 
decagon, and 756 types of dodecagon. Itis not surprising that Wakefield’s 
calculation of the latter contained slight errors, but these have subsequently 
been corrected (Domb and Sykes 1957, Rushbrooke and Eve 1959); the 
correct values are ~4) = 2412, pj,= 31754. 

The classification of polygons on a lattice into various space types provides 
complete topological information regarding these polygons, and enables 
one to determine useful additional configurational data. For example 
from the hexagons in fig. 6 we can deduce that p,,=18, and for low- 
temperature configurations P,=4, P,,=18. However, it is difficult to 
ensure that no space types have been overlooked at each stage, and the 
determination of the symmetry properties of any space type is not at all 
simple. Also for a lattice like the face-centred cubic the method is very 
cumbersome, and the polygons are much more difficult to visualize and 


draw. 
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We shall therefore describe an alternative method of counting polygons, 
which can be used as an independent approach to the problem, or as a check 
on the enumeration of space types. However, we must emphasize that 
the classification of the first few polygons on a lattice into space types 1s 
of great use; for evaluating certain types of irreducible constant the 
space-type technique is the only one which has so far proved applicable 


(§ 5.2.5). 


Fig. 7 


Octagons on the simple cubic lattice (after Wakefield). (11 types, 
total contribution 207, pg=207.) 


5.2.2. Use of the random walk problem. Path factors 


Let us consider random walks on a given lattice starting at any point of 
the lattice, and returning to the starting point after | steps. We can 
characterize each walk graphically by means of the configuration which 
the walker describes; thus for the simple cubic lattice the first possible 
return occurs after two steps by reversing the direction of the initial step, 
and this can be represented by After four steps there are three 
possibilities which can be represented by 


= /\,LJ 


(A heavy line is described in both directions during the walk; a weak line 
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is described in only one direction). Hence the total number of such 
random walks, 7,, can be split up into contributions from various 
configurations, and we can write for this lattice 


r= 6=2|—] 


(3) 
(1) 
nevmsfa} +4 [A] +8[0 1]. 
(3) Ch5) (3) 


The numbers in brackets under each configuration represent the values 
of the corresponding lattice constants. The numbers before the square 
brackets represent the number of ways in which the configuration can be 
described by a random walker who returns to his starting place after 
L steps. To evaluate this combinatorial factor for a given configuration 
we must sum the number of returns over all vertices of the configuration, 
since any such vertex can be the lattice point from which the walk starts. 
The combinatorial factor is a property of the configuration, and does not 
depend on the lattice considered. 

Expansions such as (1) are useful for calculating properties of matrices 
associated with crystal lattices; for example they have been used in 
connection with the vibration spectra of disordered crystals (Domb et al. 
1959). Explicit formulae can be derived for many of the combinatorial 
factors (some results are given in Appendix I), and 7, can readily be found 
by the use of generating functions (Appendix IT). Hence relations such 
as (1) can be used as a check on the evaluation of the lattice constants. 

For higher terms where the number of component configurations increases 
rapidly it is convenient to enumerate only silhouettes of a given configura- 
tion, and not differentiate between the number of times any line is 


traversed. Thus in r, the silhouette [ ER | 


will correspond to terms 


S34 


in our previous notation; and | ee 


corresponds to | | (ae el le 2 | 
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Appendix I contains the analysis in this notation, of 72, 74, rg and rg for 
the simple quadratic and simple cubic lattices. The coefficients outside 
the silhouettes will now be referred to as path factors. A path factor of 
order / for a given configuration is thus the total number of returns to the 
origin in a random walk of I steps on the configuration which traverses each 
line of the configuration at least once; this number is to be summed taking 
each vertex of the configuration in turn as origin. 

We have extended the analysis of Appendix I to r,) for the simple 
quadratic, simple cubic and body-centred cubic lattices but since some 50 
configurations are involved we have not reproduced the results here. In 
fact, although this analysis is of importance for matrix problems connected 
with crystal lattices, it is somewhat cumbersome as a check of lattice 
constants because of the large number of configurations which enter. 
We shall now proceed to a method in which all configurations which are 
not closed can be removed. 


5.2.3. Elimination of immediate reversals 


The problem of a random walk on a lattice in which no immediate 
reversals are allowed has been considered by several authors [ (Gillis (1955), 
Temperley (1956), Domb and Fisher (1958)]; in such a walk two successive 


steps = cannot occur, although steps like prea t can occur. 


Domb and Fisher gave a simple formula relating g,, the number of returns 
to the origin in a no-reversal walk on a regular lattice, to 7, the number of 
returns in a standard walk : 


(o—1) 


etree ive (yi ano 


BY) "l-4 


l—4)(l—5 
— of CO at lcaqg=l)2 > 2 
This formula applies to walks in which ‘improper closures’ are excluded, 
i.e. walks which start and end in opposite directions, e.g. 


ox} 


From eqn. (2) values of g, can readily be calculated from r,, and some results 
are tabulated in Appendix II; for the simple cubic lattice I= 9, 
da = 24, Yo = 264, gg = 4584. 

If we now analyse the q, into contributions from various configurations 
in the same manner as 7, in Appendix I, only closed configurations can 
enter into the analysis, and this leads to a substantial simplification. 
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Thus for q, and q, only simple closed polygons enter, and we may write 


nestes [J] 


(22) 
Hence we have verified the result of §5.2.1 that p,=22 without any 


geometrical considerations. When we come to consider q, two more closed 
configurations enter our analysis 


TO 


We shall denote each configuration by the corresponding lattice constant 
symbol. (See table A, Appendix IIT), and for gg we may write 


9g =Ag(8) Pg + Aza(8) Pra +Agn(8)Pgn + Ara(8)P- et ae ie 2) 


Here the \’s may again be referred to as path factors but they differ from the 
corresponding path factors of §5.2.2 for they represent the number of 
returns to the origin summed over all vertices in a random walk with no 
immediate reversals (and no improper closures); again each line of the 
configuration must be described at least once in the walk. Quite generally 


we may write 
DTA) Pees Ree ees miaeie a tarsus) 
#5 8<l 
where the variable / in the )’s tells us that the walk consists of / steps. 

It is clear that A,(8) = 16, A,(8) = 8 (successive squares must be described 
in the same sense since immediate reversals are not allowed). A,,(8) is a 
little more tricky. From fig. 8 it will be seen that starting from the 
vertices A, C, D, F the configuration can be described in two ways (arrows 
as shown or completely reversed); starting from B, E it can be described 
in four ways; hence ,,(8)=16. Similarly P,, can be described in four 
ways starting from A, B, C, D, E, F (fig. 9 arrows as shown or completely 
reversed), and in eight ways starting from G (an additional factor 2 arising 
from the initial choice of squares); hence A,,,(8)=32. Thus we have the 
relation 

Gg = 4584 = 16p,+ 32pg,+16p7,+8P4, + - - - (8) 
and this can be used as a check on the values of 3, P7q, Pan. But in fact 
Prq and pg, are reducible configurations, and we shall show in the next 
section how they can be expressed in terms of p,; eqn. (6) can then be 
used as a means of evaluating pg. 
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The analysis of q, is particularly useful for the face-centred cubic lattice 
in which polygons occur with odd and even numbers of sides, and con- 
figurations are difficult to visualize geometrically. We have for this 
lattice g,=0, q3=48, qy= 264, ¢5= 1680, g,=17016. 

We can immediately deduce p;, p, and p; from 

Dine 48 = 6D, 
Qa= 264= 8p, Gale eet) 
ds = 1680 = 10p;, 
and by analogy with (6) we find 
dg = 17016 = 12pg + 24p,6, + 12p5q + D3. 1 a ey 
Fig. 8 


e D 


When we extend the analysis of q, to higher terms we find that substantial 
numbers of other lattice constants enter into relations such as (6) and (8), 
and the determination of the path factors becomes quite involved. We 
will therefore now consider methods of evaluating other lattice constants, 
reducible and irreducible, and alternative sources of information regarding 
path factors. 


5.2.4. Reducible constants 


By reducible configurations we mean configurations which can be con- 
structed by placing together two lower order configurations at one vertex 
or along one side; in a terminology analogous to that of §5.1 we may say 
that reducible configurations are made up of trees with cutting points 
or stars with cutting lines. For example p,, is formed from two triangles 
with a vertex in common, and p,, from two triangles with a side in 
common, 

The method of calculation of such constants is best illustrated by 
examples, and we shall first calculate p,, and p;, in terms of lower order 
constants. We first note that since there are Np, simple I-sided polygons 
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on the lattice, and each polygon has I vertices, lp, such polygons pass 
through any point of the lattice (e.g. four squares pass through each point 
of the simple quadratic lattice). Similarly since there are NV. q/2 bonds in 
the lattice, and each polygon has / sides, 2lp,/q¢ polygons pass through 
any bond of the lattice (e.g. two squares pass through each bond of the 
simple quadratic lattice). 

We now consider two independent sets of triangles passing through 
any point of the lattice, and we form configurations by pairing together one 
triangle of each set. There will thus be a total of (33)? configurations for 
each point of the lattice; among these configurations each p,, will be 
described twice (because of the symmetry with respect to constituent 
triangles); each p;, will be described four times, twice from each of the 
vertices at the end of the diagonal; and each p, will be described three 
times, once from each vertex. Hence we may write 


OD Bash PegtiePa- os Wintoe cle | sutenee (9) 


We next consider two independent sets of triangles passing through every 
bond of the lattice, and pair together triangles as before. There are now 
(6p,/q)? configurations through each link, and hence the total number of 
configurations in the lattice is (Vq/2)(6p5/q)?=18N p,7/¢. We shall form 
Psq twice because of the symmetry with respect to constituent triangles, 
and p, three times once from each side. Hence we have 


US = 20 Os. accuse pe) pion te 10 24 (LO) 
From (9) and (10) we deduce that 
Ope Ss 
(11) 


_ We can proceed in a similar manner to determine p,,, Ys, for the simple 
quadratic and simple cubic lattices. When two squares are placed together 
with a vertex or side in common on these lattices the only configurations 
formed are p74, Pg, OT Py, and we can deduce by analogy with (11) that 

2 
Pia pe = 24, 
: (12) 


32, 
Psr= (s = =) Pa? + 2p4. 


However, for the body-centred cubic lattice other configurations can be 
formed ; if the quadrilaterals have opposite vertices in common pg, will be 
formed, and if they have two adjacent sides in common pg, will be formed. 
Both of these are irreducible constants, and we must evaluate them before 
we can determine p,, and pg, If we put together two independent sets 
of quadrilaterals passing through a lattice point, Pan will be formed twice, 
Prq four times and p, four times as before; but in addition p,, will be 
formed 12 times (6 independent quadrilaterals can be chosen from the 
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pairs of double bonds connecting opposite vertices, and a factor 2 arises 
as in p,, because of the two vertices), and p,, will be formed 18 times 
(6 times from each of the vertices at the end of the diagonal and twice from 
each of the other vertices). Hence the appropriate relation is 


16p4? = 2g, +4P7qgt+ 4pat12pg+18peq- - - + - (13) 


Putting together two sets of quadrilaterals passing through a bond of 
the lattice we will form p,, twice and p, four times as before, and p,, 
twelve times (twice from each of the 6 lines of the configuration). Hence 

322i = 2p t+ 4a t+ l2Pege - = ae 2 eS) 


For the face-centred cubic lattice configurations containing triangles 
enter, and p,,4, P7, must also be included. In the same manner as before 
we can show that to (13) and (14) we must add respectively 


dps + Op, 24peds «+. < = a ee (15) 
2P7¢ + 12Pea- 


Using eqn. (11) and substituting in (7) and (8) we can express p, in terms 
of the q,, 


12p¢=%6— 493 — 39s°(1—3/q).- - - + + - (16) 
The analysis of g, corresponding to (8) is 
7 = 14p, + 28p,, + 8497, + 146, + 28M5q- - - - - (17) 


No irreducible constants arise, and proceeding as before we can derive the 
formula 


14p,=q_+ 7a 2 es 1) (# -3) — 79394(1 — 3/q). Heth ey 


Equations (16) and (17) are valid for all regular lattices but are most useful 
for close-packed lattices. When we come to the general expression for 
Pg ((6) is valid only for the simple quadratic and simple cubic lattices) we 
must introduce irreducible constants as in (13), and we shall therefore now 
consider how these can be evaluated. 


5.2.5. Irreducible constants 


Any lattice constants which cannot be expressed in terms of lower-order 
constants we refer to as irreducible. The most important class of irreducible 
configurations are the simple polygons. After this the most frequently 
occurring class are configurations involving a double bridge, such as p,,, 
P7r7> Poe (Appendix IIT, table A). These can be determined by the method 
of space types, e.g. from fig. 6 we can deduce that p,,= 24 for the simple 
cubic lattice (note that double bridges arising from (c) are merely a 
repetition of those from (a)). But we shall now show that they can be 
calculated from the properties of random walks on the lattice and for three 
dimensional lattices this method is usually to be preferred. 
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Let ¢(x,y,z) denote the generating function for random walks on the 
lattice. Values of ¢ for various simple lattices are given in Appendix II. 
For the face-centred cubic lattice, for example, 


P=(y* +y)(e*+2) + (2+14+2)(e1+4+a)+(et+a)(yt+y), . (19) 


and this can be written symbolically 


peee ri, Maen pvmemos./ lan ha asiy 3(20) 


where the sum is taken over all points which are equivalent under the 
symmetry operations of the lattice (permutations of x, y, z, and inversions 
xa, y>y, zz), From the terms in ¢! we can determine fi, (a, 6, ¢), 
the number of walks which arrive at the lattice point (a, b,c) after 1 steps. 
As before, we can eliminate all walks with immediate reversals; Domb 
and Fisher (1958) derived the generating function for such walks, and 
from their results the following expression analogous to (2) can be deduced 
for Q,(a, b,c) after J steps: 

(8) 


2! 


Q,(a, b,c) =R,—[o(l—1) +. 1]R_.+ (ele) 2 oh, oe 


21 
For the face-centred cubic we find that oo 
fp? =12+4> xy +430? + 23> a2y2 + Sa?y?, ] 
$3 =48 + 45 Say + 24D 2? + 12322? + 18 Y22yz ee 22) 
+ 9> ay + 38> a3y?2 + 6 Sa2y22? + Sady’. 


Applying (21) we derive the corresponding polynomials for no-reversal 
walks: 


y= 4 > xy + 4327 4 2S ax?yz+ Yay?, | 
= 48 + 22S ay + 2430? + 123 2?y? + 18> 2yz ee (23) 
3 Y 
+9 Day + 3>a8y?2 + 6 x7y?2? + diady’. J 


From (23) let us now form the quantity >Q.(Q,—1)/2!, where the 
summation is taken over all points except the origin. Since }ay 
represents 12 equivalent points, Sz? 6 points, }ayz 24 points and Dixy? 
12 points, this sum is equal to 


Wiese footage! 1 241 0,12) 5182... .- . (24) 


This number represents the number of pairs of two-step bridges we can 
select from the origin to the ensemble of accessible points. From the form 
of the sum it is clear that the bridges are all distinct and can only form a 
quadrilateral. On summing over all points of the lattice as origin, and 
dividing by NV we shall form each quadrilateral four times (once from each 
vertex as origin); hence 4p,=132, p,=33. Thus we have confirmed the 
value of p, derived in § 5.2.3. 

If we now form >Q.(@2—1)(Q2—2)/3! we will enumerate all groups of 
three two-step bridges from the origin to the set of accessible points. 


P.M.S. Y 
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Summing as before we will describe each p,, twice; we thus find that for 
the face-centred cubic p,,=36. Similarly, if we form 


X@2(Q2— 1)(Q2—2)(Q2—3)/4! 
we can deduce that for this lattice p,,=9. 

When we consider configurations involving bridges of more than two 
steps (e.g. p7,) we must allow for the possibility of some steps in the bridges 
being identical. Thusif we form }Q3Q, we will not only describe pentagons 
but also the open configurations 


D 8 
eee: 
c 


the latter will be described 4 times, twice from A, and once each from B, C. 


Hence 
2302 = 10p;+ 4 | —<] ] eo epee ae 


and this can be used as a check on the constants which have been determined 
by other means. 
Similarly we may derive the relation 


03Q2(Vs —1)/2!= 2D + 4P6n + 4D 50; oy sy eo Mees 
and from it we deduce that p,,=600 (face-centred cubic). 

No general recipe can be given for calculating all irreducible constants ; 
many of those not amenable to the methods described so far correspond to 
tightly bound configurations and can readily be found by inspection of 
space types on the lattice (e.g. for the face-centred cubic p,,=2). Others 
can be found by more detailed examination of the bridges enumerated by 
(23). Thus to calculate p,, for the face-centred cubic we can examine 
all two-step bridges, and determine how many pairs are nearest neighbours 
after the first step; it will be found that of the four bridges from the origin 
to Sxy there are two such pairs, and of the four bridges from the origin 
to >? there are four such pairs. We therefore deduce that 


2p,, = 1222464690 0 0) ee ee 
and hence that p,, = 48. 


The remaining irreducible constants ps, Pgq given in Appendix ITI, 
table A can be determined from the relations 


DQ203(Q3 — 1)/2! = 2pgq+ 2Dgp + 4Pyy + 4D7_ + 4a + 271, 


+ 8p6.+ 46, + 4D 5a +2/ am | (28) 


2 Q12(Q2—1)/2! = 2pg.+ 2pggt 4p t+ 8P7q+ 12p7+ 6P7, 
+ 48p6q+ 12p6q + 45a: 
The coefficients in these relations become quite complex, and are similar 
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to path factors, but are easier to determine directly because of the absence 
of ‘improper closures’. To evaluate the second line in (28) for the face- 
centred cubic we require the next term of (23) which is 
ify = 264 + 224 xy + 152 Sa? + 136 Dayz + 114 F22y? 
+965 x%y + 36 Sat + 72> x2y22? + 52 Sa8y2z + 24 Fady3 
+ 24> atyz + 16 Saty? + 1228322 + 6 Yaty2e2 
+43 atyse4 Satyt, 9. 2. (29) 
Having shown how to evaluate irreducible constants we can now return 
_ to the extension of (18). For pg we have derived the formula 


16[ pg + 3D¢q — 6g, — 3364] = 4g — 594 — 8093 — 4(295¢3 + G2)(1 — 3/q) 
128 
— 32q493/¢ + aat(— tC my) + 8q37(3¢4 —89q3)/q?. (30) 


We have also obtained an expression for py, but this involves the irreducible 
constant »,,, which can be calculated from the relation 


> Q392(G2— 1)(Q2— 2)/3! = 29, + 2gp + 4070 peo eI) 
We then find 


18 [Pq — 697, + 27, — 1097), + 72P6aP3/q + 8P¢q(26 — 99s) | 
81 
= Jo + 20093 — 9(46%s + Y594)(1 — 3/¢) — 369593/9 + G49 (= + 36— 02) 
— 522q57/q + 274593"/q7 + 274743/97 — 144q1937/¢" 


+q3° (= — = + 15) — 48q,4/q°. Mees be Oyo) 
The determination of path factors for expressions like (30) and (32) is most 
intricate, and hence their verification is of considerable importance; 
methods which have been used for this purpose will be described in detail 
in §5.2.9. 

Relations like (30) and (32) are topological, and can be applied equally 
well to all regular lattices. We have illustrated their use by application 
to the face-centred cubic lattice since this is the lattice for which the 
algebraic technique is the most powerful, and the space type method the 
most difficult. For example we find that p,=301, 376, and this certainly 
corresponds to more than 6000, and probably to more than 10 000, distinct 
space types. 

For loose-packed lattices only even terms occur, and many of the terms 
in expressions like (30) are zero. We can give a general formula for p,, 
as follows: 

20[ P10 — 6P1oe + 2Pee + 2Py— 12M 6a(4—44/9)] 


= 49 — 10q496(1 — 3/4) — 594°(q — 2)(1 — 2/9) 
+ 10q4(1 — ga/9)(4 — 44/9). soe ee 4 (33) 
The values of the lattice constants discussed in the foregoing sections 
corresponding to closed configurations are listed in Appendix IIT, table A 
for a number of simple lattices. 
Y2 
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5.2.6. Open configurations 

In contrast to closed configurations, open configurations have at least one 
vertex to which only one line is attached. Open configurations are all 
reducible since we can always reduce the single vertex; they are very 
numerous and we have therefore not used a symbol for them but represent 
the corresponding lattice constants by a diagram, e.g. 


Darl dl—<i 


The method of calculation of open configurations is closely similar to 
§ 5.2.4, and we shall illustrate it by the example of 


[octal 


We can attach a single bond to every p, in 4(a—1) ways; if the second 
end of the bond remains free we form 


lag; 


but if it is the opposite vertex of the quadrilateral we form p,,, and when 
we sum over all lattice points we will describe each p,;, twice, once from 
each end of the diagonal. Thus 


4(o—1)p4= rete x yer et fh) 


A class of open configurations of particular interest in the theory of 
polymer solutions is that of simple chains with / bonds, the lattice constant 
of which we denote by ¢. A general theorem has been enunciated by 
Sykes which enables c, to be expressed in terms of lattice constants of 
closed configurations, and greatly simplifies the counting problem. A 
detailed proof of this theorem will be available elsewhere (Sykes 1960), and 
we shall therefore only indicate the line of approach. 

Let us add another bond to ¢,; this can be done in 20c, ways, one from 
each end. We form c¢,,, twice and other configurations which are either 
polygons or ‘tadpoles’ (fig. 10a). If we now add another bond to the 
free end of a tadpole, we form either tadpoles or closed configurations 
which are ‘dumb-bells’, ‘orthodox figure eights’, or ‘star figure eights’ 
(fig. 10b, c, d). If we take careful account of the number of times 
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the various configurations are formed, we can derive the following 
relation : 


C44 — 20¢, + 07¢_1= Ip, —(U+ 1)Dy44 


+45 [O-0 | ay [ 00 J +6> [o] > (38) 


dumb-bells orthodox star 
figure figure 
eights eights 
Fig. 10 


(a) A ‘ tadpole ’; Gyan 
(b) A‘ dumb-bell’ ; O—-) 


(c) An ‘ orthodox figure-eight ’; ex) 


(d) A ‘ star figure-eight ’. Qs) 


The sums in the last three terms on the right-hand side must be taken over 
all configurations of the appropriate type formed from (J+ 1) lines. 

It will be seen that using (35) the data given in Appendix ITI, table A 
enable us immediately to write down the first eight chains for the face- 
centred cubic lattice: 


c= 46, } 
C.= 66, 
C3 = 22c, — 121c, — 3p, = 702, 
C4 = 22c, — 121¢, + 3p, — 4p, = 7350, 
Cz = 22c, — 121¢c,+ 4p, — 5p; + 65, = 76266, 
C_ = 22c, — 121¢,+ 5p, — 6pgt 4p¢6.+ 6( Pea t+ Pen) = 786858 
Cy = 22¢, — 121¢; + 6p6— 7p, + 4p74+ 47, + 
+ 6(Drq + Pry + Pry) = 8086074, 
Cg = 22c, — 121c,+ Tp, — 8pg t+ 4(Pgj + Pax) + 4(Pso + Pen) 
+ 6(Pgq+ Pav t+ Pact Pea) = 82848522. 
For other lattices the chain series can be pushed considerably further ; in 
two dimensions special methods are available (§ 5.4.4) and for the simple 


quadratic lattice as many as 18 terms have been evaluated. The numerical 
results derived so far for a number of lattices have been tabulated by 


1 (36) 
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Fisher and Sykes (1959), and are reproduced in Appendix III, table B 
(the difference of a factor 2 arises since Fisher and Sykes consider non-self- 
intersecting walks starting from the origin whereas we consider configura- 
tions.) Numerical values of lattice constants for open configurations of 
up to six lines are also tabulated in this Appendix. 


5.2.7. Separated configurations 

In the Mayer Theory separated configurations are eliminated at a very 
early stage. However, in configuration problems on crystal lattices they 
present a more difficult problem, and can be calculated by a method similar 
to that of §5.2.4 for reducible configurations. Let us, for example, 
calculate the lattice constant corresponding to two separated triangles, 
which we denote by [p3, 3] (the definition being per lattice point as before). 
Consider all configurations which will be formed by pairing together 
triangles from two independent sets of all the Np, triangles of the lattice. 
Each p;,, will be formed twice, when the triangles have a side in common, 
each p,, twice when they have a vertex in common, and each p, once when 
they coincide; the remaining configurations will all belong to [p 5, p,] and 
will be described twice, by symmetry. Hence 


(Mp3)? = 2N Psat 2N Pee + Npg + 2N peed _ . (37) 
or [p3,)3]=$N ps" — Psa — Pee — Ps! 2. 


To determine [p5, 3,93] we must consider all configurations which 
can be formed by one each of three independent sets of all the Vp, triangles 
of the lattice. Proceeding as before we find that 


(1/N)(Np5)? = 6 [Psa P3] + 6[ Pee. P3] + 6[ ps, P3)+ 6[ ps, D3, Ps] 
+ 69. + 8P9q + Pog + bg + bPgn + OPgq + 67. + 6p7, 
+ 6p6_ + 24964 + 6P5q t+ Ds: reer rea eh 4 


The terms in (38), although lengthy, are not difficult to derive; we look 
through table A, Appendix III, and pick out all configurations which can 
be completely described by three triangles. We must now reduce 


[Psa Ps] and [16,3] as in (37): 
NPsaPs= [Psa P31 + Psi t+ Pam t+ Psqt 2Pre+ 3P79 t+ 12p6qt Wsas 
NP ePs= [Poor P3] + 2Poe + 8Pyqt 3Pgq + 2g + Pam + 2Wgq+ Prot Deo: 
(39) 
Hence we obtain 
LPs, Ps, Ps] = tN? ps® —N [peeP3 t+ PsaPs + Ps"/2] 
+ [Doe 29a + 2094 + 2g + Pam + 2gqt Wee + 2g g + Woot Seq 
+ 254+ ps3/3]. tae ih ys een atl) 


This formula can also be checked by the use of finite clusters of points, 
as will be shown in § 5.2.9, 
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Formulae for separated configurations needed for the no-field series 
for the Ising model (s = 4) are given in Appendix III, table C; they suffice 
to determine nine terms for close-packed lattices and ten terms for 
Joose-packed lattices. 

The same technique can be used for evaluating lattice constants for 
separated open configurations. However, they are more numerous, and 
we shall not reproduce any details of the calculations. For the Ising 
model with s=4 open configurations can be eliminated in many problems 
of practical interest. However, for general s and for the Heisenberg model 
these configurations always arise in the counting problems, and we have 
therefore tabulated the terms corresponding to five lines or less in Appendix 
III, table C. We have given the numerical value of the coefficient 
independent of N in an expression like the right-hand side of (40), since 
this is the value required for the appropriate term of the free energy 
expansion (§3.6.1, §5.2.10). (The numerical values are multiplied by 
l! to avoid fractions.) 


5.2.8. Homogeneous non-regular lattices 

We shall now consider the extension of the methods developed above to 
lattices in which points are all equivalent but bonds are not all equivalent. 
We have come across a two-dimensional example of such a lattice in fig. 36, 
$3. The most important lattice of this kind is the hexagonal close-packed ; 
the bonds may be divided into two separate groups of six, those lying in a 
plane of close packing, and those connecting two successive layers of 
close-packing. 

Most of the formulae and methods described previously do not use 
bond-equivalence, and retain their validity. In this category are (2) 
(a detailed justification will be given in §5.2.9), (7), (8), (9), (17), (25), 
(26), (28), (31), (35), (37) and (40). In fact the only feature using bond 
equivalence is the reduction in § 5.2.4 of a configuration with a cutting line ; 
hence an equation like (10) remains valid only if the numbers of triangles 
passing through each type of bond is the same. 

The generation of walks on non-Bravais lattices is not as simple as on 
Bravais lattices, and needs slightly more general treatment. Lattice 
points which constitute the unit cell must be considered separately from 
those within the unit cell, and instead of a simple generating function like 
(19), we are led to the traces of the nth powers of matrices whose order is 
equal the number of points in a unit cell. For the hexagonal close-packed 
the matrix is 2x2, since we must differentiate between odd and even 
numbered close-packed layers. We may write for the position after 
I steps 

y= 1490, T 12 +24)0 4°C,T (2 +24 PU 
10,7 (¢+24802V +... ae eer, (41) 
Here 7 is the generating function in a close-packed plane, (i.e. of the 
triangular lattice), 
Ta=xtut+ytyt+ey +a, ape ae a (42) 
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and z is the direction perpendicular to this plane; U represents a jump 
from a close packed layer to an adjacent close packed layer, and V the 
reverse jump, 

Uaamtyi( tatty), 

V =x 28y-V2(1+2+y), 3 oe ee ee 

UV=3+T. 

Domb and Sykes (1957c) noted that the number of returns to the origin 
in a random walk of 1 steps was exactly the same for cubic and hexagonal 
close-packed lattices (Fisher and Sykes have both given independent 
proofs of this interesting result). Hence from (7) p3, p, and p; are the 
same for both lattices. 

Because of the difficulty of reducing configurations with a cutting line 
the method of § 5.2.5. is the most useful for the hexagonal close-packed 
lattice, and it can also be applied equally well to reducible configurations 
like p;,,, Pep. Owing to the lower symmetry of the hexagonal close-packed 
lattice the classification of equivalent points is more complicated (as is 
also the space-type classification). 

From (41) we readily deduce that 


Y= (eta tyty tay tatty) + ey P(z+z24\(l+at+y), 
Xe=124+4(e+att+tyty*+aey7* +a) 
+ 2(xy + aly) + wy + a-24y + why? + xy-*) 
+ (09 + 0-9 + 2 + yy? + arty 2 + oy?) 
+ x28yll2(z 4 2-1)[4(1 +214 y-1) 
+2WUe+ytur+y +a + xy) 
+ (xy + arty + at) 
+ (27 +27)[34+ (eta *+yt+y1+ay+ay)]. 
ae (44) 
Groups of equivalent points have been bracketed together; although the 
expressions are long they have simple geometrical interpretations, as will 
be seen by drawing them on the triangular lattice. Forming 
%9192(M2— 1)/2! 
we find that p,;,=36, as for the face-centred cubic lattice. However, on 
forming }Q.(Q@2—1)(Q2—2)/3! we obtain 
2Neq = +.6+ 4.64 4.64+2=74. «| Ait eats eh OD 
Hence p,,,=37 and differs from the value of 36 for the face-centred cubic; 
from $Q.(Q2—1)(@2—2)(Q.—3)/4! we find that p..=9 as for the face- 
centred cubic. 

We may proceed further as for the face-centred cubic lattice. In fact 
the general result can be established that the number of triangles, quadri- 
laterals, and pentagons passing through each type of bond is the same; 
hence the reductions of § 5.2.4 remain valid for the hexagonal close-packed 
lattice. We then find from (30) (p,,=2, as for the face-centred cubic 
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lattice) that pg=42066. Similarly p,,=51 for the hexagonal close-packed 
lattice, and hence p,=301 358. 

There are more ‘tightly bound’ configurations on the hexagonal close- 
packed than on the face-centred cubic, and in this sense it may perhaps 
be regarded as slightly more ‘close-packed’. 


5.2.9. The use of finite clusters of points 


The configurational formulae derived so far in this chapter have in no 
way depended on lattice structure, and apply equally to finite clusters of 
points with connecting bonds. For such clusters the enumeration of 
configurational constants is often quite elementary, and hence they can 
be used as a check on topological formulae like (11) (17) or (40). We can 
readily extend our lattice classification to finite clusters. A cluster will 


Fig. 11 
ee |X 
(2) 
le OS 
(6) 

Hy YY 
(1) (c) (2) 


Some finite clusters of points. (a) Close-packed clusters; (b) Loose-packed 
clusters; (c) Other homogeneous sets: (1) Cube, homogeneous and 
regular; (2) Prism, homogeneous and non-regular. 


be described as regular if all points and bonds are equivalent, and 
homogeneous if all points (but not necessarily bonds) are equivalent. For 
finite clusters we shall use the term configuration constant instead of lattice 
constant, which we shall define per site as before. Thus, if the cluster 
consists of NV points, it will contain Np, triangles, Np, quadrilaterals, etc. 

We shall first make use of clusters with a high degree of symmetry. 
A close-packed cluster of N points is a set of N points every pair of which is 
connected by a bond. For four points this corresponds to a tetrahedron, 
as shown in fig. 1la, and for six points to a hexagon with all diagonals 
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connected. A close-packed cluster of V points has coordination number 
(N—1) which we denote by (¢+1) in accord with our practice for 
standard lattices; we shall often refer to this cluster briefly as CP(J). 

We can also consider random walks on clusters, and can derive a formula 
for r, the number of returns to the starting point after / steps. To do this 
we do not make use of generating functions, but instead establish a set of 
difference equations which can usually be solved quite easily. For a 
close-packed cluster the procedure is particularly simple, since all (NV —1) 
points other than the starting point 0 are equivalent. Denoting one of 
the equivalent points by P the difference equations are 


r141(0) =(o+ 1)r(P), \ Weer 10 
1144(P) =7,(0) + or(P). 


From this it is easy to show that r,(P) is the coefficient of z’— in the expan- 
sion of [l1—ov—(o+1)a?]- and hence that 


n= TS [e+ + (-1)1. ese rae a 


oO 


Domb and Sykes (1957 a) observed that the formulae (2) for q, remained 
valid for homogeneous clusters of points. The argument given by Domb 
and Fisher (1958), which is based on the use of generating functions, is no 
longer valid in this case, but the following alternative derivation, due to 
Fisher, provides a clear insight into the essential topological nature of the 
formula. We can rearrange (2) so as to give 7, in terms of q, as follows: 


U(1—1) 


71 =G + olq,_»+ 0° — 


Q-at...+B(c) . . . (48) 


where the q, series terminates at qg, if / is odd, and at g,if/is even. Each 
term then has a simple topological interpretation, the totality of random 
walk returns being sub-divided into classes each of which contains exactly 
r immediate reversals. The term in qg, corresponds to a walk with no 
immediate reversals. The term in q,. corresponds to a walk with exactly 
one immediate reversal; this may be specified by the first step of the 
reversal which can occur at any one of the / steps of the walk (the occurrence 
at the /th step is interpreted as an improper closure); also there are o 
possible directions for this reversal, and when the reversal has been 
removed we must be left with a non-reversal walk of (J—2) steps. The 
term in g,_4 corresponds to a walk with exactly two immediate reversals 
which can occur at any pair of the / steps of the walk, and gives rise to the 
factor ‘C,o* (typical contributions of various terms are illustrated dia- 
gramatically in fig. 12). Similarly the term in q,_», corresponds to exactly 
r immediate reversals which give rise to the factor ‘O,o". The final term 
Bc) arises only if J is even and corresponds to walks with no closed 
configurations, which we may describe as walks on a Bethe lattice (§ 4.6) ; 
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it can be calculated by putting 95, q4,... successively equal to zero in (2) 
and deriving the corresponding r,, PRG 


If we now apply (2) to (47) we can readily determine the values of q, 
as polynomials in o: 


1 =9, 

J2= 9, 

q3=0(o+1), 

qa=0(o+1)(c—1), ieee (49) 
95=9(o+1)(o?— 80+ 2), 

dg = 9(o + 1)(o? — 307 + 5a — 2), 

97=0(0 + 1)(o4 — 30° + 70? — 80+ 8). 


Consider now a configuration with v vertices on a close-packed cluster. 
The v vertices can be chosen in YC, ways; if we use the symbol 7,,” to 
denote the corresponding configurational constant (defined per site as 
before) we may write 


Pix” =4z"(o + 1)o(o—1)...(o-—v +3), 0s ee OU) 
Fig. 12 
(2) (0) (c) 
Analysis of random walk returns. Typical contributions to rz. (a) One 


immediate reversal; (b) Two immediate reversals; (c) Involving no 
closed configurations. 


where the a,,” are pure numbers independent of o which can easily be 
determined by the symmetry of the configuration. Thus a,6=1/12 since a 
polygon can be described in either direction from any ofits vertices; similarly 
Mee? = 1/8, A5q*= 1/4, a32=1/6. Hence we can calculate configurational 
constants for close-packed clusters of various sizes, and use them to test 
topological formulae derived previously. 

In particular, if we substitute into (5) from (49) and (50) for close-packed 
clusters with general oc, and compare the polynomials in o, we can derive 
relations which must be satisfied by the path factors A,,,(1). Itis convenient 


for this purpose to write (49) in the form 


Q=u%(o+1)o+B(o+1)o(o—1)+y(o+1)o(o-1)(o—2) +... mk 


320 C. Domb on 


where the constants «,, f,,y,...are easily evaluated numerically. Taking 
account of the vertex classification we may re-write (5) in the form 


g= > Agra?) aoe ee 


238<l v=3 


using (50) and (51) and comparing both sides, we find that 


eS Nox? Vasa? ae 1 


%;8<l 

>» s08(l)dsa* = By e ° . . ol oe (53) 
@58<l 

2s Asx? (gn? =p 
L538 


the sums being taken over all configurations with 3, 4, 5... vertices respec- 
tively. We thus obtain (/— 2) relations to be satisfied by the path factors 
of q. For 1=6 these relations suffice to determine the path factors 
uniquely; for />6 there are not sufficient relations to determine the path 
factors but we obtain (J—2) checks on the correctness of the values 
determined by the direct method. 

For checking topological formulae on loose-packed lattices it is convenient 
to introduce a loose-packed cluster of N points (NV even) which consists of 
two equivalent sub-clusters of V/2 points, no points in the same sub-cluster 
being connected, but each point of a sub-cluster being connected to all 
points of the other sub-cluster (examples are given in fig. 116). It can 
readily be shown that for a loose-packed cluster LP(V), 


du=07[o'2+(-1)P/(c+1), . . .-. . . (84) 


and we can proceed in a manner closely similar to that for close-packed 
clusters. 

Any homogeneous set of points can be used similarly. For the vertices 
and edges of the cube, which form a regular set, it is easy to derive the 
formula 


r,=3(14+-3")/4 (Peven) . 2 =) aleeeoml 


and hence q, can be found from (2). The triangular prism (fig. 11c) is 
an example of a homogeneous non-regular set. 

It is interesting to observe that the determination of path factors is 
related to certain classical topological problems. Thus Eulerian paths 
which arose out of the problem of the Kénigsberg bridges, are paths on a 
graph which describe each line of the graph once only (Berge 1958, p. 106). 
The problem which we have considered is more complicated since we wish 
to know the number of independent closed paths of J steps, each line of the 
graph being described at least once. Similarly the problem of Hamiltonian 
paths and circuits, the specific requirement of which is that no point may 
be visited more than once (Berge 1958, p. 162), is related to the deter- 
mination of p, and ¢, for the appropriate configurations. Sykes has used 
formula (35), and the specific two-dimensional methods to be described in 
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§ 5.4.4 to enumerate all polygons and chains on a regular dodecahedron. 
Some features of this enumeration correspond to the classical solution of 
the problem by Hamilton. 

When we come to consider non-homogeneous clusters it is useful to 
define 7;, q, as averages, taking all vertices of the cluster in turn as origin. 
The determination of 7, is not too difficult, but no simple formula such as 
(2) exists for gq. The evaluation of q for such graphs would represent 
considerable progress in the direct determination of path factors, and an 
approach by the authors using symmetry properties to reduce the 
appropriate difference equations has led to useful preliminary results. 

Several of the formulae derived in the previous sections of this chapter 
do not require the cluster to be homogeneous, and can be applied quite 
generally. or example the formula (40) for [p3, p3, 3] can be applied 
to the cluster represented by the configuration p,,; for this cluster 
[Ps Ps, P3]=9, P3=3/7, Deo= 2/7, Yo,= 1/7, and all other terms in (40) are zero. 
It will be verified that the equation is satisfied. 


5.2.10. Series expansions for the simple Ising model in terms of lattice 
constants 
The lattice constants listed in Appendix ITI arise in all high-temperature 
expansions of interacting systems on crystal lattices. Usually a further 
calculation of mean values of interactions for these configurations must 
be undertaken before the series can be specified, and this calculation will 
be considered in §5.3. However, for the Ising model with spin $ the mean 
value problem is particularly simple, and has already been discussed in 
3.4.1. 
In the absence of a magnetic field the partition function is given by 
eqn. (34) of §3.4.1, and we can write down the first few terms of p(r) as 
follows : 


+ Pent 3Peqt Ps); 
p(9) =4N2p33 — N*(DsaP3 + Ps" — PePs — PsPa) + N(— 29. 2Poz— 4P oo 
— 29n+Po—Psa— Psv— 2Pse— 2Pap — 2Psq— Pro + Pre 
= 29 + 27g — 2D7n— Pov + 8Poat 2Psa+ BPs): 
; (56) 
For loose-packed lattices another term can be added: 


(10) = N26 94—N (2Pi0n + 2P 100 — P10 + Pow + Dom + 2D ge 
+ 2Ds + 2D1q)- 5 a e (57) 
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To derive the partition function per spin we know from §$3.6.1, 3.6.2 
that we must put V =1 in (56) and (57); for the logarithm of the partition 
function we must take the coefficient of N. It will be noted that the 
coefficients of the latter series are all linear in configurational constants, 
and no configurations with cutting points enter; this result exactly 
parallels the Mayer theory. 

The disappearance of product terms on taking the logarithm is an alge- 
braic property, and does not depend on N being large ; it will therefore also 
be valid for finite clusters of spinst. We can write quite generally 


InZ,/=NIn2+W In(coshK)+N > aw’ (w=tanhK) (58) 
1=1 

where the a, are the linear combinations of lattice constants given by the 
coefficients of NV in (56) and (57). For a finite cluster the series for Z,/ 
must terminate after a finite number of terms, and this means that the 
corresponding combination of constants must be zero for the cluster. 
In fact we can use the partition functions of finite clusters to obtain some 
information about the series Saw’. Because of the linearity of a, in the 
lattice constants this series can be put in the alternative form. 


DPeafea) - - + - + + + + (69) 


where f,,.(w) 1s a power series in w. Substituting the partition function 
for a triangle the only lattice eit which enters is p,; hence 
f,(w)=m(1+w).. . es GAN 
Similarly for the cluster fase aa to p;,, the only terms which enter 
are pz, p, and p,;,; hence we deduce that 
In (1 + 203 + w4) = 2f3(w) + fa(w) +foa(), 
Fre(to) sein (1+ 208 + 4) ee ye 
Se (1+ w)?(1 + wt) 
This method of deriving the functions ier (w) is reminiscent of that used by 
Yvon (1945), Rushbrooke and Scoins (1955), and referred to in §§ 4.5.1, 
4.5.2. The present approach can also deal with non-homogeneous con- 
figurations. It can be used to provide high-order approximations in 
closed form to the no-field partition function; and with its help separated 
configurations and configurations with cutting points need not be considered 
for the no-field series. 
The series for the magnetic susceptibility in zero field is obtained from 
eqn. (155) §3.6.2, using the relation 
g2 
Xo= — ON /OH* = kT —, (In Z?) Slatin ek sem (OS) 
(/=Free energy) (H = Magnetic field) 
Remembering that 
7=(1—p)/(l+yp)=[1—- eee 2 mH /kT)]/[1+exp (—2 #H/kT)] 
we must pick out the terms in 7? in the above expansion, and Heal arise 


ee er 
+ This result can best be established rigorously by the use of pseudo-lattices 
as in § 4.5.1. 
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from configurations with two odd vertices. Some examples of these 
configurations, which we shall refer to as magnetic configurations, are given 
in fig. 13. It is then easy to derive the formula 


we co 
X\Oues Tail 1 oe 2 >, byw! . . . . ° ° (63) 
where 6, is the coefficient of V in the sum of magnetic configuration constants 
containing / lines. 

For the susceptibility coefficients b, Sykes derived a theorem analogous 
to (35) for chains which enables them to be expressed in terms of closed 
configurations. The relation is 

6,1 — 20b, + o7b,_, =Sum of constants of 
closed configurations of 
order (J+ 1) or less. : th on (64) 


= 
ae /\ 


(separated ) 


Some examples of magnetic configurations with six lines, 


The configurations which contribute to b, are considerably more complex 
than the chains c,; hence the proof of this theorem is also somewhat 
complicated and is given elsewhere (Sykes 1960). 

With the aid of (64) we can write (63) in the form 


UZ 


2 oo 
Mim ee (1—ouy-4| 1—(o— 1)— out + > an vars mt(OD) 
1=3 


where the d, have the following values: 
ds = — 6p3, 
dy = — 8p,, 
d;=6p3— 10p5+ 85a 
dg = 6pg + 8p4+ 28p5q— 126+ 8(Pea t+ Pov + Poe)» 
dy = 10p5 + 36p5q+ 16 p54 + 30D, + 16p6,— 14p, ; (66) 
+8(Drq+ Prot Pret Prat Pret Prog) + 4470 

dg = — 63+ 84 — 28P5q + 12 p64 56D 6q + 24265 + 8P60— 48 Pea 

+ 32(DyqgtPrp + Dre) + 16(Prat Pre t+ Prs) + 8P79— 16pg 

+ 8(Deqg+ Psp t+Pset+ Psat Pset Part Pao t Pant Psj + Pax 

+ Dg + Pap + SPs + 6ss) 
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For loose-packed lattices we can add another term : 
dy = 8P7q+ 16(Pget+ Pan + 2Par) + 8( Po; + Pox + Pa + Pom)- (67) 

It will again be noted that the coefficients d, are linear combinations of 
the lattice constants p,,. The formulae can be applied as before to finite 
homogeneous clusters, which can be used to determine the coefficients of 
particular configuration constants in d,; the application to non-homo- 
geneous clusters needs more careful consideration. 

No general rule has been given for determining the weight factors of 
lattice constants in (64) or (66) analogous to those for chains in (35), and 
in fact these factors must be determined specifically for each configuration. 
However, Sykes subsequently showed that a further transformation was 
possible in which the number of constants occurring on the right-hand 
side of (66) was substantially reduced, and a clear definition of weight 
factors was possible. If H(w) represents the internal energy, then the 


series Siu! in (65) may be put in the form 
3 
2wH(w) 
J 


DY dw'= + 
3 


where the g, are linear sums of lattice constants of order / only which are 
zero for 1<5. The detailed proof of (68) and specification of g, will be 
given in the paper by Sykes (1960) referred to above. 

Formula (68) has enabled the suceptibility series to be substantially 
extended, particularly for two-dimensional lattices for which E(w) is 
known exactly, and 15 terms have been evaluated for the quadratic lattice, 
12 for the triangular lattice, and 24 for the honeycomb lattice. Several 
useful physical conclusions on the form of y, near the Curie point have 
resulted from the use of (65) and (68). In the case of an antiferromagnetic 
(for which w must be replaced by —w), it has focused attention on the 
analogy in critical behaviour between y)(w) and H(w) (Fisher and Sykes 
1958). 


+8(1+0)? S gw,  . . . (68) 
In 


5.3. Mean Values of Configurations 
The partition function for the Ising model with general s can be written 


8 
Zy'= > [I exp(4Ks,,s,;)[] exp (2Ls,,) (L=«H/2skT). (69) 
8yj=—8 6,9 i 
Here we have taken K=J/4s?kT' in accordance with §2.2, so that the 
maximum and minimum interaction energies remain +J; the maximum 
magnetic moment remains », and we have maintained agreement with 
our previous results for s=}. 

For the Ising model with s=} the relation § 3.4.1 (32) means that each 
bond occurs only once in an expansion of the partition function such as 
§ 3.4.1 (33). For s=1 a somewhat analogous treatment is possible, but 
for general s the exponential must be expanded as follows: 

7 \2 c\r 
(4K) , GK) 


exp(4Ks_;s.;) =] + 4K (s.;8,;) + D1 (8.8,;)° bee ai ery ae (S-iSz3)"+ (70) 
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When we substitute the expansion (70) into (69), the term in (4K) will 
still consist of contributions from lines, but multiple occurrences of bonds 
must be allowed. (Typical configurations are shown in fig. 14.) Also 
each configuration will be multiplied by an appropriate product of Sai 
and s,; there being one s,,8,, for each bond; this product must be summed 
over values of s,, equal to —s, —(s—1)...(s—1), s. We thus see that each 
term in the coefficient of (4K)! consists of two factors, one representing 
the number of occurrences of a configuration on the lattice, which we — 
denote by C, and the other which we may call the mean value for the 
configuration, and which we denote by {C;} (this will be defined precisely 
later for particular models). If for the configuration C, the multiplicities 
of the bonds are a, £, y... the coefficient of (4K) in (70) is 

OC OR Visage Ne eet ee dare tem it) 


Fig. 14 


Ly.) 


Typical configurations with mutiple bonds. 


the sum being taken over all configurations of / lines. 

We may use the tables of lattice constants in Appendix III to determine 
C, but we must include a factor to take account of possible change of 
symmetry; this factor is easily found by inspection, and has the values 
3, 4, 6 for the configurations shown in fig. 14 (the appropriate values of 


C, being thus 33, 454, 6 [ a | ) : 


For the Heisenberg model s,,s,,; must be replaced by 
S8j( = Sai5nj + SyiSyq + SaiS 29) 
where s denotes the vector spin operator. The components of s; do not 
commute with one another, and the form (69) is no longer correct and 
must be rewritten 


Zy* =(exp| 4KYss,+2033q |), Been Se cigh (42) 
459 a 


where the brackets { ) denote the trace. The second term represents 
the z component of total spin and commutes with the first term. Hence 


(72) may be rewritten 
Z,% = (exp (1K ¥ss;) [[ exp (2L8,;)) ee eke 2.50.7) 
ij i 


since the various s,, commute among themselves. The exponential can 
be expanded in powers of 4K, as before, and contributions of various 
configurations sorted out. Formula (71) remains valid but now in 


P.M.S. Z 
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calculating the mean value {C;} care must be taken to consider separately 
all possible permutations of the factors s;s; which make up the configura- 
tion. The contributions from these permutations will differ because of 
the non-commutation of the factors, and it is this feature which substan- 
tially increases the difficulty of calculations for the Heisenberg model. 


5.3.1. Ising model general spin s 
Let us define 
& 
tb= . > -exp(2Ls,)=p +e s+... eee. 
re . (74) 
p=exp (—2L)=exp(— #H/skT). 
Then in any configuration each free vertex gives rise to a factor ft), and 
the first term in the expansion of (69), corresponding to a configuration 
with no lines, is 4%. Similarly a vertex which is the meeting point of a 
lines gives rise to a factor 
t= St-§ + (8 —1)4u-“ YD + 2. + (—8 +1) 48-9 + (—8)us. . (75) 
We shall take out the factor ¢,¥ and write 7, =t,,/t). The mean value for a 
configuration with vertices of multiplicity a, b, c...is then defined by 


{Oy} = eres ae » fat ics Seen see ee 
and 
Ligh HTP | ae pS (4K) > Cyrat pT - .-/alBly!.. ] $< toca aan 
i=1 configurations 


with li 

In the absence of a magnetic field ade zero for odd a, and only con- 
figurations all of whose vertices are even contribute. For the suscep- 
tibility in zero field configurations with not more than two odd vertices 
contribute. To evaluate this susceptibility we must expand each 7, as far 
as the second term in H, and pick out the coefficient of H? in (77). The 
coefficients in these expansions are given by successive derivatives of ty 
[=(u-§—pSt)/(1—y)] at w=1, and can readily be obtained from the 
Taylor expansion of t) about ~=1. We find that 


hi 3 E+ 2 
xX 8X —6 
n= t+ 15 Zs al X=s(s+1), 
OX (are 
Teric te a5 “3 nLite, 
48X2— 76X +30 7, 
1 a 5 sf ae | L= mH |2skT, (78) 


3X2_3) 
oe Se ee eC 


3X2—3X +1 40X83 — 110X2 iA — 42 
T= a apatiant + a ; 
5 
2X [35X38 — 70X24 63X — 21 
= pbieertote Acer Eh |. 


T= ea 


3 105 


— 
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Using (78) and the configurational data of Appendix III series expan- 
sions for the partition function and susceptibility in zero field can be 
written down fairly readily. The expressions are somewhat complicated 
for general s, and are therefore not reproduced here. Some results for the 
face-centred cubic lattice are given in the paper by Domb and Sykes 


(1957b), where the physical conclusions resulting from change of s are 
also discussed. 


5.3.2. Heisenberg model spin s=4 


The major additional difficulty associated with the Heisenberg model is 
the determination of mean values of non-commutating spin operators for 
the various configurations. Although in principle the basic problem is 
the same for all s, we shall consider first the case s=4, since simplifying 
features then exist. The mean values for the first four terms in the 
expansion were first derived by Opechowski (1937), and a serious error 
was corrected by Zehler (1950). Subsequent additional calculations 
were performed independently by Rushbrooke and Wood (1955) and 
Domb and Sykes (1956). 

By analogy with (76) we define the mean value {C;} of a configuration 
by the relation 

(Ay) 
to {Cy} = > ((6;9;)"(6;6;.)°(6,.0,)" es Il PSY (wHo,/kT)) 


permutations 


(79) 
ty = (ut + pi. 
Here the o; are the Pauli spin matrices, and the sum is taken over the / 
lines and v vertices of the configuration ; 
(Ay) 


permutations 
denotes the average over all the possible different orders of the factors 
(o,0;), (6;6,),.-.. Equation (77) then remains valid with 7,7,7,... 
replaced by {C;} as given by (79). 

The direct approach to the evaluation of (79) consists in enumerating 
all permutations for a given configuration, and reducing each permutation 
using the commutation relations for the o;._ It is convenient to use suffix 
notation r,s,¢ for the coordinates x,y,z and to employ the summation 
convention; we may thus write 


GC Ore tC Gysks «ee - .- (80) 


Each product can then be reduced to terms which are linear in o,, (for 


each j). For example, for the configuration A we must reduce 
2 3 
(Sp + b€puF ur) (Ops + 1€ psy v2) (Os4 + 4 €stw% ws) = 8 — 25 1% v2 — 20 y1F v3 + 20 y29 v3- 
e e e e (81) 
For this product we do not need to consider permutations since the trace 
of a product of operators is unchanged by cyclic permutations, and points 


Z2 
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2 and 3 can be interchanged by symmetry; in the more general case the 
number of permutations is appreciably reduced by these properties. We 
must now multiply (81) by exp [ #H(o., + 0,9+0,3)/kT'| and take the trace 
of the resulting product. It is easy to see that only o, terms will give any 
contribution; for any particular vertex the absence of a o, gives a factor 
(u-"2+4 2) and the presence of a o, gives a factor (u~?—j4?). Hence 
for this configuration 
{O}=3=27? [r=(1—p)/(1+p)), = 2 2 esas 

Rushbrooke and Wood (1955) showed that a number of elementary 
relationships could be used to reduce configurations before applying (80). 
Thus if we denote 0,0, by (12), we have 

(12)2=3 — 2(12), 

23)[(12)4 (13)]=(12)+ (13), + 2) ee 
)(23) + (23)(12) = 2(13). 
Nevertheless the labour involved in deriving higher terms is large and it is 
desirable to devise independent means of checking the mean values. 

Domb and Sykes (1956) suggested an alternative approach which could 
be used to derive mean values, and could serve asacheck. This alternative 


method calculates the averages of permutation operators P;;, rather than 
spin operator o,o;, using the Dirac relation 


Py =4(1 +6,;). * * . S ° . 4 (84) 


( 
(12 


It is much easier to determine the resultant permutation which is the 
product of a number of permutation operators; if this resultant permuta- 
tion is expressed in the form of cycles, each of these cycles can be treated 
independently, and the averages readily determined. Thus for the 


1 
configuration /\ considered above, 
2.8 


I eS 
PuPuPa=(; : >) = (123), ~\ i Oe 


Taking o,;, 0,2, 0,3 a8 diagonal, the only diagonal terms in a cyclic permuta- 
tion correspond to states which remain unchanged by the operation of the 
permutation, i.e. those for which all elements of the cycle are the same, 
+ or —. Hence the only diagonal terms in (85) are represented by 
(+, ++)(4+, —-)(-, ++)(-, ——) and 


(PypPogP 3) exp [H (0,1 + 029+ O,3)/KT"]) = 9? + U2 + lin Uae 
(86) 


To determine the mean value of the P;; interaction as defined by (79) we 
must divide (86) by (u-!? + ."2)3, and we obtain the value 


(Vapel) (2. S02. cian 0 Tae eg 
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The resulting mean values are tabulated by Sykes (1956) for all con- 
figurations of up to six lines, and enable six terms of the complete partition 
function in a magnetic field to be written down. The zero field partition 
function and zero field susceptibility agree with the values derived 
independently by Rushbrooke and Wood and can be written as followsy : 


In Zy% = N In2-+ Se,K/21!, ; 
12 
of 8(o + 1)/2, 
ég= —3(o+1)4+ 189, (38) 
é,= —3(o + 1)(40 + 1) + 72p4— 72s, 
é, = 60(o + 1)(20 + 1) —360p,(20 + 1) — 480p, + 360p,, 
€= 12(o + 1)(440?— 100+ 13) + 216p,(440 + 3) — 288p,(200 + 3) 
—~ 3600p, — 59045, + 2160p, — 2880 pg, + 2887%¢,, 


J 


kT x eT 
=1 K'/24, 7 
Bae ta 24 / 
j= (e+), 


ja=2(o+1)(o—1), 
J3 = (e+ 1)(607 — 120 — 2) — 60p,, 
 jg=(o+ 1)(2408— 7202 + 480 + 40) — 240p,(20 — 1) — 320p,, 

J5= (+ 1)(12004 — 48003 + 6400? + 1600 — 104) 
— 12p,(3000? — 6400 — 140) | (89) 
— 16p, (2000 — 160) — 2000p, + 1440p,., 

jg =(o + 1)(7200° — 36000 + 696003 — 4064.0? — 51680 — 1232) 
— 12p,(240003 — 82800? + 35040 — 816) — 16p,(18000? 
— 4992c) — 20p;(12000 — 1128) + 8p,,,(2160c + 3888) 
— 144005 + 9792p6q 
+ 921676, + 295685, + 307264. | 


As an additional check on the mean values Domb and Sykes evaluated the 
exact partition functions for finite clusters such as 


Pas Pear Co 


and expanded them in the same manner as (77) (using the mean values 
defined by (79)); if different interactions are used for different bonds 


K! a" 
(e.g. TAN ) 
K 


the mean values for configurations with repeated lines can be derived 


pm e 
+ The subsequent calculations of Rushbrooke and Wood 1958 for general s 
enable the term e¢, to be added. 
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directly. For homogeneous clusters direct verification by substitution 
is possible in formulae such as (88) and (89). 

The physical conclusions to be drawn from the application of formulae 
(88) and (89) have been discussed in the papers by Domb and Sykes and 
Rushbrooke and Wood. For the Heisenberg model it seems clear that 
additional terms in these formulae could be very useful; however, the 
methods described above do not seem capable of extension without very 
great labour. Domb and Sykes (to be published) have therefore recently 
investigated an approach rather like that of (59) in which the free energy 
is expressed as a sum of contributions from finite clusters. Repeated 
factors and permutations would then be avoided, but open configuration 
cannot be eliminated for the Heisenberg model, and must be taken into 
account. 

We shall illustrate the method by writing down the first three terms of 
the expansion. The finite configurations required are 


IN eae 


the lattice constants of which we denote by p,%, po, ps, ps, pz; we 
denote the corresponding free energies by f,%, f,®, f;%, f,°, fz. The free 
energy for a general lattice can then be written in the form 

P= N[pyPb,® + pops + ps hs + pyh3 + ahs + --- ]- (90) 
But this formula applies equally to any cluster. Hence applying it 
successively to the above configurations we deduce that 


f= $y, 
fr = 24, + $0, 
faP = 34, + 34,0 + 64, - >. en a ee 


| fx = 34, + 2b. + $9 
eS 36,0 ot 36.) oe 3; 


from which the $’s can readily be obtained in terms of the f’s. High- 
temperature expansions for finite clusters can be obtained from the traces 
of powers of matrices, and it is possible that computers could hence be put 
to effective use for this problem. 


5.3.3. Heisenberg model general spin s 


The extension of high-temperature series expansions for the Heisenberg 
model to general s was carried out by Brown and Luttinger (1955) who 
derived four terms of the susceptibility expansion. A comprehensive 
attack on the problem was subsequently provided by Rushbrooke and 
Wood (1957, 1958) who obtained two additional terms of this expansion 


(incidentally correcting a small error in a calculation of the fifth term by 
Brown (1956)). 
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The mean values for any configuration for general s depend on the 
multiplicities of the vertices, as for the Ising model (§5.3.1.), and will be 
functions of the 7’s; but the traces of products of spin operators must be 
averaged over all aoeatls permutations as in the previous section. The 
properties of spin operators for general s can be handled in the standard 
manner by taking s, as diagonal and introducing s+=s,+is,, s-=s, —8y 
(Mandl 1954). In re of these operators, the Hoimenbers interaction 
can be written 

SS 558 gt $(8;"S;. 18; 8," )e. pee gen els 25 (92) 


Operators corresponding to different vertices commute, and in any product 
of operators corresponding to a single vertex, the only diagonal terms are 
those with equal numbers of s+ and s~ operators. Also 


sts—= s(s “ 1) ai hs, = as sts, ae 8st woe ist, 
s-st=s(s+1)—fhs,—s,?,  s-s,—s,s-=hs-, 


and with the aid of these relations the diagonal terms in a product of 
operators can be reduced as polynomials in s, When we multiply by 
exp (2Ls,,) and take the trace by summing over the (2s +1) possible values 
of s,,, each s,,“ gives rise to a term 7,, as for the Ising model. In this manner 
we could ane the complete partition function in the presence of a magnetic 
field. However, we are mostly interested in expansions for the partition 
function and susceptibility in zero field, and Rushbrooke and Wood 
(1958) showed that if the terms in these expansions were considered 
directly, a considerable simplification arose because of the symmetry of 
the spin operators in 2, y, z. 

For the no-field series by analogy with (79) the required mean values 


are given by 
(Av 
(2s+1)%O},= >  C(55))*(5)8;,)°(8,8)”...). (98) 


permutations 


Rushbrooke and Wood (1958) give an example of a particular calculation 


for the configuration 
3 
yay 


There are twelve permutations of the product (93) : 


( 
(12) (12) (18) (23) Ava eet ) 
(12) (12) (23) (13) (12) (23) (13) (12) (23) 
(12) (13) (12) (23) (18) (12) (12) (23) (18) 
(12) (23) (12) (13) (23) (12) (12) (18) (28) 
Because of the invariance of the trace under cyclic permutations only three 


of these products need be considered, (12) (12) (13) (23), (12) (12) (28) (13) 
and (12) (13) (12) (23), and of these the first two must be equal by symmetry 


(13 
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inland2. Using suffixes r, s, t, u...for x, y, 2 as before, and summation 
convention, the right-hand side of (93) may be written 


$1281 ,SopS (S9/81wS3uS2v83y + 8178 o754/53/51 1820820830 | 


= 4[2 (81,5181) (SorS2p820)(SsuSar) + (S151 S14) (SopS2482u) (S3x53 v) |. on 


Now it is easy to show that the only products of two spin operators giving 
non zero trace are 


(8,2) = (8,2) = (8,7) =4(s*?)=4F(28+1)X [X=s(s+1)] . . (95) 
and the only corresponding products of three spin operators are 
(878482) = — (8788y)= F(2s+1)X. . . . ~ (96) 
Hence we find that (94) reduces to 
$16(85y8,)"<82" )]; 
and that for this configuration, from (93) 
1 
{Cho= soy 54 X3, . . . . - . . (97) 
Fig. 15 
360 360 
60 15 


Number of independent permutations to be considered in various configurations 
for the Heisenberg model (after Rushbrooke and Wood). 


All configurations must be treated in this manner, and for the higher 
terms the work is very lengthy. Figure 15 represents a number of con- 
figurations, and the corresponding number of permutations to be taken 
into account. Rushbrooke and Wood list all non-zero traces, such as 
(95) and (96), of products of up to seven spin operators. They also enunciate 
a theorem which sometimes reduces the work entailed by the loosely bound 
configurations, such as simple polygons. 


Theorem I. If a connected configuration C, has somewhere three 
consecutive single lines (7j) (jk) (km), the points j and k being not otherwise 
involved in the configuration, and if O,_, is the configuration obtained by 
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removing the centre one, (jk), of these lines, and making j and k coincide 
(the other points of the diagram remaining unchanged) then 


AC ol SA XI ye ee (98) 
where 64;,,=1 if the points 7 and m coincide and 0 otherwise. 

Rushbrooke and Wood (1958) then proceed to show that even for the 
zero-field susceptibility expansion only traces similar to those in (93) need 
be evaluated. To see this we must replace a; by s, in (79) and expand {C}} 
as far as the second order in H. The first order terms i in H are zero since 
there is no spontaneous magnetization; and the second order terms 
involve averages like 

(Av) 
{(s,8;)*(S;S;,)4(s 181)” Sen Sig) E ° a O (99) 
permutations 

where p and q are any two points of the lattice. Again all possible per- 
mutations must be considered, but only of the s,s; in C, and not of s,,,,8,,. 
The configurations involved in (99) may be Feprosented by a set of / lines 
and two crosses (corresponding to » and q). The two crosses may lie on 
the configuration C; or be partially or wholly separated from it; they may 
also coincide in which case they will be referred to as a double cross. 
Rushbrooke and Wood then give three theorems which show that double 
crosses and separated crosses can be ignored; and that averages for a 
configuration C, with two crosses can be replaced by averages for a new 
configuration C),, without crosses. We shall use the symbol C;* to 
represent a configuration with a single cross, C;,* a configuration with a 
double cross, and C;**+ a configuration with two single crosses. Mean 
values are defined in the same way as (93), but v on the left-hand side 
corresponds to the number of vertices of the new configuration with crosses 
added. We then have 


Theorem II. 
SO ar eA ge rant a em ks; 's (100) 
This follows readily from the relation 
a a )” 5 + Bey) = €(SS;)*(S8x)9(S48))” - - eR 
= ¢(S,$;)%($j5;,)9(S,8))” - +» Sap") = ne #(S,8;,)9(S,S))”---Sp°), - (101) 


since S,,” is Brake and avila to X. pee the double cross is separated 
there is an extra (2s+1) in the trace, but this is cancelled by the extra 
vertex. As a consequence of this theorem double crosses, make no con- 
tribution to the susceptibility, since the coefficient of N in the total number 
of configurations C,,+ is zero. 


Theorem III. 
LG ee acar cema tine alts ee i! <2) (102) 
The reduction of such a diagram yields only odd powers of s, and hence the 
trace is zero. This is related to the vanishing of the first order terms in. 
the expansion of the free energy with respect to 1. 
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Theorem IV. If O,*+ has crosses at points i and j of a configuration 
and C,,, is the configuration obtained by replacing the crosses by the line 
ij (whether i and j are neighbouring lattice points or not) then 


{Crt} = HOR lon a ee eee 
This follows from the relation 
<(sjs i (S)S;,)9(S48))” - « « Szp8zq) = €(SS;)*($j8x)9(SZ81)” - - - SypSyq) 
= ((88;)*(SjS,)/(S,5))” - a an = §{(S8;)*(S;8;,)4(S,8))” - - - (SpSq) )- 
(104) 
The right-hand side is one of the permutations of {C/,,}) and each permu- 
tation can be treated similarly since s,,,s,, can be fitted in equally well at 


any point in the product... In counting the number of C;** we must 
remember a factor 2 arising from the cross-terms in the product 


Tl exp (2Ls,)). 


A typical example of C;+* and C;,,, is shown in fig. 16. 


Fig. 16 


Cas Cri 


A typical example of C7++and C,,,. Susceptibility configurations (with crosses) 
can be replaced by normal configurations. 


As a consequence of theorems III and IV configurations which fall into 
two parts on the removal of a single line can be ignored. For on replacing 
the line by crosses at each end it is clear that the trace of the configuration 
vanishes. 

Rushbrooke and Wood determined the zero field mean values of all 
configurations required for six terms of the susceptibility series (88) 
(which are polynomials in X.) As a result they derived six terms of the 
series (88) and seven terms of the series (89) for general s. The results 
are somewhat complicated, and we therefore do not reproduce them here. 


5.4, Low Temperature Lattice Constants 
5.4.1. Closed configurations 
We have discussed briefly the nature of low-temperature configurations 


in §5.1.2, and introduced the notation to be used. We have also seen 
in §5.2.1. that the classification of configurations into space types enables 
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low-temperature lattice constants to be determined. In fact for low- 
temperature expansions tightly bound configurations like P,, play a very 
significant part, and the method of space types becomes particularly 
important; the analytic methods employed for high-temperature lattice 
constants are correspondingly less useful. However, algebraic relations 
can be derived between high- and low-temperature lattice constants which 
enable the latter to be calculated. 

In a low-temperature configuration none of the vertices unconnected 
by lines may be nearest neighbours, whereas for high temperature con- 
figurations no such restriction exists. For close-packed configurations 
in which all pairs of vertices are connected the high-temperature and low- 
temperature lattice constants are the same; thus p,;=P3, peq=Peq- 
Close-packed configurations with five vertices (10 lines) do not exist on 
any of the lattices we have considered. A configuration of five vertices 


and nine lines exists on the hexagonal close-packed lattice P,,,; 


the high-temperature lattice constant would enter only at the tenth term 
of the susceptibility series. For the face-centred cubic lattice the maximum 
number of lines for a configuration with five vertices is eight. There are 
two such configurations P,, and P,, and for both of these we have 
Psq=Psq Psu=lPsy,; however, the latter was not tabulated among the 
high-temperature lattice constants since it enters the series (56) and (66) 
only at a much later stage. 

Let us now consider the relation between p;, and P;,. If we examine 
all the p,,, we find that we include all the P;,, when the unconnected opposite 
vertices are not nearest neighbours, and all the P,, when they are nearest 
neighbours; each P,, will occur six times since there are six independent 
Psq O11. a tetrahedron. Hence we have 


ore Oe tae, (105) 


Similarly, if we examine all the p, we shall include P, once, P;, once, and 
P,q three times. Hence 


aie Roan a ee... - = (106) 


In Appendix IV we have tabulated low-temperature lattice constants 
of closed configurations arising from tive spins or less for the lattices 
considered previously. For larger numbers of spins the lattices are best 
treated individually. Thus for the simple quadratic lattice constants 
have been evaluated for configurations of up to 10 spins, and for the 
simple cubic and body-centred cubic lattices for configurations of up 
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to eight spins; however, for the face-centred cubic lattice the number of 
configurations which enter increases rapidly, and is quite large even for 
six spins. 


5.4.2. Open and separated configurations 

The method of evaluating open configurations is closely similar to that 
described in §5.2.6, but when we attach a single bond, we must take 
account of the nearest neighbours of its freeend. Thus if we attach a single 
bond to P, we can form 


fat) 


as in § 5.2.6, but also P,,,, Pep, Pz, Psy. (However, not P,, since opposite 
vertices of P, cannot be nearest neighbours.) Hence instead of (34) we 
have 


4(a — 1)Py= Boi) + 6P,,+2P,.,+ 6P,,+4Pe,- a (107) 


Similarly low-temperature lattice constants for reducible configurations 
(e.g. P5,, P,.) can be calculated directly by a method analogous to § 5.2.4. 
For separated configurations we follow the procedure of §5.2.7. Thus 
by putting together a single spin and a triangle we derive the relation 


N?P,;=N[P3,@]+N [ TANS ] +3NP,+2NP,,+4NPea, 


and hence it is easy to show that 
[Ps, @]=N*P,—30P,+2P;,+8Peg - - - - (£108) 


We shall see shortly that by means of certain algebraic transformations 
all configurations which are not multiply-connected can be eliminated in 
deriving low-temperature series for the simple Ising model. Nevertheless 
this procedure introduces complications in the actual derivation of terms 
of the series, and other configurational problems on crystals involve lattice 
constants of all configurations. We have therefore reproduced in Appendix 
IV the lattice constants corresponding to all configurations consisting of 
up to five overturned spins. 


5.4.3. Low-temperature series expansion for the Ising model 


We have already discussed the form of the low temperature series 
expansion in § 3.6.1 as follows ($3 eqn. (143)): 


A(u, 2) = 1+ fiz) + w2fol2) +... +ulfilz)+ -.. 


We also derived the first few terms for the simple quadratic lattice. These 
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terms can be written in general form in terms of low-temperature lattice 
constants as follows: 


Fil) = 22, 
fale) =| 4 ga A 


2 2 
Js(2) = 22% {Pet a -8P, ga | <= 1) 


Z ) 


2 
+8P, fe (5 -P,)| . 


fyle)=242 {Poe + Pegi! + [Py + 3P9(q— 2) —4P5q—12P yg) 
tf ere — P,(12q—21)—4P,+ 6P,,+28Peq |2-8 
=e E (— 169? + 25q — 10) + 3P3(6q¢—7)+6P,—4P,,— 27P gq |e 
+ E (8q? —‘7q + 2) —3P3(4q—3)—4P,+ P54+ 12P ya fe? 


+| — fi (16g? —39+2)+ Pyl8q—1)+ Py 2Poa \. (109) 


Here the final term in the polynomial /,(z) corresponding to 24 arises from 
t separated spins, the term before this to ¢ spins with one bond, and so on, 
the first term corresponding to the most tightly bound configuration with 
the maximum number of bonds. 

It was shown by Domb (1949) that only a portion of the coefficients in 
the polynomials f(z) need be determined, and the remainder can auto- 
matically be calculated by algebraic manipulation. We shall illustrate 
the method by considering the simple quadratic lattice, but the results 
apply equally to any other of the lattices considered; if we start by 
determining most tightly bound configurations we will find that only 
configurations with ¢ spins and ¢ or more bonds need be considered. For 
the simple quadratic lattice we have 


A=14 pu + p?(2u3 — 2u4) + u3(6u4t— 140° + 8u8) + w4(ut+ 180° —77u8 
+98u?—40u8)+... (w=2?). cant Se ent (110) 


This is best regarded as a low-density expansion in p valid in the neighbour- 
hood of DE (§3 fig. 38). Ifwe put w= 1 in (110) we find that the A=1+ yp, 
and.this is the result we should expect in the limit of infinite temperatures 
($3, eqn. (156)). Hence f,(1)=0 for all t>1 and we have one relation 
between the coefficients of f,(z). 

We can now proceed to find an expansion nearw=1. Putting (1—u)=¢, 
and expanding A as a power series in ¢, the coefficient of ¢ is 


— Qu + 2yu27—2yF+2u4—... 0. ew (ETT) 
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and we may assume that the sum of this is —2y/(1+). Similarly the 
coefficient of ¢? is 

pe— Gu? + 16u2—3lpA+ ..., 2. 2 LES} 
‘and in this series second differences of coefficients are constant; this 
indicates that the sum is (~—3yu?+*)/(1+ )%. We are thus tempted to 
postulate a general expansion of the form 

S Polvo 

A=(1+p)+ at pra (en) ae (113) 

where the ¢,(j) satisfy the relation 
tele) =p (y). se ee (AI) 

In fact the expansion (113) is a consequence of §3 eqn. (156), when we 
remember that ¢,(7) is a polynomial in 7? whose maximum power is 7””. 

Thus we see that from (111) and (112) we can deduce the values of the 
first two rai at z=1, f(1), f@(1) for all ; In general from a 
knowledge of f,(u). . .f,_,(w) we can ee aA ), f(1), . ..f(1), and hence 
t relations efesrt the coefficients of the polynomial f,(u). 

If we consider the expansion of In A we have seen by analogy with 
the Mayer theory (§ 4.5.1) that separated configurations can be eliminated, 
and if we write as in § 3, eqn. (146): 

In A(p, 2) = 1gi(2) + 7 go(z) + ...w'g,(z) + «- 
the g,(z) will all depend linearly on the lattice constants (although the 
terms of f,(z) given in (109) depend linearly on the lattice constants, f,(z) 
would bring in P,? from separated triangles and similarly for higher terms). 
We may write the first four terms as follows: 


nO) =2| $23 G+) 1) ] 
qq—- 


fai rhe — SP, J++ ( (3P,—@2)z2 


+] 
abe ane 


SS 
— 
A 

I] 

xX 
wo 
Q 

——— 

y 


ENE 0 Peat 12+ Pez 194+ [P+ 3P;(q—2) —4P,, —12P.q]z 
q—1)(4q-—5 
+ |“ us —3P,(4q—6)—4P,+6P,,+28P,q |2-° 
8q?— 10q+3 
+| - ct eee + 18P,(q—1)+6P,—4P,,—27P,4 le-# 
4q?-—q+l1 
an eee — 6P3(2q—1)—4P,+ Py +12P¢5q |2-2 
1 8q°+ 6¢+7 
+|- ;- deen +3qP3+P, -2P xa ok (IBY 
However it will be seen that the general BERN is rather complicated. 
If we now introduce the variable f=(z-*—1) as in the Mayer theory 


(§ 4.5.1) we might expect an appreciable ene ee We shall express. 
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the result in terms of high temperature lattice constants since these arise 
naturally in the Mayer development. We find that 


9x(2) = 24, 

Jo(2) eel: ’ 

ale) =2| pg pee Ag) pe a+], 

Ja(z) = 244 {puah *+ Dad? + Paf*— 3psf? + eae ee! ps 


+| Mrs Ps |fP— 2q(q—1)f*— aq f?+ +taf—3}. 


We have written the contributions in such a manner that they can readily 
be identified with appropriate configurations. Thus for g,(z) all connected 
configurations with four vertices or less enter as follows: 


Pea Psa Pa Ps YS ; ENE ? EX > Up ? 


Proceeding with the Mayer theory and expressing In A as a power series in 


(116) 


0 
oer hers (In A) 
in the form I 


In A(ay, 2) = oyhy(2) + 047ho(z) + ...0,2,,(2) +... | 
we find that 
hy(z)=1 
ha(e)= — 8; = 4-5 f 
a (117) 
293, 492" 
h3(z) = — 3 + 2 =$+¢f?—2psf? 
91 91 
3g 18939, _ 20g.° _ Pe iareseesee 
ASS ery tea ea tare Baer EA) 


+ 9p3(2f4 +f *)— 39,4 f4— 3p5qf > — 3peaf °. 

A comparison between (117) and (109) shows strikingly the great 
simplification produced by the changes of variable in accordance with the 
Mayer theory; in (117) only lattice constants of multiply-connected 
configurations enter. Comparing (117) with §4 (82) we see that 


Iysal@)=— Bale) - AA Clee Hoan (LTS) 


Several of the terms contributing to h,,(z) can be derived from the results 
of Rushbrooke and Scoins (1955). Thus from § 4 (87) it is clear that the 
term independent of lattice constants is equal to 1/u. Similarly expanding 
§4 (90) in powers of «, we find that 

4H (04,2) = — doef t oy 3f?— oy GfP+of)t+-.. 6 « (119) 
which checks with (117). Terms in p; and p, can be derived in the same 
manner from §4 (91). This method cannot be applied to lattice constants. 
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of non-homogeneous configurations, although it seems likely that the 
more general formulae of Yvon $4 (105), (107) and (111) could be used. 
Alternatively, direct calculation can be undertaken from § 4 (83) starting 
with the given configuration having one system at each vertex, and then 
considering the possibility of multiple occupation of the vertices. 


5.4.4. Special methods for two-dimensional lattices 


All the methods so far described can be used to derive irreducible con- 
stants on two-dimensional lattices but it is usually easier to exploit the 
dual-net property mentioned in § 3.4.2. 

The method is best illustrated by an example. The simple quadratic 
lattice is self-dual and every connected Jow-temperature configuration of 
t spins and r bonds corresponds to an arrangement of / =qt—2r bonds 
with a connected interior on the dual lattice. For small values of / all 
such arrangements must be polygons, and for /=10 the only possible 
partitions are t=4, r=3; t=5, r=5; t=6, r=7 and from Appendix IV 
we see that these correspond to the low-temperature configurations in 
fig. 17a. We thus deduce at once that there are 18 


oY =14N Pale pe OM OO. 5 


(ii) (ili) (iv) 


decagons of area 4, 8N decagons of area 5 and 2N decagons of area 6 and 
that pj,=28. It is not necessary to specify the individual space types 
and the low-temperature configurations provide a systematic way of 
enumerating polygons on a plane lattice. For larger values of / not all 
structures with a connected interior are simple polygons but the exceptional 
cases are not difficult to allow for (e.g. fig. 17) and the method has been 
used to derive p14, i, and p,, for the simple quadratic lattice. Ina similar 
way polygons on the triangular lattice may be obtained by enumerating 
the corresponding duals on the honeycomb lattice. 

The dual lattice also provides a convenient way of evaluating other 
high-temperature lattice constants. For example, the number of polygons 
with a single bridge can be determined by considering the possibility of 
crossing single bonds in the low-temperature configurations of 
fig. 17a; it is easy to see that there are 42N from (i), 12N from (ii) 
8N from (iii) and none from (iv), and these can be sorted into various 
categories. Similarly polygons with double bridges arise from the 
possibility of crossing certain pairs of bonds in the low-temperature 
configurations. These can only occur with polygons of types (iii) and (iv) 
there being six possible double bridges in each case. 


5.4.5, Low-temperature antiferromagnetic expansions 


As remarked in §3.6.1, low-temperature antiferromagnetic expansions 
can be derived for lattices which can be decomposed into two equivalent 
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sub-lattices. The enumeration of disordered states is more complicated 
since the contribution of any state of overturned spins now depends on 
the manner in which the spins are distributed over the two sub-lattices. 
Thus for example on the simple quadratic lattice there are 4N(N —5) 
ways of choosing two separated overturned spins and in the case of a 
ferromagnet the contribution of these would be 4N(N —5)u?28. For the 


Fig. 17 


(0) 
(a) Decagons on the simple quadratic lattice classified by means of low- 
; temperature configurations on the dual lattice; (6) A connected low- 
temperature configuration whose dual is not a simple polygon. 
— — — — Original lattice. —— Dual lattice. 


antiferromagnet we must distinguish between the }N(jN—1) ways 
of choosing two spins on the same sub-lattice and the 3N(3N —4) ways 
of choosing two spins on different sub-lattices. The contribution of 
the first will be 4N(4N — 1)u2y8 for the sub-lattice whose spins are aligned 
with the applied field and }N(4N —1)u-*y® for the sub-lattice whose spins 


P.M.S. ZA 
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are aligned against the applied field. The contribution of the second type 
will be }N(4N —4)y8. The sorting out of the low temperature configura- 
tions given in Appendix IV into sub-lattices is tedious but straightforward 
and the series given by Brooks and Domb (1951) for the simple quadratic 
lattice is readily derived. The corresponding series for the simple cubic 
lattice has been given by Wakefield (1951). 

One might again expect a substantial simplification to arise from the 
application of the Mayer theory, and an approach of this kind has been 
initiated by Temperley (1959). Alternatively it seems very likely that 
the method of Yvon (§ 4.5.2.) could be usefully applied to this problem in 
a manner similar to § 5.4.3. 


Notation for § 5 


Configurations classified by p points and / lines, c=q—1, g coordination 
number of lattice. 


Typical high-temperature lattice constant of / lines p,,. 


Simple closed polygons p,, simple chains ¢, (If we wish to record the 
number of vertices v of the configuration we write p,,”). 


Typical low-temperature lattice constant of / lines P,,. 
Lattice constants of separated configurations [p)., Pmy| [Piz Pmy]- 
r, number of returns to the origin in a random walk of / steps. 


q, number of returns to the origin in a no-reversal random walks of / 
steps (excluding improper closures). 


A,2(l) path factor of configuration p,., in random walk of / steps. 
P(x, y,z) generating function for random walks on a lattice. 

(x,y,z) generating function for no-reversal random walks on a lattice. 
(a, 6,c) numbers of walks arriving at (a,b,c) after J steps. 

Q,(a, b,c) number of no-reversal walks arriving at (a,b,c) after 1 steps. 


1—z : ; : 
w=tanh K = peer variable for high-temperature expansions. 


6,,d, coefficients in susceptibility expansion for Ising model (s=4) 
[eqn. (63), (65) ]. 

C, number of occurrences of a general configuration (with multiple 
bonds) on the lattice. 


{Cj} mean value of this configuration. 
Z,* partition function for Heisenberg model. 
t= sty? + (8— 1)eu-e D+... (—8 + 1 apne  fe ( me 8S)", Tg=balto. 


Jj, coefficients in susceptibility expansion for Heisenberg model 
(s=)leqns. (88), (89)]. 


SAz), 9,(2), ¢$,(7) as in § 3. : 
u=2,C=1—u, o,() defined by (113), h,(z) by (117), H,,(04; 2) by § 4 (86). 
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APPENDIX I 


ANALYSIS OF RETURNS TO THE ORIGIN IN A RANDOM WALK 


(a) Simple quadratic lattice 


fe =4 22 [— 
nose =2 [SJ +4[o]+0(0] 
400 =2 [—] +2[A]+ 0[-~] +12 [>] 
+aa(O]+12 [T)}+2[O] 
4900= 2 [—] +20 ~] +32 : 2 lhe 
+16 [~<]+ 48] ]+264[ + n2[T ]+I6 [CE | 
rio (a) +6l sic p sl hk} 


=~ 
@ 
~—. 
= 
NO 
ie) 
— 


| 
ave (CO) sfc} fe 3} 


(2) 


(b) Simple cubic lattice 


= AeA 
-2(- =|+4[~}+0[0 D1 | 
r, = 1860= fie -] +12 cone: fw] +12 [> ]+48 [0] 
Ha +12[ <>] 
aia sJexfv]+ 7 Psi 


(363) 


mite 

+16 [XS ~]eel x x<}+269O}+e[ ET of a 
+ 16 [4] +16 iso SN slcb hos lk 

[CO 


cts Ee ob 


+48 
(18) (207) 
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(c) Combinatorial factors (which represent the total number of returns to 
the origin having the given configurational pattern, these returns being 
summed over all vertices of the configuration) 

atte 1)!(7g+73—1)! 


ip S 
lll p ea pptritret+rs—1)'\(p+si+ $24 83— 1)! 


= = oy ate) rt %3—1)! 


[= 2(ry+7ret+1rs) +4. 
APPENDIX II 


r, AND q, FOR SIMPLE LarTicEs 
Generating functions $(x, y, 2). 
in [A(x, y, 2) } 
Simple quadratic 
o=tar+u+y"+y), 
r= (L/L) Nt (even). 


(r, is the coefficient independent of x, y, z 


Triangular 
p=t¢(ut+a+y t+y+ay+a-ly), 
= bs uh . gern. rg eh ee 
3, t 8 lt! 5=5 [(t+q)/2!}? q!(s—q)!’ 


$,t=0,1...1; (28+t)=1; (t+ ) even. 
Simple cubic 
o=$(e +u+y1+y+z21+42), 
r= ~ acca r,s,t=0,1...v,r+s+t=v, l= 2p (even). 
Body-centred cubic 
P= EU +2)(y*+y)(2-1 +2), 
r= (1!)3/((31) 1} (Leven). 
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Face-centred cubic 
=Sl(y* +y)(e1 +2) + (2 1+2)(a41 +a) + (x ++a\(y41+y)], 


Le ay ne (ocean its (eid)! 
T= Ree aS) ce ee ere (rd) Um § 
 Arlsltl (Gers) Plae+o) Plae+oye Tet h 
r; 8,t=0,1...1, 7, 8, ¢ all odd or all even. 
Honeycomb} 


ig ‘ 1/2)1)2 
@=(ltat+y)'a1 +a+y)t2, n= 2, oy ne 
leven r+s+t=1/2. 
Diamond+ 
$=[(wtat+elyty)PP[(et+at+ety+y4)p? 
S a2" > (28)! (1 — 2a)! 
r= [(1/2)!? 2, )*[U2—s) 


l even. 


Hexagonal close-packed. See §5.2.8. %, q, as for face-centred cubic. 


NUMERICAL VALUES OF 1, g, FOR SIMPLE LATTICES 


Loose-packed lattices 


Honeycomb S.Q. Diamond S.C. B.C.C. 
ry 3 4 4 6 8 
"4 15 36 28 90 216 
1 93 400 256 1860 8000 
Ts 639 4900 2716 44730 343000 
Tio 4653 63504 31504 1172556 16003008 
Ty. 35169 853776 387136 32496156 788889024 
Io 0 0 0 0 0 
da 0 8 0 24 96 
%s 6 24 24 264 1776 
Ys 0 216 48 4584 67584 
dio 30 1520 480 88560 2455200 
Os 66 12080 3432 1750800 94735920 
Close-packed lattices 
Triangular dete py Triangular Pa as 
fo 6 12 Gaie0 0 
7 oe -12 48 Ge send? 48 
ts 90 540 qa 24 264 
rs 360 4320 V5 60 1680 
Le 2040 42240 I 444 17016 
fo 10080 403200 Iq 1680 151872 
r, 54810 4038300 Is 7464 1401720 
To 290640 40958400 Io 35040 13238400 


+ These non-Bravais lattices require special treatment. The generating 
functions quoted will give 7 correctly. 
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APPENDIX III 
TABLE A 


Lattice constants of closed configurations 
(} denotes irreducible configuration) 


Constant Triangular pp ant CP(6) CP(7) 

o=q—1 5 11 4 5 
Pst /\ 2 8 10/3 5 
Pat | 3 33 15/2 15 
pst 4), 6 168 12 36 
Psa iY] 3 36 15 30 
Det <a 15 970 10 60 
Peat 0 36 10 30 
Per S e 384 60 180 
Pec Pea 204 15 45 
Peat xX] 0 ” 5/2 5 
Prt CD, 42 6168 0 360/7 
Pra ES 12 966 30 180 
Day <<) 30 2400 60 360 
Pr hb 6 192 60 180 


Pre 


Past 


Per 
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Constant 


SOOVOSSORGI TS. 


TABLE A (continued) 


Triangular 


30 


123 


96 


60 


24 


12 


66 


Face-centred 


cubic 


2196 


2040 


600 


48 


42069 


16464 


11616 


3888 


3132 


2016 


1056 


12096 


CP(6) 


15 


30 


30 


10 


30 


120 


60 


347 


CP(7) 


90 


180 


180 


30 


90 


360 


360 


180 


180 


720 


360 


180 
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Table A (continued) 


Constant Triangular Bs ese CP(6) CP(7) 
Di aS 24 4968 0 90 
Pai ae 144 22176 0 180 
Pex or: 186 22632 0 90 
Psi ae 12 1032 30 180 
Pam Ne 18 1176 30 180 
ry oy 0 528 60 360 
Peg <> 0 24 30 90 
Dart <@> 0 9 5/2 15 
Dae q> 0 384 30 180 
Pot e@ 380 301376 0 0 
Pon OM 180  ~° 79920 0 0 
Por S 102 58032 0 0 
Poe SN 48 6432 0 90 
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TABLE A (continued) 


Face-centred 


Constant Triangular 
cubic 


CP(6) CP(7) 


2 744 0 15 


&»> 0 2352 0 180 
Pos <b ; 0 1848 0 180 


2 160 20 120 


Pont 0 192 0 60 
Content eae, Se Le 
o=q-1 3 5 7 3 

Pat Ee 1 3 12 9/2 

Pet Cs 2 22 148 12 

Peat 0 0 12 6 

Dra <> 2 18 192 36 

pet eS 7 207 2736 9 

Pest <> 0 24 480 36 

Den — 2 30 618 18 

Part <e>> 0 0 y de 
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TABLE A (continued) 


Cantata) Se 
Patt Ke» 0 0 24 18 
Dos a“ 8 192 4632 18 
Dok me) 12 324 6576 72 
Det a 0 20 436 12 
Pom Vane 0 0 456 72 
Prot @) 28 2412 61896 0 
Proa Co 8 576 19584 0 
Prov CS) 0 12 304 12 
Pret <> 0 0 120 0 
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TABLE B 


Lattice constants of open configurations 


ey 


ys | 


31 


3:2 


4-1 


4-2 


4-3 


4:4 


5:2 


5:3 


5-4 


5:5 


Simple 
quadratic 


XBR AASRXARY> > 


bo 


18 


36 


50 


Triangular 


15 


20 


69 


24 


15 


252 


309 


42 


90 


36 


108 


Simple 
cubic 


15 


20 


75 


15 


300 


363 


48 


Body- 
centred. 
cubic 


28 


56 


196 


70 


1176 


1324 


288 


56 


351 


Face- 
centred 
cubic 


66 


220 


702 


240 


495 


6780 


7350 


1248 


2328 


1080 


2496 


792 


352 


or 
or) 


5:7 


5:9 


or 
= 
oS 


FO OAPEAQSSSXX 


6-1 


6-3 


6-4 


6-5 


6-6 


6-7 


6:8 


Simple 
quadratic 
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TABLE B (continued) 


— 
bo 


18 


92 


100 


142 


16 


24 


0 


Triangular 


228 


1104 


1020 


1365 


96 


24 


18 


75 


132 


54 


186 


104 


Simple 
cubic 


300 


300 


1404 


1452 


1767 


96 


240 


Body- 
centred 
cubic 


1960 


1764 


7656 


7944 


8980 


828 


1728 


720 


1944 


Face- 
centred 
cubic 


19620 
15786 
69924 
68508 
76266 
7632 
696 
648 


5808 


11400 


12744 


7296 


6-9 


6-10 


6-11 


6-12 


6-13 


6-14 


6-15 


6-16 


6-17 


6-18 


6-19 


6-20 


6-21 
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quadratic 


Meee 


TABLE B (continued) 


Simple 


i) 


12 


38 


184 


520 


268 


Triangular 


24 


792 


192 


288 


480 


102 


330 


966 


3240 


8796 


4836 


Simple 


cubic 


150 


600 


1308 


5616 


13656 


6876 


Body- 
centred. 
cubic 


28 


1960 


5880 


10656 


45936 


103872 


52152 


353 


Face- 
centred 
cubic 


20304 


2880 


47760 


11784. 


21648 


25872 


924 


37800 


87948 


163926 


627912 


1400328 


719748 
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TABLE B (continued) 


Body- Face- 


Simple qyiangular ™P!€ centred _centred 
quadratic ng cubic Sates. cubic 


84 1312 2260 17240 226780 


50 1242 2178 19860 286638 


6-23 owe 28 984 1356 12280 199476 
YM OAT A 


390 5973 8463 60028 786858 


Lattice constants of open configurations (simple chains) 


The values of ¢, the total number of non-intersecting chains of / lines, for 
various lattices. 


cusdiatig: © Toanealen, ee 
Cy 2 3 mT. 4 6 
Co 6 15 15 28 66 
Cs 18 69 75 196 702 
C4 50 309 363 1324 7350 
Cs 142 1365 1767 8980 76266 
mp 390 5973 8463 60028 786858 
Cy 1086 25941 40695 402412 8086074 
Cg 2958 112065 193983 2675860 82848522 
Cy 8134 482067 926943 17826340 = 
Cio «22050 2066583 4 en a 
Cy 60146 na = a ~=a 


C1. «162466 a = es 2h 
Cs 440750 = pm a ve 
Cig «21187222 ie e Ra a 
C45 3208298 hed << Ta ce 


Cy_e «8622666 a we _ fr. 


Separated configurations. 
of NV) where / is the number of lines) 
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TABLE C 


(To avoid fractions we tabulate 1! (coefficient 


etre Triangular 
/ fe ee = 33 
UG) ietor? its 
/ i iy, ee +864 
iN if 0 —720 
J i —1248 —9504 
NG vy, — 5520 — 33336 
aX i i 413872  +85176 
/ fi i / 6252 —36846 
Vex IN 2976  —19260 
LJ y, —1440 —6840 
[V/ / 0 — 56160 
x< / —1920 42480 

— 67200  —732960 


INL | 


Simple 
cubic 


—33 


— 1440 


— 10080 


—37512 
+ 90072 


— 38286 
— 20556 


— 7200 


— 46800 


— 924480 


355 


Body- Face- 
centred centred 
cubic cubic 
— 60 —138 
— 3696 — 13320 
+2192 +7752 
0 — 6336 
— 38976 — 231840 
— 135264 — 743472 
+ 319392 + 1766736 
— 132456 — 720756 
— 74688 — 421920 
— 40320 — 169920 
0 — 1249920 
— 302400 — 3196800 
—5011200 —44144640 
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TABLE C (continued) 


: Simpl Body- Face- 
sua Triangular ee centred centred 
quadratic cubic habne cahak 


/\/\ / 94080  —913680  —1121040 —5650560 —48113280 


x / if 0 453280 0 0 + 1036800 


ENS / / 468640 +812160 +908640 +4864320 +44160480 
UG / if 4302160 +2883240 +3372840 +16829280 +142173360 


JX if i ib _ 677280 —6267840 —7201440 —35373120 —300739680 
ENED 4235488 +2197872 +2359152 411313216 494744224 


ES JAN A 4327840 +3329280 +3705480 +18639360 +161395200 


iS ifs 0 — 24480 0 0 — 498240 
JAWS KES 424960 —317520 —354240 —1752960 —18123840 


TONE /\_ +111840 1119600 = 1324080 —6106560  — 58263840 


Separated configurations required for the first nine terms of the Ising 
no-field series 
[Ps Ps] = Nps"/2 —(Pect+ Pat $Ps) 
[Ps Pal=NPsPa— (Pret 3P7q + Pov + 2P5a) 
(pa, Pal = Npy?/2 — (Pen + 3Pgr + Pas + Pra t+ Pret 3P6a + 3Peat $a) 
[Ps P5]=NPsP5— (Pag + 2Pgs + Pry + Pre + Poo) 
(Ps, P6l=NPsP6— (Poa t+ 2Poet+ Por t+ Pog + Pea t+ Pse + Pro) 
[Pa P5]= NpsPs— (Pon + Poc + 2Poz + 8P9qg + 3P gn + Psy t+ Paet 2Pap 
+ 2pgq t+ 2P7¢+ 2P7n+ Peon) 
[P3> Pool = NPsP6e— (2Pge + BPga+ BPgqg + 2Pg + 2Pam + 2Pgq + Pre + 2ec) 
[P3, Ps, Ps |= N*p39/6 —-N(PecPs + PsaPs t+ P3"|2) + (Poet 2Poat 2Pog + 2g 
+ Pam + 2Pgaq + 2Pre + 2P 7g + 2P Ge + SPea+ 2P5a + Ps/3) 
For loose-packed lattices we quote an additional term 
[Pa P6]=NPaP6— (Prod + 3Pr105 + 8Pr00 + Po + 2Pom + 2Pge + 2D + 22a): 
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APPENDIX IV 
Table of Low-temperature Configurations (Evaluated for N=1) 


Simple Simple Body- Face- 
Configuration PS j centred ‘Triangular centred 
quadratic cubic : 
cubic cubic 
a +2 +3 +4 43 +6 


8 +15 +28 +9 +42 
ff ry —14 —33 — 60 —27 —114 
eae +8 +18 +32 +16 + 64 


y 
g 


vo HN x 


c 


0 0 0 +3 +24 


Pe +1 +3 +12 0 bes 
+120 
+14 +63 +148 +27 +282 

+4 +20 +56 +2 +44 


2B 
P.M.S. 


358 C. Domb on 
APPENDIX IV.—Table continued 


Simple Simple Body- Hace- 
Configuration apes P centred Triangular centred 
quadratic cubic : ; 
cubic cubic 
e . 
0 0 0 —22 —192 


f if ~21 =e —194 —57 —513 
UN ° — 56 — 228 ~568 —108 — 1080 


Uf ee +98 +360 +868 +249 +2166 


— 40 —137 —322 — 106 — 856 


av 
e 
Es 
e 
t e 
a 
i=) 
Oo 
<= 
rm) 
th 
bo 
ns 


é 


@ 
fons) 
Oo 
Oo 
oO 
| 
Or 
rx 


“ 
Kl > 
° 
° 
Pa 
bo 
° 
) 


+24 


0 0 0 +6 +168 


oe 
° 
° 
° 
° 


Configuration 
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APPENDIX IV.—Table continued 


Simple 
quadratic 


+ 
oO 


| 
—_ 
—_ 


PbUNRREUYS 


INS +34 
NE +20 
<i fl 


Simple 
cubic 


+48 


—57 


+ 204 


+15 


Body- 
centred. 
cubic 


+216 


— 288 


Triangular 


+12 


+3 


—39 


— 90 


— 168 


+81 


359 


Face- 
centred 
cubic 


+240 


qt 72 


+48 


+792 


+576 


+96 


— 696 


—87 


— 1680 


— 3672 


+1902 


+828 


+9 
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APPENDIX IV.—Table continued. 


gimol Simpl Body- Face- 
; imple a tred ‘Triangular centred 
Configuration quadratic cubic sari é eaibag 


+222 +4080 


© 
© 
© 
o 
oO 


> @ — 44 — 372 — 1360 — 30 — 1404 

ING — 166 — 1227 — 3860 —405 — 9162 
e 

/ —160 —1077 — 3536 — 441 — 9480 


+349 +2130 + 6834 +1203 + 22773 


ute 


+450 + 2925 +9600 +1143 + 23910 


2 


—706 — 4041 — 12852 — 2370 — 42366 


ns 


eoece +225 +1185 +3672 +789 +12957 
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CALCULATORS 


RECIPROCALS OF THE 
INTEGERS FROM 
1000 TO 9999 


By T. H. Redding 
M.Sc.(Lond.), A.M.I.Mech.E., A.F.R.Ae.S. 


FOR USERS OF 


T 


as 


ye 


Arranged for use with mechanical 

calculating machines to facilitate the 

evaluation of quotients and com- 
pound fractions 


“= ray : 


z . oo 


“ 


With an Appendix on mechanical barrel- 
setting calculators 


These tables list the 10° multiples of the 
reciprocals of the integers 1000 to 9999 
and are displayed as 1000 entries on each 
of nine double-page spreads which are 
thumb-indexed (1000... , 2000... , 3000 
4 ..+)++. etc) to facilitate rapid manipulation 
sh one hand only. The entries are direct reading and are correct to the nearest integer 
) the third place of decimals. 

Although suitable for general use, the tables have been prepared and arranged in a 
Janner particularly suitable for use in conjunction with mechanical computing machines 
f the ‘‘ barrel-setting ’’ type in the evaluation of quotients and compound fractions. 
jally the tables find application in effecting the evaluation of quotients on adding and 
ing machines when these incorporate a mechanism for automatically evaluating simple 
ducts. 

In the first instance it was anticipated that an Appendix dealing with the advantages to 
gained by the use of the reciprocal function in conjunction with such machines—and 
ticularly with barrel-setting machines—would suffice to ensure the best use being made 
e tables. The present 22-page Appendix, however, goes further than this, since the 
ssion of the tables is prefaced by a description of barrel-setting machines and their 
lest) method of operation in the evaluation of products, quotients and compound 
tions. Thus, the description of each calculating procedure is followed by a numerical 
imple with accompanying diagrams to illustrate the keyboard-displays at the beginning, 
end, and at intermediate stages in the calculation. The machine-evaluation of series 
Aprising the sum of a number of product or quotient-terms is also discussed. Generally it is 
ed that the Appendix will serve as an introduction to mechanical methods of computation 
that it will materially assist prospective purchasers in the choice of a machine. 

A valuable guide to all users of Computers. Bound in blue publishers case, lettered on 
he front, with an 18 page table index. Printed on good quality paper for constant use. 
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